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1. INTRODUCTION 

If we are given the age distribution of a population on a certain date, we may require to 
know the age distribution of the survivors and descendants of the original population at 
successive intervals of time, supposing that these individuals are subject to some given age- 
specific rates of fertility and mortality. In order to simplify the problem as much as possible, 
it will be assumed that the age-specific rates remain constant over a period of time, and the 
female population alone will be considered. The initial age distribution may be entirely 
arbitrary; thus, for instance, it might consist of a group of females confined to only one of 
the age classes 

The method of computing the female population in one unit’s time, given any arbitrary 
age distribution at time ¢, may be expressed in the form of m+ 1 linear equations, where 
m to + 1 is the last age group considered in the complete life table distribution, and when 
the same unit of age is adopted as that of time. If 


n,, = the number of females alive in the age group x to x + 1 at time ft, 

P, = the probability that a female aged x to x +1 at time ¢ will be alive in the age group 
x+1tox+2 at timet+1, 

F,, = the number of daughters born in the interval ¢ to ¢+ 1 per female alive aged z tox+1 
at time t, who will be alive in the age group 0-1 at time ¢+ 1, 


then, working from an origin of time, the age distribution at the end of one unit’s interval 
will be given by 


m 

PE = mou 
Py Noo mm Waa 
Pin = Nq) 
Py Na = Ns) 


Pa-1%m-1,0 = Mma 
Biometrika 33 
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or, employing matrix notation, Mn, = n,, where ny and n, are column vectors giving the 
age distribution at ¢ = 0 and 1 respectively, and the matrix 


; - F, Me wes : _ ta << Fn Fn] 
: eh 
os 
2 gD Re er ere rn areie ner ros 18 2 SO Se. eae 0<P,<1; F209. 
Pry . 
P, 
| Fa = 








This matrix is square and consists of m+ 1 rows and m+ 1 columns. All the elements are 
zero, except those in the first row and in the subdiagonal immediately below the principal 
diagonal. The P, figures all lie between 0 and 1, while the F, figures are by definition neves- 
sarily positive quantities. Some of the latter, however, may be zero, their number and 
position depending on the reproductive biology of the species we happen to be considering 
in any particular case, and on the relative span of the pre- and post-reproductive ages. If 
F,, = 0, the matrix M is singular, since the determinant | M | = 0. 

Since Mn, = n,, and Mn, = M;,, = Np, etc., the age distribution at time ¢ may be found by 
pre-multiplying the column vector {m9 %49 Ngo --- Mmo}, 1-€. the age distribution at t = 0, by 
the matrix M‘. Moreover, it will be seen that with the help of the jth column of M‘ the age 
distribution and number of the survivors and descendants of the n;_, individuals, who were 
alive at t = 0, can readily be calculated. Thus, n,_, 5 times the sum of the elements in the jth 
column of M' gives the number of living individuals contributed to the total population at 
time ¢ by this particular age group. 


2. DERIVATION OF THE MATRIX ELEMENTS 


The basic data, from which the numerical elements of this matrix may be derived, are given 
usually in the form of a life table and a table of age specific fertility rates. To take the P, 
figures first; if at ¢ = 0 there are n,) females alive in the age group 2 to x + 1, the survivors of 
these will form the x + 1 to x +2 age group in one unit’s time, and thus P,n,») = nz,;. Then 
it is usually assumed (e.g. Charles, 1938, p. 79; Glass, 1940, p. 464) that 





a etl 
P, = a? 
zr 
z+1 
where L, = l, dx, 
z 


or the number alive in the age group x to x + 1 in the stationary or life table age distribution. 
This method of computing the survivors in one unit’s time would be exact if the distribution 
of those alive within a particular age group was the same as in the life-table distribution. 
The F, figures are more troublesome, and in the numerical example which will be given 
later they were obtained from the basic maternal frequency figures (m, = the number of 
live daughters born per unit of time to a female aged x to + 1) by an argument which ran 
as follows. Consider the n,. females alive at t = 0 in the age group x to x + 1, and let us sup- 


P. H. LESLIE 185 


pose that they are concentrated at the midpoint of the group, x + 4. During the interval of 
time 0—1 some of these individuals are dying off, and at ¢ = 1 the n,,,, survivors can be 
regarded as concentrated at the age x+14. Although these deaths are taking place con- 
tinuously, we may assume them all to occur around ¢ = 3, so that at this latter time the 
number of females alive in the age group we are considering changes abruptly from n,» to 
Nz+11 = PN. Then during the time interval 0—} these n,) females will have been exposed 
to the risk of bearing daughters, and the number of the latter they will have given birth to 
per female alive will be given by the maternal frequency figure for the ages r+ } to x+1. 
This figure may be obtained by interpolating in the integral curve of the m, values, and thus 
expressing the latter in } units of age throughout the reproductive span instead of in single 
units. The daughters born during the interval of time 0—} will be aged 3-1 at ¢ = 1, the 
number of them surviving at this time being determined approximately by multiplying the 


: 
l.dx according to the given life table. Similarly, 
t 


each of the Pn,» females during the interval of time }-1 give birth to m,,,_,,,, daughters, 
the survivors of which form part of the 0—} age group at ¢ = 1. The survivorship factor is 


appropriate m,,, figure by the factor 2| 


} 
in this case taken to be 2 Lda. 
0 


Combining these two steps together we obtain a series of F, figures, which may be defined 
as the number of daughters alive in the age group 0-1 att = 1 per female alive in the age group 


xtox+latt=0. Putting 
' , 
k,= 2 dz, k= 2| L.dz, 
0 + 


then F, = (kam, yg + hy Pz mys, 2414), 


m 


€ \" " — 
and > F229 = %1; 
z=0 


the total number of daughters alive aged 0-1 at ¢ = 1. 


3. NUMERICAL EXAMPLE 


In order to see whether the P, and F, figures obtained in this way from the basic data give 
a reasonably accurate estimate of the population in one unit’s time, a numerical example 
was worked out for an imaginary rodent population, the species chosen being the brown rat, 
Ratius norvegicus. Full details of the basic life table and fertility table which were used are 
given in an appendix, together with a short account of the genesis of these tables and the 
methods employed to estimate the rate of natural increase (7) and the stable age distribution. 
Compared with man, the fertility of this imaginary rat population was relatively very great; 
thus, the gross reproduction rate was 31-21 daughters and the net rate (Ry) 25-66, the life 
table used being a reasonably good one. The inherent rate of natural increase was estimated 
to be 0-44565 per head per month of 30 days, and the stable age distribution was so overladen 
with young that the proportion of females in the post-reproductive age groups was negligible. 
Some 74-45 % of the females were younger than 3 months, at which age breeding was 
assumed to commence. 

By definition the Malthusian age distribution is stable; that is to say, once a population 
subject to the given rates of fertility and mortality achieves this form of distribution, it 
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continues to increase e” times every unit of time and the proportions of the population alive 
in each group remain constant. Thus, in the present example, given 100,000 females dis- 
tributed as to age in the stable form at t = 0, the number alive in each age group in 1 month’s 
time can be immediately valculated by multiplying each element in the original distribution 
by 1-561505. This ‘true’ age distribution at t = 1 is compared in Table 1 with that obtained 
by operating on the original distribution with the P, and F, figures, which are given in 
Table 5 of the Appendix. 

The agreement between the true and estimated age distributions is remarkably close. It 
might be expected that the principal errors would occur in the early age groups, since the 











Table 1 
(1) | (2) (3) (4) 
(Units of Population Expected popula- | Population at t=1 
30 days) at t=0 tion at t=1 Estimated by 

Age group Stable age Col. 2 x operating on col. 2 
distribution 1-561505 with the matrix M 

0- 37,440 58,463 58,374 

1- 22,595 35,282 35,455 

2- 14,417 22,512 22,519 

3- 9,227 14,408 14,406 

4- 5,903 9,218 9,218 

5- 3,775 5,895 5,895 

6- 2,413 3,768 3,768 

7 1,542 2,408 2,407 

8- 984 1,537 1,537 

9- 627 979 980 

10- 399 623 623 

1]- 254 397 396 

12- 161 251 251 

13- 101 158 159 

14- 64 100 99 

15- 40 62 62 

16- 25 39 39 

17- 15 23 24 

18- 9 14 14 

19- 6 9 8 

20- 3 5 5 

Total 100,000 156,151 156,239 




















The span of the reproductive ages is from 3 to 21 months. 


P, figures are based on the stationary age distribution which is clearly very different from the 
stable form. However, as will be seen from Table 1, the biggest error from this cause is due 
to the first P, which overestimates the number alive in the 1-2 age group at ¢ = 1 by some 
0-5 %. The F, figures underestimate the number alive in the 0-1 group by 0-2 %, and the 
total population is overestimated by 0-06 %. On the whole these results are satisfactory 
and, judging from this example, it would seem that the matrix M operating on a given age 
distribution should give a reasonable estimate of the population in one unit’s time, provided 
that the unit of time and age chosen be not too coarse as compared with the life span of the 


species. The degree of cumulative error which is introduced by continued operation with the 
matrix will be considered later. 
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4. PROPERTIES OF THE BASIC MATRIX 
The matrix M is square and of order m + 1; it is not necessary, however, in what follows to 
consider this matrix as a whole. For, if z = kis the last age group within which reproduction 
occurs, F, is the last F,, figure which is not equal to zero. Then, if the matrix be partitioned 
symmetrically at this point, A. 
u-[5 ol: 


The submatrix A is square; B is of order (m—k) x (k+1); C again is square consisting of 
m — k rows and columns, the only numerical elements being in the subdiagonal immediately 
below the principal diagonal. The remaining submatrix is of order (k+1) x (m—k) and 
consists only of zero elements. Then in forming the series of matrices M*, M*, M‘, etc., 


=> ee int 


The submatrix C is, however, of such a type that C™—~* = 0, so that M'‘, t>m-—k, will have 
all its last m—k columns consisting of zero elements. This is merely an expression of the 
obvious fact that individuals alive in the post-reproductive ages contribute nothing to the 
population after they themselves are dead. It is the submatrix A which is principally of 
interest, and in the mathematical discussion which follows, attention is focused almost entirely 
on it and on age distributions confined to the prereproductive and reproductive age groups. 


The matrix A is of order (k + 1) x (k+ 1), where x = kis the last age group in which repro- 
duction occurs, and written in full, 


RRRAR ... Ra HI 


See eee eee eee eee eee eeeoeeeesesesseses 








This matrix is non-singular, since the determinant | A| = (—1)*+?(P,P.P,... P._.F,). 
There exists, therefore, a reciprocal matrix of the form 








eC P;* : : rm ‘ 
; pe : 
A- = . . . Py 
: ; ‘ : ae Py 

| Fg} —(Ph)' HR -(AA)'H -(AATA :. - (Rea) Re 
Thus, given an initial age distribution n,) (x = 0,1, 2,3, ...,k) at ¢ = 0, in addition to the 
forward series of operations Any, A2no, A®no, ..., etc., there is also a backward series A~'no, 
A-*n , A- no, ..., etc. There is, however, a fundamental difference between these; for, 


whereas the forward series can be carried on for as long as we like, given any initial age dis- 
tribution, the backward series can only be performed so long as n,, remains 2 0, since a 
negative number of individuals in an age group is meaningless. Apart from this limitation, 
it is possible to foresee that the reciprocal matrix might be of some use in the solution of 
certain types of problem. 
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5. TRANSFORMATION OF THE CO-ORDINATE SYSTEM 

Hitherto an age distribution n,, has been regarded as a matrix consisting of a single column 
of elements. For simplicity in notation, this column vector will now be termed the vector 
£ and different £’s will be distinguished by different subscripts (£,, €,, etc.). We may picture 
an age distribution as a vector having a certain magnitude and related to a definite direction 
in a vector space, the space of the £’s. The different age distributions which may arise in the 
case of any particular population will be assumed to be é’s all radiating from a common origin. 
The numerical elements of a € vector are thus taken to be the co-ordinates of a point in 
multi-dimensional space referred to a general Cartesian co-ordinate system, in which the 
reference-axes may make any angles with one another. At this point in the argument another 
type of vector will be introduced, which in matrix notation will be written as a row vector, 
and which will be termed the vector 7. There is an intimate relationship between this new 
type and the old, for, associated with each vector &,, there is a uniquely determined vector 
Nq, and vice versa. The inner or scalar product, 7, ,, is the square of the length of the vector 
£,. Either we may picture each of these vectors as associated with a different kind of vector 
space, the space of the £’s and the dual space of the »’s, which are not entirely disconnected 
but related in a special way; or, alternatively, we may regard them as two different kinds of 
vector associated with the same vector space. The relationship between 7 and ¢ is precisely 
the same as that between covariant and contravariant vectors in differential geometry. 

If we pass from our original co-ordinate system to a new frame of reference, and the 
variables y and & undergo the non-singular linear transformations, 

n= 6H, -€=H"y, |H|+0, 

it can be seen that since the variables are contragredient, 7§ = dy, so that the square of the 
length of a vector remains invariant. Moreover, since the result of operating on a vector £, 
with the matrix A is, in general, another vector £,, where €, and &, are both referred to the 


original co-ordinate system, it follows that in the new frame of reference which is defined by 
the linear transformations given above, the relationship 


ee 

Ag, ™ Sb 

v¢) _ -—1,, — 2 
becomes HAH" y, = YW; 
or By, = Yo- 


Thus, in the new frame of reference the matrix B = HAH~ operating on the vector y, is 
equivalent to the matrix A operating on the vector €, in the original frame. 

It is convenient, for the purposes of studying the matrix A and of performing any numerical 
computations with it, to transform the variables y and € in the above way, choosing the 
matrix H so as to make B = HAH~‘ as simple as possible. For B! = (HAH-) = HA'H-, 
and since A is non-singular, by the reversal law, (HAH~)-! = HA-1H-1. Thus, if f(A) is 
a rational integral function of A, f(B) = f(HAH~) = Hf(A)H-; and the properties of 
matrix functions f(A) can be studied by means of the simpler forms f(B). Moreover, the 
matrices A and B have the same characteristic equation and, therefore, the same latent 
roots. For B—AI = H(A—AI) H~ and, forming the determinants of both sides, 

|B-AI|=|H||A-AI||H|-, 


so that the characteristic equation is 


| A-Al| =| B-Al| =0. 


in 
de 
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If, in the present case, the transforming matrix is taken to be 
(FFF... F.) : ie nha cea ; J 
; (FF F....2) : 
: ‘ (PP, «.. Fe») 
, é a, So Se4 
‘ i: ‘ |. ee 


in which, it is to be noted, the only numerical elements lie in the principal diagonal and are 
derived entirely from the life table, then 





L 








Fy PF, P,P,Fy P,P,P,Fy ... (PyP,Py--- Pea) Fe 
1 : 
1 : 
B=HAH" =| | 1 ; 
1 
Bis 1 | 





Comparing this matrix B with the original form A, it can be seen that the latter has been 
simplified to the extent that the original P, figures in the principal subdiagonal are now re- 
placed by a series of units, and the matrix A has been reduced to the rational canonical form 
B = HAH“ (see Turnbull & Aitken, 1932, chap. v). In this way any computations with the 
matrix A are made easier, and we may work henceforward in terms of ¢ and y vectors 
together with the matrix B, instead of with the original y and £ vectors, and the matrix A. 
Any results obtained in this new system of co-ordinates may be transformed back again 
to the original system whenever necessary. It is evident that by suitably enlarging H the 
original matrix M may be transformed in a similar way. 

This linear transformation of the original co-ordinate system is equivalent biologically 
to the transformation of the original population we were considering into a new and com- 
pletely imaginary type which, although intimately connected with the old, has certain quite 
different properties. Thus, it can be seen from the transformed matrix B that the individuals 
in this new population, instead of dying off according to age as the original ones did, live 
until the whole span of life is completed, when they all die simultaneously. This is indicated 
by the P, figures being now all equal to unity; an individual alive in the age group x to x+ 1 
at ¢ = 0 is certain of being alive at ¢ = 1, excepting in the last age group of all where none of 
the individuals will be alive in one unit’s time. Accompanying this somewhat radical change 
in the life table, there is a compensatory adjustment made in the rates of fertility so that the 
new population has the same inherent power of natural increase (7) as that of the old. This 
follows from the fact that the latent roots of the matrices A and B are the same, and, as 
will be shown later, the dominant latent root is closely related to the value of r obtained by 
the usual methods of computation. Insomuch as the transformation is reversible and 
A = HB3H, it can be seen that by changing H we could transform the canonical form B, 
if we wished, into another matrix in which the P, subdiagonal might be a specified set of 
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figures derived from some other form of life table. But, for our present purposes, the canonical 
form B, in which all the P, figures are units, offers advantages over any other matrix of a 
similar type owing to the greater ease with which it can be handled. 


6. RELATION BETWEEN THE CANONICAL FORM B aND THE L,m, COLUMN 


The actual computation of the matrix B by way of the steps indicated in the theoretical 
development is by no means difficult, although it is a somewhat tedious process, particularly 
if the matrix is of a large order. The numerical elements in the first row of B for the brown 
rat are given in Table 5 of the Appendix. These values were obtained from the F, and P, 
figures which have already been used in the numerical example in § 3 and which will be found 
in the same table. Further reflection suggested, however, that instead of first of all obtaining 
A and then transforming to B, a short cut could be taken which would save labour and 
which also would tend to eliminate some of the small cumulative errors arising in the longer 
method. 

The series of values P,, P,P,, Py P, Py, ..., (PoP, Pe--. Pes) by which the individual F, 
figures are multiplied in order to obtain the first row of B, is essentially a stationary age 
distribution. For, since by definition, 





| rey 
P, = _ 
(P, P,P, ... P.) = es, 
0 


1 
where L, = | l.dx. Hence the required -series of multipliers is given by a stationary age 
0 


distribution in which only one individual is alive in the age group 0-1. Now, the F, figures, 
as defined in § 2, already contain within them some allowance not only for the probability 
of survival during the first unit of life, but also for the fact that some adult individuals in 
each age group are dying off during the interval of time 0-1. The process of multiplying F, 
by (P,P, P, ... P,_,) is thus analogous to the formation of the Z,m, column, by means of 
which the net reproduction rate is estimated. The chief difference between the first row of B 
and the L,m, distribution is that in the former the maternal frequency is expressed as 
between the ages of x+ } to x+ 14, instead of between x to x+1 as in the latter. If each 
element (P,P, P,... P._,F,) of the first row of B is regarded as centred at the age of x+ 1, 
the sum, mean and seminvariants of this ‘distribution’ may be estimated and compared 
with the values which are obtained from the Z,m, column in the process of calculating r by 
the usual methods. In the present numerical example the results of this comparison were 
as follows: 














Parameter L,m, column | First row of B 
| | 
Sum (R,) | 25-65786 | 25-6603 
ean 9-60604 9-5948 
Ms | 14-14397 | 14-1839 
m, 22-15696 21-9358 
m,—3mi | — 117-6480 | — 117-920 





After allowing for the small cumulative errors which might be expected to occur in the 
calculation of the matrix elements, there is a substantial agreement between the respective 


a 
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estimates. This agreement strongly suggests that if we had wished to pass immediately to 
the matrix B without going through the laborious process of calculating the F, and P, 
figures, the elements of the first row could have been obtained by forming a new L,m, 
column in which the age group limits were shifted a half unit later in life. This could readily 
be done by interpolating in the integral curve of the Lm, values for the ages x + 3. This 
method of forming the first row of B has been adopted in other instances, when the matrix A 
was not of any immediate interest. It proved to be relatively quick and certainly less 
laborious than the method of first establishing A and then transforming to B which was the 
one used in the present numerical example. 


7. THE STABLE AGE DISTRIBUTION 


The result of operating on an age distributon y, with the matrix B is, in general, a different 
distribution y,. But, in the special case when the relation between the two distributions is 


such that By, = Avan 


where A is an algebraic number, then y, may be said to be a stable age distribution appro- 
priate to the matrix B. For the sake of brevity it will be referred to as a stable y. Similarly 
for mnitinl sow vectors, if $,B =Adbu 
then ¢, is said to be a stable ¢. 
The matrix equation defining a stable y may be written as k + 1 linear equations, of which 

the ith is - k+l 

2; bin; — An; = 0, 

j=1 


where n; (i = 1, 2,...,4+41) are the co-ordinates of the stable y, and b,; the element in the 
ith row and jth column of B. Eliminating the n; from this system of equations, we obtain 
the characteristic equation of B, namely, 


| B—AI| = 0; 
and, expanding this determinant in powers of A, we have in the present case, 
Ak+1— FoA* — PFA‘ — PP, BA*? - ... — (Py Py... Pre) Fy A— (PoP, --- Py_i) F, = 9. 


The k+ 1 roots A, of this equation are the latent roots of B, and corresponding to each dis- 


tinct A, there is a pair of stable vectors, ¢, and y,,, determined except for an arbitrary 
scalar factor. 


Once a latent root A, has been determined, it is a comparatively simple matter to find the 
appropriate stable y,, and ¢, vectors. Thus, it is easily shown that the stable y, is the column 


vector {AKAE-1Ak-2 A, 1}. A short method of estimating ¢, is the following. Suppose, to 
take a simple case, that 
abed 


: 
and let y, (x = 1, 2,3, 4) be the elements of the stable ¢, appropriate to the root A,. Then 
$B =[ay:t+¥. by:+¥3 CYity, dy] 
=[ Acti Auta Aas AuYal: 
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° ge 6 ; C+i ts 
By equating similar elements and putting y, = 1, y, = d/A,, y3 = —*, etc., it is easy to 


. 
see how the required row vector can be built up. Having in this way obtained the stable 


y and ¢ vectors for the matrix B, they may be transformed to the appropriate stable € 
and » for the matrix A by means of the relations 
n=$H, §= Hy. 

The characteristic equation of the matrix B, when expanded, is of degree k +1 in A, and 
once B has been obtained this equation can immediately be written down, since the numerical 
coefficients of A*, A*—1, A*-2, etc., are merely the elements of the first row taken with a negative 
sign. Since there is only one change of sign in this equation, only one of the latent roots 
will be real and positive. Excluding the rather special case when the first row of B has only 
a single non-zero element, and taking the more usual type of matrix which will be met with, 
namely, that for a species breeding continuously over a large proportion of its total life span, 
it will be found that the modulus of this root (A,) is greater than any of the others, 

|Ay|>|Ag|>]Ag|>--->]Agual, 
the remaining roots being either negative or complex. 

This dominant latent root A,, which will be 21 according as to whether the sum of the 
elements in the first row of B is 21, is the one which is principally of interest. Since it is 
real and positive, it is the only root which will give rise to a stable y or £ vector consisting 
of real and positive elements. It is this stable £, associated with the dominant root A, which 
is ordinarily referred to as the stable age distribution appropriate to the given age specific 
rates of fertility and mortality. Since 

Af, = AiEs, 
it can be seen that the latent root A, of the matrix A and the value of r obtained in the usual 
way from w 
| el m,dx = 1, 

0 
are related by log, A, = r. 
From the mathematical point of view, however, the negative and complex roots of the 
characteristic equation are of importance in the further theoretical development. Moreover, 
as will be shown later, the stable vectors associated with them are not entirely without 
interest. Two main cases then arise: when the remaining roots are all distinct, and when 
there are repeated roots. For the present it will be assumed that the latent roots of the 
matrix are ali distinct. 


8. PROPERTIES OF THE STABLE VECTORS 
Before proceeding further it is necessary to mention briefly the reasons why the methods 
given above for the computation of the stable y and ¢ vectors were adopted, apart from their 
simplicity in practice. If the k+1 distinct roots of the characteristic equation are known, 
we may form a set of k + 1 matrices f(A.) by inserting in turn the numerical value of each root 
in the matrix [B—A,, J]. The adjoint of f(A,) is 
F(A,) = Tis-” 1} and f(A,) F(A,) = 0. 

It may be shown that the stable y,, appropriate to the root A, can be taken proportional to 
any column, and the stable ¢, proportional to any row of the matrix F(A,) (see e.g. Frazer, 
Duncan & Collar, 1938, chap. m1). Moreover, F(A,) is a matrix product of the type W¢, 


10) 


on 


ial 
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where the y vector is given by the first column and the ¢ vector by the last row of F(A,), 
each divided by the square root of the element in the bottom left-hand corner; and the 
trace of the matrix is equal to the scaiar product dy. Now [B—A, J] is a square matrix of 
order k + 1 with only zero elements below and to the left of the principal subdiagonal, which 
itself consists of units. The product of k such matrices, which gives F(A,), will have therefore 
a unit in the bottom left-hand corner. Since the stable ¢, and y, vectors obtained by the 
methods suggested in §7 have respectively their first and last elements = 1, it follows that 


Vaba ae F(A,), Pa Va = trace F(A,). 
The stable vectors may now be normalized. If the scalar product, ¢, 7, = 2*, say, then 


J, 
Pa Va et 1 


lz[Jz| 


From now on it will be assumed that the stable vectors appropriate to each of the latent 
roots have been normalized in this way. 


These vectors have the following important properties: 


(1) The k+1 stable y are linearly independent. There is thus no such relationship, with 
non-zero coefficients c, as 


CY tCgWat sgt... + Cpe = 9. 
(2) The scalar product of a stable w, ww, with the associated vector of another stable 
] Y, Va 
We, Wr is zero, i.e. $,%e=0 (a+b). 


The normalized stable y thus form a set of k+1 independent and mutually orthogonal 
vectors of unit length. 





3) Any arbitrary yy—v,,, say—can be expanded in terms of the stable y, thus 
: Y ¥—¥Vz 8a} I y 


Wn = CW t Cag Cg Wat... + Cp Wer 


where the coefficients c may be either real or complex. Similarly an arbitrary vector ¢, can 
be expanded in terms of the stable ¢. 


9. THE SPECTRAL SET OF OPERATORS 


The matrix product y,¢, of the normalized stable vectors associated with the latent root 
A, will be termed the matrix S,. From the relationships which have already been given, it 
can be seen that S, is merely the adjoint matrix F(A,) of the previous section after each 
element in the latter has been divided by the sum of the elements in the principal diagonal ; 
in other words it is the normalized F(A,). In the case of all the latent roots being distinct, 
there are thus k + 1 matrices S,, and these S, form a spectral set of operators with the following 
properties: k+1 
S=§, §S&=-0 (a+b), YS =f. 

a=1 
Moreover, if f(.B) is a polynomial of the matrix B, we have by Sylvester’s theorem (Turnbull 
& Aitken, 1932, chap. v1, §8) ; 


so that the matrix 


B= A,8,+AgSe4+..- + Ags Segr, and Bl = ALS, +ALS.4+ ... +AggaSpar- 
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If the latent roots in the expansion of B are raised to a high power, the term associated with 
the positive real root predominates over all the others, so that when ¢ is large, we have 
approximately B= ALS. 


In any particular case the power to which B will have to be raised in order that this equation 
should be approximately true, will depend both on the order of the matrix and on the relative 
magnitude of the dominant root as compared with that of the remaining roots. 

At this point it is possible to attach some biological meaning to one of the ¢ or 7 row vectors, 
which in the first place were introduced into the theory for reasons of symmetry, and which 
were defined solely in terms of their mathematical properties. If at a given moment a trans- 
formed population has an arbitrary age distribution y,, and the sequence By, B*y,, ..., By, 
is formed, it can be seen that when ¢ is large and y, is expanded in terms of the stable y, we 
have approximately BY, = c,ALy,. 


Thus, a population with any arbitrary age distribution tends ultimately to approach the 
stable form appropriate to the given rates of fertility and mortality, provided that these 
age-specific rates remain constant. This theorem is, of course, well known; and it is clear 
that the achievement of the stable form of age distribution associated with the dominant 
latent root is very unlikely to occur in practice, except in the case when the initial distribu- 
tion is already of that form or exhibits only small departures from it. Now, it has already 
been shown that the sums of the columns of a matrix BY provide a measure of the contributions 
made to the population at time ¢ per individual alive in the respective age groups at t = 0. 
When ¢ is large, the matrix B' is equivalent to the matrix S, multiplied by a scalar factor. 
From the way in which this latter matrix was constructed by the outer multiplication of 
y, and ¢,, it is evident that the sums of the columns of S, are proportional to the vector ¢,. 
Thus, transforming back again to the original co-ordinate system, the stable 7, associated 
with the dominant latent root provides a measure of the relative contributions per head 
made to the stable population by the individual age groups. 


10. Repvuction oF B TO CLASSICAL CANONICAL FORM 


From the k+1 stable y a matrix Q can be constructed, whose columns are the stable ¥ 
arranged, reading from left to right, in descending order of the moduli of the roots with which 
they are associated. Corresponding to every pair of complex roots, u+iv, there will be in 
this matrix a pair of columns consisting of complex elements, the one column being the 
conjugate complex of the other. Some of the columns associated with the negative roots 
may be purely imaginary owing to the normalization of the corresponding y and ¢ vectors. 
In a similar way a matrix U may be formed, whose rows, reading from above down, are the 
stable ¢ arranged in the same order. Since the stable ¢ and y are normalized, and ¢, y, = 0 
for a+b, Ug =I, 


and, therefore, U and Q are reciprocal matrices. By premultiplying and postmultiplying © 


respectively with U and Q, the matrix B may be reduced to the classical canonical form C, 
in which the only elements lie in the principal diagonal and consist of the latent roots 
arranged in the order prescribed above. This reduction of B to a purely diagonal form by 
means of the collineatory transformation UBQ = C is, however, only possible in the type 
of matrix we are considering, when the latent roots are all distinct. 


nanan anlUreet,lC Or lCUC KM CféCoe 
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The expansion of an arbitrary y, in terms of the stable y, 


Vn = Wy tegpet--- +p Wer 
may be written in matrix notationas y= Qc, 


where c is the column vector {c,¢,C,...¢,,,}. Similarly, the expansion of the vector ¢, 
associated with y, may be written ¢ = dU 


where d is the row vector [d,d,d,...d,,,]. This is again a transformation to another co- 
ordinate system, but this time the reference axes are at right angles to one another. Since 
the variables transform contragrediently, dc = dy = 7. At this point it is necessary to 
make some assumption as to the relationship between the elements of the vectors d and c. 
Since these elements may be either real or complex, it will be assumed that 


d=@', 


where the row vector Cc’ is the transposed conjugate complex of the column vector c. Hence, 
the square of the length of a vector referred to this orthogonal co-ordinate system is given 
by ¢’c, a number which is essentially real and non-negative. (The assumption that d = c’ 
will be found, in the particular case studied here, to lead to values of c’c which, although 
real, may be negative.) 


11. THE RELATION BETWEEN ¢ AND VECTORS 


Since 7, = Qc,, and the associated ¢, = ¢, U, it may be seen from the relations given in the 
two previous sections that z= 0’ Uy, = Gy 7 
—_— ed =z 


The matrix G = U'U is symmetrical and all its elements are real numbers, those in the 
principal diagonal being necessarily positive in sign. It therefore remains unaltered after 
transposition. Since the elements of the vector y,, which is by definition an age distribution 
transformed by the matrix H, are also necessarily real, we may write 


¢: a 2G. 
The role of the matrix G is therefore the same as that of the double covariant metric tensor 
Jmn in the tensor calculus. It transforms any y vector into its associated ¢ vector. This 
process is reversible, the reciprocal matrix being given by G-! =.QQ’. 
The magnitude of a vector y,, is defined by the equation 


x = (¥Gy,)t, 
where the square root is taken with a positive sign. If we have two vectors y, and y,, both 
radiating from the common origin, the angle between them is given by 

cos @ = vey > 

xy 

from which it follows that when ¥,Gy, = ¢,~, = 0, the two vectors are at right angles to 
one another, and when cos@ = | their directions are the same. If, in the last equation, we 
take y, to be the stable vector y, associated with the dominant latent root A,, then knowing 
the magnitude of a vector BY, and the angle which it makes with the y, axis, we can obtain 


a graphical representation of the way in which a particular age distribution approaches the 
stable form. 
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The matrix G also defines the angles between the reference axes of the co-ordinate system. 


If we introduce into the vector space of the w’s a system of reference axes defined by the 


unit column vectors e, ~ ss 0} 


e ={f0 10 .. 0 
e, ={f0 0 1 ... 0} 
Ennw 6. 6°. 9 
then the distance from the origin of the unit point e; is 
Oe; = vyi5» 
and the angle between any two of the co-ordinate axes is given by 
cos 0;; = ij 
ae (Giz 933) 


where, in both cases, g;; is the element in the ith row and jth column of G. 


By transforming back to the original co-ordinate system, the metric matrix associated 
with the vector space of the é’s will be found to be 


G, = HG,H. 


Hitherto we have been cliiefly concerned with an operator B which, acting in the vector 
space of the ’s, has the power of transforming a vector y, into what is in general a new vector 


y,. We may now have reason to inquire how the associated vectors ¢, and ¢, are related in 
the vector space of the ¢’s,. whenever 


By, -_ Vy 
Since y¥,=G"¢, and y, = G4, 
we have GBG"¢, = $1, 


and hence, by transposing, 
¢,G"B’G = ¢,,. 


Thus, the matrix which transforms ¢, into ¢, is not the same as that which transforms yy, 
into y,. In order to distinguish these two operators, they will be referred to as B,; and B, 
respectively. Ir the few numerical examples which have been worked out, the matrix B, 
differed greatly from the rational canonical form B,, and consisted of (k + 1)? real elements, 
some of which were negative. In addition to the relationship 


Bs = GB, G, 
it was also found that in the case of distinct latent roots 
By.= A, S,+A,S,+AgS3+ eee + Apis Spit 


where the S, matrices are the spectral set of operators defined in §9 and A, is the conjugate 


complex of the latent root A,,. It may be seen from this expansion of B, that the necessary’ 


condition for B, = B, is that all the latent roots of B should be real, the one positive and the 
remainder negative. (It is to be noted that we are dealing here with the case of distinct 
latent roots; it would appear that even if all the A were real, B, + B, in the case of repeated 
roots.) Unless, however, the matrix B is of a small order, it is unlikely that this condition 
would be fulfilled, since in the more usual-sized matrix we shall be dealing with in the case 
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of human or other mammalian populations, some of the roots will almost certainly be 
complex.* 

It seems unlikely that the equations given in this section will have very much practical 
application at the moment; they have been included merely to fill in the picture of the 
relationship between the two types of vector. For all ordinary purposes no one would choose 
to work in terms of ¢ vectors and the operator B,, instead of the more obvious y vectors 
and the more simple matrix form B,. Nevertheless, since it has been necessary to assume 
that there are such vectors as 7 or ¢ associated with every £ or y, and since these vectors 
play such an important part in the mathematical theory, the question naturally arises as to 
what significance must be attached to them from the biological point of view. Have they in 
fact any real meaning at all? Or must they be regarded purely as mathematical abstractions? 
At the end of §9 it has been suggested that the row vector associated with the stable £, 
appropriate to the dominant latent root is a measure of the contributions made to the stable 
population per individual female alive in the respective age groups of the initial distribution; 
but this is a special case and the interpretation offered here is not applicable, even in a wider 
form, to 7 vectors in general. It may well be, of course, that the latter as a class have no 
concrete meaning: and that in seeking to define them in terms of some property or character- 
istic of an age distribution one is merely attempting the impossible. But the fact of one 9 
vector having been defined in non-mathematical terms, even though on further consideration 
some revision may be needed of the actual definition given here, suggests that impossible 
may perhaps be too final a word to use in this connexion. 


12. CASE OF REPEATED LATENT ROOTS 


When any of the latent roots other than the real positive dominant root are repeated, a 
number of the relations given in the previous sections no longer hold good and certain 
equations must therefore be modified. Suppose a root A, has a multiplicity s, and consider 
the matrix f(A,,) such as, to take a simple example, 


a—A, 6b c d 


‘ i .—--* 0 
f( a) = 0 l ra da 0 
0 pe a 


Then, since the determinant | f(A,) | = 0 and at least one of the first minors of order 3 is not 
equal to zero, the above matrix has rank 3 and, therefore, nullity 1. Hence it can be seen that 
f(A,) of whatever order it may be has nullity 1. Since f(A,) is thus simply degenerate, there is 
only one stable y appropriate to the s equal roots A, (see e.g. Frazer, Duncan & Collar, 1938, 
chap. II). 

Certain consequences immediately follow. Since the matrices Q and U cannot be con- 
structed in the way given in § 10, the reduction of B to a purely diagonal matrix by means 
of the collineatory transformation U BQ can no longer be carried out. Neither is the expan- 
sion of B in terms of the spectral set of S,, matricés, nor the expansion of an arbitrary y, 

* The interesting theoretical case of the matrix A or B having a number of its latent roots real and positive, 
with the remainder real and negative, is outside the scope of the present study. The necessary conditions for 
this to be true would involve a number of the F, figures becoming negative, a case not considered here, but which 


biologically might be held to correspond with the destruction of eggs, or the very young, by certain age groups, 
e.g. as observed by Chapman (1933) in experimental populations of the flour beetle, Tribolium confusum. 
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in terms of the stable y, possible in the actual forms given in §§9 and 8. We may, however, 
proceed as follows. 


The matrix Q is essentially an alternant which has been postmultiplied by a diagonal 
matrix N, the elements in the latter being given by the reciprocals of the scalar factors | z | 


by which the stable vectors were divided in the process of normalization (= A = 1) , 
so that Q = XN, where 

r A‘ ao 
a we? 4 a 
X=| : : : 
Ay Ry ae? 
aoe ae ee 


When a root A,, is repeated s times, s of the columns in X become the same, and therefore 
the matrix becomes singular. In place_of this alternant matrix we have the confluent 
alternant form (see Turnbull & Aitken, 1932, chap. v1) in which the s columns (s = 0, 1, 2, 3, 
...,8— 1) corresponding to the repeated root A, are written 








yk EAk-2 ke(e — 1) AR? 
a "a 2! 
t-1 (LE -2 (k—1)(k-2) yy. 
AK-1 (k—1) Ak ee : 
. 3A, 
2A, 1 
ra 1 0 
l 0 0 


the column s being obtained from column 0 (the non-normalized stable y,,) by the operation 
8 
( = ) / s!. This confluent alternant form of X is non-singular and therefore a reciprocal 


a 


matrix can be determined (see Aitken, 1939, § 50). The general classical canonical form of B 
obtained by the collineatory transformation, X-!'BX = C, has corresponding to the repeated 
root A, a diagonal submatrix: 


See eee eee eeereeeeeeeeseeseeese 








The matrix product of a column of X with the appropriate row of the reciprocal X-! 


gives as before S, and the k+1 8, form a spectral set with the same properties as those 
defined in § 9, except that 
k+1 


B+ > A,S,,. 
a=1 
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In place of this expansion of B in terms of the S, matrices, we have in the case of repeated 
roots the confluent form of Sylvester’s theorem, for details of which reference may be made 
to Frazer, Duncan & Collar, 1938, chap. m1. Apart from this modification, however, we may 
obtain by inspection of the S, the factors by which the respective columns of X must be 
divided in order to express this matrix in a form comparable to that of Q in § 10. Similarly, 
when the respective rows of X~! are multiplied by the appropriate factors, the matrix U is 
found and hence G = U’U can be constructed. (It is to be noted that (X—1)’ X- is neither 
equal to, nor directly proportional to U’'U.) An arbitrary y, can be expanded in terms of 


the column vectors of Q, though in the case of only one column associated with the repeated 
root A, does the relationship By, = A,y, hold. 





13. THE APPROACH TO THE STABLE AGE DISTRIBUTION 


A stable age distribution appropriate to the matrix B has been defined mathematically by 
the equation By = Ay, 


and it has already been shown that since only one latent root of B is real and positive, only 
one of the stable y will consisi of real and positive elements. But, in addition to this Mal- 
thusian age distribution, it is also of some interest to inquire whether any significance can 
be attached to the remaining stable y associated with the negative and complex roots of the 
characteristic equation. 

Any age distribution y,, the elements of which are necessarily 20, may be expressed as 
a vector of deviates from the stable y/, associated with the dominant latent root, and we may 
therefore write the expansion of y, in terms of the stable y as 


(Yr—C1Wi) = CoWa t+ Cg Wat... + Cpa Vers = Var 


where the coefficients c are given by the vector c = Uy,. Thus, the way in which a particular 
type of age distribution will approach the stable form may be studied by means of the 
vector Wy. 

Among the terms occurring in the right-hand side of this expression there will be, corre- 
sponding to each negative root, a single term c, y, which will consist of real elements alter- 
nately positive and negative in sign. (Even if the normalized y, is imaginary this term will 
consist of real numbers, since in this case c, will also become imaginary.) Moreover, corre- 
sponding to every pair of complex roots there will be a pair of terms (c,,¥,,+¢,W,) which 
taken together will also give a single vector with real elements. This follows from the fact 
that c,,, is the conjugate complex of c,, owing to the way in which the matrix U is constructed. 
Then, apart from the scalar c, which must necessarily be > 0, some of the remaining coeffi- 
cients Cg, Cy, ..-, C,,, in the expansion of yz may be zero. The first and most obvious case is 


when they are all zero, and the age distribution y, is therefore already of the stable form. 
But, if either Va = Ca¥ 
ad ava 


where y, corresponds to a negative latent root, or 


Va _ Cn¥m +Cy Vn 


where y,, and y,, are associated with a conjugate pair of complex roots, then it follows that 


the age distribution y, will, as time goes on, approach the stable form in a particular way 
defined by either 


Bg = CaX'Wa or BYyg= CmA Yin + CrAW ys 


Biometrika 33 
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in which A‘ for a pair of complex roots w+iv with modulus r may be written in the form of 
r'(cos Ot + isin Ot). Thus, the negative and complex latent roots of B serve to determine a 
number of age distributions which are of some interest owing to the fact that they will 
approach the Malthusian form in what may be termed a stable fashion. 

Since |A,|>|A,|>|A3|>-..>|Agyi|, the vector of deviates yy will tend towards zero 
as t->0o whenever A, = 1. Thus, in the case of a stationary population, any y, will converge 
to the stable form of age distribution. But if A, > 1, there is a possibility of one or more of 
the remaining roots having a modulus 21, e.g. | A,|21. In the latter case there may be 
certain age distributions with c, + 0 for which the amplitude of the deviations from the stable 
form tend either to increase (| A, | > 1), or to remain constant (| A, | = 1). From the practical 
point of view, however, we may still say that a population with such an age distribution 
approaches or becomesapproximately equal to the stable population, since A{ is much greater 
than A‘ when ¢ is large. 


14. SPECIAL CASE OF THE MATRIX WITH ONLY A SINGLE NON-ZERO F, ELEMENT 


The interesting case of the matrix A having only a single non-zero element in the first row has 
been illustrated in a numerical example by Bernardelli (1941).* This author has also used 
a matrix notation in the mathematical appendix to his paper, and the form of his basic 
matrix is the same as that referred to here as M or A. It is not clear, however, from the 
definitions which he gives whether he regards the elements in the first row of his matrix as 
being constituted by the maternal frequency figures (m.,) themselves, or by a series of values 
similar to those defined here as the F, figures. He refers to them merely as the specific 
fertility rates for female births. ; 

In discussing the causes of population waves, Bernardelli describes a hypothetical species, 
such as a beetle, which lives for only three years and which propagates in the third year of 
life. He assumes, for the sake of argument, that—to employ the terminology used here— 
P, = } and P, = , and that ‘each female in the age 2-3 produces, on the average, 6 new 
living females’. Assuming for the moment that he is here defining a F, figure, we may, write 
this system of mortality and fertility rates as 


00 6 001 
A=|} 0 0], B=HAH-=]1 0 oO}. 
0240 0 1 


The characteristic equation expanded in terms of A is A?— 1 = 0; and the latent roots are 


3B 


3? , all three being of equal modulus. The matrix A has the interesting 
properties A4§*=A-", A'= I, 


therefore 1, —}+ 


so that any initial age distribution repeats itself regularly every three years. Thus, as 
Bernardelli shows, a population of 3000 females distributed equally among the three age 


* At this point I should like to acknowledge the gift of a reprint of this paper, which was received by the 
Bureau of Animal Population at a time when I was in the middle of this work, and when I was just beginning 
to appreciate the interesting results which could be obtained from the use of matrices and vectors: also a personal 
communication from Dr Bernardelli, received early in 1942, at a time when it was difficult to reply owing to the 
developments of the war situation in Burma. Although the problems we were immediately interested in differed 


somewhat, this paper did much to stimulate the theoretical development given here, and it is with great pleasure 
that I acknowledge the debt which I owe to him. 
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groups becomes a total population of 6833 at ¢ = 1; of 5166 at ¢ = 2; and again 3000 dis- 
tributed equally among the age groups at t = 3. Unless a population has already an initial 
age distribution in the ratio of {6: 3:1}, no approach will be made to the stable form associ- 
ated with the real latent root, and the vector of deviates £, will continue to oscillate with a 
stable amplitude, which will in part depend on the form of the initial distribution. Although 
this numerical example refers specifically to a stationary population, it is evident that a 
similar type of argument may be developed in the case when | A|> 1 and A? = AJ. 

We have assumed here that his definition of the fertility rate refers to a F, figure. But, if 
we were to interpret the words quoted above as referring to a maternal frequency figure, 
namely that every female alive between the ages 2-3 produces on the average 6 daughters 
per annum, then the results become entirely different. For, deriving the appropriate F, 
figures by the method described in § 2, the matrix is now 


oe Ss 0 4 3 
A=]} 0 0], B=HAH“=|1 0 O}, 
0 i 0 0 1 0 


and the latent roots are 1, — $+ 4%. The modulus of the pair of complex roots is 1/v2, which 
is < 1, so that every age distribution will now converge to the stable form associated with the 
real root. Thus, to take the same example as before, 3000 females distributed equally among 
the age groups will tend towards a total population of 4000 distributed in the ratio of 
{6:3:1}, and it was found that this age distribution would be achieved at approximately 
t = 23. During the approach to this stable form periodic waves are apparent both in the age 
distribution and in the total number of individuals, but these oscillations are now damped, 
in contrast with the results obtained with the first type of matrix. 

This simple illustration serves to emphasize the importance which must be attached to 
the way in which the basic data are defined and to the marked difference which exists between 
what are termed here the m, and F, figures. Nevertheless, apart from the question of the 
precise way in which the definition of the fertility rates is to be interpreted in this example, 
the first type of matrix with only a single element in the first-row does correspond to the 
reproductive biology of certain species. Thus, in the case of many insect types the individuals 
pass the major portion of their life span in various immature phases and end their lives in a 
short and highly concentrated spell of breeding. The properties of this matrix suggest that 
any stability of age structure will be exceptional in a population of this type, and that even 


if the matrix remains constant we should expect quite violent oscillations to occur in the 
total number of individuals. 


15. NUMERICAL COMPARISON WITH THE USUAL METHODS OF COMPUTATION 


From the practical point of view it will not always be necessary to estimate the actual values 
of all the stable vectors and of the associated matrices which are based on them. Naturally, 
much will depend on the type of information which is required in any particular case. In 
order to compute, for instance, the matrices U, Q and G, it is necessary first of all to deter- 
mine all the latent roots of the basic matrix. The ease with which these may be found depends 
very greatly upon the order of the matrix. Thus, in the numerical example for the brown 
rat used previously in §3, the unit of age and time is one month and the resulting square 
matrix A is of order 21. To determine all the 21 roots of the characteristic equation would be 
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a formidable undertaking. It might be sufficient in this case to estimate the positive real 
root and the stable vector associated with it. On the other hand, it is possible to reduce the 
size of the matrix by taking a larger unit of age, and in some types of problem, where extreme 
accuracy is not essential, a unit say three times as great might be equally satisfactory, which 
would reduce the matrix for the rat population to the order of 7 x 7. It is not too difficult to 
find all the roots of a seventh degree equation by means of the root-squaring method (Whit- 
taker & Robinson, 1932, p. 106). But the reduction of the matrix in this way will generally 
lead to a value of the positive real root which is not the same as that obtained from the larger 
matrix, and it is therefore necessary to see by how much these values may differ owing to 
the adoption of a larger unit of time. 

Another important point which must be considered is the following. By expressing the 
age specific fertility and mortality rates in the form of a matrix and regarding an age dis- 
tribution as a vector, an element of discontinuity is introduced into what is ordinarily taken 
to be a continuous system. Instead of the differential and integral calculus, matrix algebra 




















Table 2 
| | 
Age group | ‘True’ Matrix Age group | ‘True’ Matrix 
(units of | stable age stable age (units of stable age stable age 
30 days) | distribution | distribution 30 days) distribution | distribution 
—_—__—_—_| re Sec 7 Prnas a 
0- 37,440 37,362 12- 161 160 
1l- 22,595 22,644 13- 101 101 
2- 14,417 14,444 14- 64 63 
3 9,227 | 9,238 15- 40 40 
4- 5,903 5,906 16- 25 24 
5- 3,775 3,775 17- 15 15 
6- 2,413 2,412 18- 9 9 
7- 1,542 1,540 19- 6 5 
8- 984 982 20- 3 3 | 
9- 627 626 | 
10- 399 398 “aa 
1l- 254 | 253 Total 100,000 100,000 | 

















is used, a step which leads to a great economy in the use of symbols and consequently to 
equations which are more easily handled. Moreover, many quite complicated arithmetical 
problems can be solved with great ease by manipulating the matrix which represents the 
given system of age specific rates. But the question then arises whether these advantages 
may not be offset by a greater degree of inaccuracy in the results as compared with those 
obtained from the previous methods of computation. It is not easy, however, to settle this 
point satisfactorily. In the way the usual equations of population mathematics are solved, 
a similar element of discontinuity is introduced by the use of age grouping. Thus, in the case 
of a human population, if we were estimating the inherent rate of increase in the ordinary 
way, we should not expect to obtain the same value of r from the data grouped in five year 
intervals of age as that from the data grouped in one year intervals. The estimates of the 
seminvariants would not be precisely the same in both cases. Nevertheless, the estimate 
from the data grouped in five year intervals is usually considered to be sufficiently accurate 
for all ordinary purposes, and there is little doubt that if we merely require the inherent rate 
of increase and the stable age distribution, these methods of computation are perfectly 
satisfactory when applied to human data. But, in the case of rodents, and probably also 
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other species with high gross and net reproduction rates, it will be found that even a 4th 
degree equation in r with the coefficients based on the seminvariants of the L_m, distribution 
is, in many examples, not sufficient to give an accurate estimate of the rate of increase, and 
it is necessary to arrive at a better value of r in a somewhat roundabout way. Here, the 
determination of the positive real root of the characteristic equation for the matrix, once the 
latter has been established, may be even quicker than finding a solution from the L,m, 
column by a method such as that described in the appendix. 

In order to compare the values of r obtained from the characteristic equation of the matrix 
with those obtained from the L,m, column, both methods were used in the numerical 
example for the brown rat, and a comparison was also made between the values when the 
data were grouped in 1 month and in 3 month age intervals. In addition, the stable age 
distribution appropriate to the positive real root of the matrix was also calculated in both 
cases. The results were as follows. 

(a) One month unit of grouping ; matrix of order 21 x 21. Using the method of computation 
indicated in the appendix, the value of r was estimated to be 0-44565 per month of 30 days. 
The positive real root of the characteristic equation was A, = 1-56246, whence r = 0-44626, 
a value which differs from the former only in the fourth decimal place. The appropriate age 


distributions, expressed per 100,000, are given in Table 2, the ‘true’ stable being obtained 
from 


z+1 
n, = 100,008 { e**l_ dz. 
z 


The agreement between these distributions is very good, although the one derived from 
the matrix shows certain small rhythmical departures from the ‘true’ distribution par- 
ticularly in the earlier age groups. The maximum difference between them in this region, 
however, is not greater than 2-2 per thousand. Since the matrix stable distribution is pro- 
portional to the columns of the matrix H-4S,H, which in turn is proportional to A‘ when ¢ 
is large, this agreement between the two distributions also indicates that the cumulative 
errors which might be expected in forming the series A*, A*, A‘, etc., owing to the P, figures 
being based on the life table age distribution, are not very serious. Judging by this example 
it seems that satisfactory estimates of the inherent rate of increase and of the stable age 
distribution may be made from a large order matrix. 

(b) Three months unit of grouping ; matrix of order 7 x 7. Clearly there are several ways in 
which a large order matrix may be condensed into one of a smaller order. The method which 
was used in the present instance was to construct the first row of the condensed canonical 
form B by interpolating in the integral curve of the original Z,m, column (1 month units 
of grouping) for the ages 4-5, 7-5, 10-5, etc., and taking the first differences of the seven 
values thus obtained. Since interpolation was not very satisfactory in the earlier part of 
the integral curve—the differences converged rather slowly in this region—the elements 
were expressed to only three places of decimals. (Some preliminary transformation of the 
integral L,m, figures might have been better in this case.) The characteristic equation, 
expanded in terms of A, was found to be 


A? — 1-756A® — 6-899A5 — 7-203A4 — 5-344A3 — 3244.2 — 1-110A—0-102 = 0. 


It will be seen that the sum of the coefficients, Ry = 25-658, which is necessarily the same 
as the original net reproduction rate owing to the way in which the coefficients were derived. 
For interest, the seven roots of this equation were then determined by the root-squaring 
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method, using 4-place tables of logarithms and Barlow’s tables of squares (Whittaker & 
Robinson, 1932, p. 110). The approximate values of the roots, arranged in descending order 
of their moduli, are 
A,= 4016, 
A, = — 1-032, 
AgA, = — 00215 + 0-67621 (mod. = 0-6765), 
A;Ag = —0°5245 + 0-3486i (mod. = 0-6298), 
A, = — 0-135. 


Thus, apart from the positive real root, there are two negative and two pairs of complex 
roots. It is interesting to note in passing that in this example the modulus of the second 
latent root is > 1 (vide §13). From the value of the dominant root we find r = 0-4634 per 
head per month of 30 days, an estimate of the rate of increase which is 1-71 % per month 
higher than that from the large order matrix. 


Table 3. Stable age distributions 




















Age group / wr From condensed Summation of 
(in months) aeaataead oe “*tade matrix ‘true’ distribution 
H- 75,960 75,762 74,452 
3- 18,055 18,267 18,905 
6- 4,514 4,519 4,939 
9- 1,120 1,111 1,280 
12- 273 267 326 
15- 64 61 80 
18- 14 13 18 
Total 100,000 100,000 100,000 














The net reproduction rate given by the new L,,.m,, column which was obtained by working 
in units of three months, was 25-6162, a figure somewhat lower than the original one of 
25-6579. The rate of increase, estimated in a similar way to the former example for one month 
age units, was r = 0-46034, again a higher figure, though one of much the same order as that 
obtained from the condensed matrix. 

The appropriate age distributions are given in Table 3, together with the ‘true’ distribu- 
tion of Table 2 summed in three month age groups. 

Compared with the last column, both of the stable distributions for the data grouped in 
three month age intervals are tilted towards the younger age classes, so that the number of 
immature females (< 3 months of age) is overestimated, while the remaining age groups are 
underestimated. The distributions derived from the integral and from the matrix are again 
of much the same order, and the differences between them and the last column, although not 
very great, are quite marked. ' 

The four estimates of the inherent rate of increase which have obtained from these 
numerical data may be compared in the following table. 

In both cases the estimates from the I,m, column and from the matrix agree very well: 
for a given unit of grouping both methods would seem to give comparable results. The 
differences between the estimates made by the same method are much greater, and the effect 
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of increasing the unit of grouping, and in this way shortening the labour, is to increase the 
value of r quite appreciably. Whether, or not, we should regard these estimates for the three 
months age grouping as satisfactory would depend on the degree of accuracy required in any 
particular calculation. It must be remembered, however, that the basic numerical data are 
of rather an extreme type in this example. It is doubtful whether any naturally living rat 
population would have so good a life table and so high a degree of fertility as that assumed 
for this imaginary population. In fact it was for these very reasons that these data were 








L,m, column | Matrix Difference 
RE RNA Ae pile ale: ele Reh: + piel — 
| 
| 1 month age groups 0-44565 0-44626 0-00061 
| 3 ee | 0-46034 0-4634 0-0031 
| | 
| : eel e mises Ri ee ee <i ta = La i 
| Difference 0-01469 0-0171 
ered Lian dtl A RRBs aM Aan RAs tae Ste ho ew 


chosen as the basis of the numerical calculations in this work. For, if it could be shown that 
the two methods of computation gave comparable results in this case, it was felt that an 
even better agreement should be obtained in less extreme examples and more particularly 
with data relating to populations whose rate of increase is nearer to the stationary state. 
Although in this example the larger unit of grouping leads to rather unsatisfactory estimates 
of the rate of increase and the stable age distribution, it seems probable that, for the reasons 
just given, the differences would be less for instance in the case of human data. Hence, the 
question of the unit to be adopted is likely to become of less importance in the type of data 
more commonly met with, though it would be necessary to work out an example for such a 
population in order to check this point. 

There is, however, one way to avoid this difficulty of the working unit for populations with 
a high relative rate of increase. For example, returning to the numerical data used here, 
supposing that it was necessary in some particular problem to have a fairly high degree of 
accuracy in the results, but that the work involved in manipulating the large order matrix 
of 21 x 21 was too excessive. It might be sufficient in the case we are imagining to know the 
age distribution of the population in three month age groups at some particular time in the 
future, which we will take to be a multiple of three. Then, once the real latent root (A,) of 
the large matrix and its associated stable vector have been determined, it is possible to 
construct a small order matrix of 7 x 7 which has A? as its dominant root and therefore the 
same real stable age distribution as the larger matrix, only expressed in three month instead 
of in one month age units. It is convenient to carry out the calculation in terms of the 
canonical form B and of y vectors. Having determined the dominant latent root and the 
stable vector for the larger matrix, the elements of the first three rows of B® are then written 
down and summed in columns. This can be done very quickly in the present example, where 
reproduction does not start until the age of 3 months, for the third row of B® is the same as 
the first row of B; the second row is merely the first row of B shifted one age group to the left; 
and similarly again for the first row. The sums of the columns are then weighted with the 
number alive in the appropriate age group in the stable population (y, vector), and 
by summing the weighted column totals in groups of three and taking the weighted 
mean, we obtain the elements of the first row of a 7 x 7 matrix which has A? as its domi- 
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nant latent root. Thus, the characteristic equation of the original matrix condensed in this 
way was 


A? — 1-5056A® — 6-4694A5 — 7-2047A4 — 5-5371A3 — 3-4537A? — 1-3451A — 0-1447 = 0, 


and, out of idle curiosity, all the seven roots were extracted in order to compare them with 
those of the previous example of a condensed matrix. The estimation in this case was carried 
out to a higher degree of accuracy. The results were: 


A, = 381452, 
A, = — 1-02526, 
AsA, = —0-5905 +0-3782i (mod. = 0-70125), 


AsAg = 00280 +0-6879i (mod. = 0-68847), 
A, = —0-15876. 


The dominant root of the original matrix was 1-56246 and the cube root of A, is 1-56248. 
The remaining roots may be compared with those given for the previous example. The two 
negative roots are very similar and, in the second case, the two pairs of complex roots appear 
to have changed places, the real part of one pair becoming positive instead of negative. 
Although the cube root of A, is equal to the dominav:* root of the original matrix, it is un- 
fortunately not true that a similar relationship holds for the remaining six values of A. 
There is, for instance, no negative latent root > 1 for the larger form in this actual example. 

This point, however, raises an extremely interesting question. For a given series of data 
a finite matrix of a relatively small order may be constructed, as in the first example given 
here of a condensed matrix. Supposing that the order of this matrix is increased step by step 
and that in each case the latent roots are found. Then, in this approach to an infinite matrix, 
how do the latent roots behave and what relation does the array of roots in each case bear to 
those of the preceding steps? For the purposes of co.aparison it will be necessary to express 
the roots in terms of some suitable unit of time, e.g. per month or per year. So far as the 
real positive root is concerned, it seems likely that the series of individual roots will approach 
nearer and nearer to a limiting value. For the root A, is the ratio N,,,,/N,, or the number of 
times the stable population has increased at the end of the interval of time h. Then, expressing 
A, in the chosen unit of time, we have A, = (A,)', and taking logarithms; 


log. Ni.n—log. N, 


log, A; = h ’ 





so that, when the interval of time becomes very small, corresponding to a matrix of a very 


large order, and h->0, the right-hand side of this equation approaches the limit os =p, 


the true instantaneous relative rate of increase of the stable population. This argument is 
put forward with a certain amount of diffidence; it is only too easy for the biologist to over- 
look some flaw which will be immediately obvious to the trained mathematician. But, even 
if it were a valid argument for the behaviour of the dominant root, it can hardly be extended 
in this form to the case of the remaining roots; and thus the main question is left unanswered. 
From the point of view of the biologist, it would be interesting to know whether with an 


increase in the size of the matrix the array of secondary roots tends to coalesce round 
certain values of A*, 
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15. FURTHER PRACTICAL APPLICATIONS 


If we wish merely to estimate the inherent rate of increase and the stable age distribution 
appropriate to some system of age specific fertility and mortality rates, there is evidently 
little to choose between the matrix and the ordinary methods of computation. The advan- 
tages of expressing the basic rates in the form of a matrix are more clearly seen in con- 
sidering the type of problem such as the following. Let us suppose that a species of mammal 
at a certain season of the year invades a fresh environment where there is an ample food 
supply, a freedom from predators, and plenty of space to accommodate any rapid increase 
in numbers which might take place. Under these conditions it might be assumed for theo- 
retical purposes that some age specific rates cf fertility and mortality would remain approxi- 
mately constant over a period of time. The age distribution of these immigrants then 
becomes of some importance owing to the effect which it must necessarily have on the future 
course of events. For this initial distribution must clearly be very different from that which 
would ultimately be established in the case of a species, such as a rodent, with possibly a 
very rapid rate of increase, since nestlings will not be represented in it and young individuals 
may be present in only relatively small numbers. Supposing then that we have a number 
of such populations subject to the same age schedules of fertility and mortality, but differing 
in the age distribution of the original immigrants, we may have reason to enquire how far 
the development of ti.sse populations is affected over a limited period of time by the varying 
form of this initial distribution, assuming for simplicity that no further waves of immigration 
occur. 

If an estimate of the number and age distribution of the female population at successive 
intervals of time is alone required, the answer for any form of initial distribution is readily 
obtained once the series of matrices M, M?, M°, ..., M' have been constructed. But, in 
addition, we may require to know the changes which might be expected to occur in the birth 
rate and death rate, and also, for example, in some such rate as the percentage of adult 
females pregnant, a figure which is one of the simplest measures we have of the degree of 
fertility among wild populations. Again, in a species like the wild rat we never know the 
exact age of individuals caught in the field, and thus the only measure of the form of the 
female age distribution is the percentage of immature females, provided, of course, we are 
sampling the complete population. Some method is therefore required for calculating such 
rates at successive intervals of time. 

Once the age distribution of the female population at time ¢ is known, an estimate of the 
expected number of female births per unit of time may be obtained by operating on the age 
distribution with the maternal frequency figures. Thus, in matrix notation we may write, 
the number of female births equals m,M'E,, where &, is the initial age distribution and the 
m, figures are treated as a row vector. Similarly the estimated number of deaths per unit of 
time may be obtained with the help of the age specific death rates (D,). The relative rate of 
increase calculated in this way is not necessarily exact, but it may be sufficiently accurate 
for our present purposes. As an example of the degree of error involved in this method, we 
may compare the values of the stable birth rate and death rate, as given in the appendix, 
with those derived from the matrix stable distribution in Table 2 by operating with the m, 
and D, figures. The latter were in this case computed from the stationary age distribution 
and the d, column of the life table (D, = d,/L,). The results were as follows: 
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By operating with 
‘True’ values m, and D, on 
from appendix matrix stable 
distribution 
Birth-rate (b) 0-51265 0-51257 
Death-rate (d) 0-06700 0-06154 
b-—d=r 0-44565 0-45103 

















The rate of increase is overestimated by about 5-4 per thousand per month, the principal 
error being in the death rate. This discrepancy is due to the fact that the number of deaths 
under 2 months of age is underestimated by applying age specific death rates, which are 
based on the stationary age distribution, to the stable population grouped in one month 
intervals at these ages. The difference between these distributions happens to be quite 
marked in this example. The degree of error, however, is not very great; and in the type of 
problem we are considering, when the age distribution of a population may take any form, 
this seems to be the only practical method of estimating the rate of increase.* 

Supposing, then that in the case of the rat population used here as a numerical illustration, 
we wished to estimate the number of females, the birth rate and death rate, and the per- 
centage of immature females at monthly intervals up to—say—7 months from the origin 
of the time scale, when the initial immigration is assumed to take place. Since the jth column 
of the matrix gives the age distribution of the survivors and the surviving descendants per 
individual female alive in the age group j—1 to j at ¢ = 0, the sum of the elements in this 
column gives the number of times the original population in this age group has increased, 
or decreased, at time t. The percentage of immature females may be obtained once the sum 
of the first three elements in the column is known (reproduction begins at the age of 3 months 
in this example); and the number of births and deaths per unit of time may be found by 
operating on the column with the m, and D, figures. Each of these totals, of course, will have 
to be multiplied in the end by the number of females alive in this age group at ¢ = 0. Since 
the initial age distribution may be of any form under the conditions of the problem, it will 
be necessary first of all to calculate these four totals for every column of each of the seven 
matrices M'. Now, to add up the elements forming each column of a matrix is equivalent to 
premultiplying the matrix by a row vector of units; the sum of the first three elements may 
be obtained by premultiplying with a row vector of which the first three elements are units 
and the remainder zeros; and similarly the numbers of births and deaths are found with the 
help of the row vectors m, and D,. Thus, the operations which it is necessary to perform on 

* Another similar method, which avoids the actual calculation of the number of deaths by means of the age 
specific death rates, is suggested by the following relationship. If the transformed age distribution y = H&, 
where H is the matrix defined in § 5, is operated on with a row vector which consists of the L,m, figures (Ap- 
pendix, Table 4), and the resulting scalar is divided by the sum of the elements of y, an estimate is obtained 


of the relative rate of increase of a population with an age distribution & in the original co-ordinate syst: m. 
This follows from the properties of the transformed population discussed at the end of § 5 and from the relation- 


ship between the first row of the canonical form B = HAH- and the L,m, column (§ 6). In the transformed ° 


population the death rate =0, and the maternal frequency is given by L,m,. Thus, by transforming back again 
to the original co-ordinate system 

_(L,m,] HE 

~ (tHe 


where [1] defines a row vector of units. Taking £ as the matrix stable distribution of Table 2, and calculating 
the row vectors [L,m,] H and [1] H, the value of r was estimated by this method to be 0-44468. 
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each of the columns of M‘ may be written as the matrix R, which will consist of m + 1 columns 
and » rows, the number of the latter depending on the number of operations. Then the 
required totals for the matrix M‘ will be given by 


Z = RM‘ = RMMM...M, 


and it is easy to see how the Z matrices may be built up in succession withont calculating the 
actual matrices M*, M?, M4, etc. Once the series of Z matrices have been constructed, we 
can obtain from Z€, the necessary figures from which the required rates at time ¢ for a 
population with an initial age distribution £ may be calculated. Moreover, if we wish, the 
contributions made, for instance, to the total number of births or deaths by any particular 
age group in the initial distribution can also be determined. 

The computations in this illustration have been greatly simplified by the assumption that 
the system of age specific fertility and mortality rates remains constant. In the case when 
the basic matrix M is ehanging with time and the age distribution at time ¢ is given by 
M,... M;M,M, £,, some of the rows of R will also be varying. Hence the series of Z matrices 
could not be built up without first computing M,M,, M,M,M,, etc. The latter, however, 
may often be of interest in themselves. For, if each column of M‘—or M,... M,M,M, in the 
case of a variable matrix—is multiplied by the number of females alive in the appropriate 
age group at ¢ = 0, the complete age structure of the population at time ¢ is represented in 
the form of a two-dimensional array. Since the sum of the elements in each row is the total 
number alive in the age group x to x+ 1 at time t, the number contributed to this total by 
each age group at ¢t = 0 is given by the individual entries. 


APPENDIX 
(1) The tables of mortality and fertility 

The basic life table and fertility table which have been used in the numerical part of this 
study are given in Table 4. The adult ,, figures from the age of 2 months onwards are based 
on the mortality observed among the females of a domesticated brown rat stock housed at 
the Wistar Institute, Philadelphia. According to the data for 26 generations of this laboratory 
stock published by King (1939) it appears that out of 1384 females alive at the age of 2 months 
(60 days), 1337 were alive at 12 months, and 984 at 20 months. This information gave three 
points on the 1, curve, supposing that these survivors could be regarded as ordinates at these 
exact ages. In order to interpolate for other ages, a logistic type of curve was fitted to the 
data, the values of the constants being chosen so that the curve passed through these three 
points. The 1, values in Table 4 are given by 

0-85156355 
== T+ o-00101065erI0* fF 722. 

Although the original data did not extend beyond the age of 20 months, by which time the 
vast majority of the females had ceased breeding (King, 1939), this J, curve was extrapolated 
to later ages, whenever necessary, simply for the purposes of this theoretical study. 

The degree of infant mortality assumed here, namely 15 % between birth and the age of 
2 months, is entirely arbitrary; it represents a moderate degree of loss at these early ages. 
Some care, however, was taken to weld the infant mortality smoothly on to the remainder 
of the 1, curve, and it was assumed that the number of deaths according to age (d,,) decreased 
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geometrically between birth and the age of 2 months. The actual calculations for these age 
groups were carried out in units of 1/8 of a month and the resulting J, curve was integrated 
by means of Simpson’s rule. The same method of numerical integration was used also for 
the adult part of the life table in order to obtain the L, figures. 

The fertility table is partly artificial and was constructed in the following way. The gross 
reproduction rate of these domesticated brown rats was estimated from the data published 
by King (1939) to have been just under 10 litters for the later generations, when the stock 
was thoroughly adapted to life in the laboratory. The frequency of litter production according 
to the age of the mother has been found by the author (unpublished observations) to be 
represented closely by a Pearsonian type I curve im the case of certain litter fertility tables, 
for example in the cross-albino rat, the vole and some human populations with a high degree 
of fertility; and, moreover, the values of £, and f, were very similar for all three species. 
The actual equation for the curve used here may be written 


y = y'2*(q— 2), 

and the range was assumed to be from 3 to 21 months, which for grouping purposes represents 
the span of the reproductive ages observed by King in this Wistar strain of brown rats. The 
ordinates of the integral curve of the above equation were taken from the-tables of the in- 
complete Beta function and, in this way, a column was formed which gave a gross reproduc- 
tion rate of 10 litters. The individual entries were then multiplied by the mean number of 
daughters per litter according to the age of the mother, which was recorded by King, and thus 
the m, figures in Table 4 were obtained. The gross reproduction rate is 31-21 daughters and 
the net rate 25-66. These tables of fertility and mortality were originally constructed in 
order to determine the relative rate of increase and the type of stable age distribution which 
might be expected in a brown rat population living under more or less optimum conditions. 


(2) Calculation of the rate of increase 


Some difficulty was experienced in obtaining a satisfactory estimate of the rate of increase 
(r) from the usual solution (Dublin & Lotka, 1925) of the equation: 


ao 
} e=1_m,dz = 1. 
0 


The 4th degree equation in r with the numerical coefficients based on the seminvariants of 
the L,m, distribution, the estimates of which are given in § 6, was 


4-90199r4 + 3-69283r3 — 7-07198r? + 9-60604r — 3-2448498 = 0, 


and the real root was found to be 0-42447. This value of r was, however, clearly too low and 
a better estimate had to be obtained in a rather roundabout way, since it was thought that 
the use of higher moments than the fourth would be unsatisfactory in the present example. 

If the force of mortality represented by the original life table is increased by a constant 
factor (r) which is independent of age, the new life table is 


i = eI, 


and the net reproduction rate will be given by R, = 2Lim,, where L;, are the integrals of 


the new J, curve. Clearly, the greater r is taken to be, the smaller R, becomes. Then, suppose 
that the relation between R, and r is given by 


log, R, = a+ br+cr*+dr>+..., 
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Table 4 
Age (x) Life table Fertility table 
Units of or, 
30 days 
1, L, m, L,m, 
0 1-00000 0-46544 _ — 
0-5 0-88706 0-43489 — — 
1 0-85882 0-42725 - oe 
1-5 0-85176 0-42534 — — 
2 0-85000 0-84973 a -— 
3 0-84945 0-84910 1-1342 0-96305 
4 0-84871 0-84824 2-0797 1-76408 
5 0-84772 0-84708 2-6596 2-25289 
6 0-84638 0-84553 2-8690 2-42582 
7 0-84458 0-84344 2-9692 2-50434 
| 8 0-84217 0-84063 2-9535 2-48280 
9 0-83893 0-83687 2-8143 2-35520 
10 0-83459 0-83184 2-6114 2-17227 
ll 0-82881 0-82515 2-2455 1-85287 
12 0-82113 0-81629 2-0533 1-67609 
13 0-81098 0-80463 17971 1-44600 
} 14 0-79768 0-78940 1-5561 1-22839 
15 0-78039 0°76975 1-2175 0-93717 
16 0-75821 0-74472 0-9548 0-71106 
17 0-73018 0-71342 0-6610 - 0-47157 
} 18 0-69548 0-67512 0-4043 0-27295 
19 0-65354 0-62953 0-1846 0-11621 
| 20 0-60434 0-57696 0-0435 0-02510 
21 0-54859 oa —- — 
Lee, 
Total. — — 31-2086 25-65786 
where a = log, Ry, and the constants ), c, d, etc. are to be determined. It was assumed in the 


present instance that a 4th degree polynomial in r would be sufficient, and four new life 
tables were constructed taking r to be in turn 0-1, 0-2, 0-4 and 0-5. The Zi, integrais were 


| obtained by Simpson’s rule for the reproductive ages and the four values of R, calculated. 
The equations for finding the values of the constants were: 


0-0001le + 0-001d + 0-01c + 0-16 = — 0-8934974, 
0-0016e + 0-008d + 0-04c + 0-26 = — 1-6708610, 
0-0256e + 0-064d + 0-16c + 0-46 = — 2-9792984, 
0-0625e + 0-125d + 0-25c + 0-56 = — 3-5490542, 


whence b = —9-617235, 
c= 17-371816, 
d = —5-698291, 
e= 2-062332. 


Inserting these values of the constants in the equation for log, R,, the value of r for which 
log, R, = 0 was found to be 0-44565, 


The stable birth rate was estimated in the usual way from 


1 fe 
57 foe 1dz. 
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The integrals were computed by means of Simpson’s rule, treating the age groups 0-2 
separately from the remainder of the life table, the units adopted being 1/4 month for the 
early ages compared with 1 month for the 























later. It was found that Table 5 
21 ——_—_——-————_ 
—0°445652, = 1-0506 
I, e l,dx = 1-95064, Matrix A Matrix B 
b = 0-51265, Age 
group 
and hence the death rate, x P PF Elements of 
ad ” first row 
d = 0-06700. 
The value of r is so high in this case that 0- 0-94697 0 0 
the error in the estimate of b due to 1- 0-99665 0 0 
: 9 2- 0-99926 0-3964 0-3741 
neglecting the ages from 21 onwards would 3_ 0-99899 1.4939 1-4089 
only be in the last figure. The stable age 4- 0-99863 2-1777 2-0517 
a" . : , ; . 5- 0-99817 2-5250 2-3756 
distribution is given in Tables i and 2 of 6 099753 92-6282 92-4682 
the text. Owing to the way in which the 7- 0-99667 2-6749 2-5059 
: . . 8- 0-99553 2-6018 2-4293 
value of r was determined, it will be found 9- 0-99399 92-4419 92-2698 
that the birth rate of the stable popula- ® = 2-1865 2-0202 
. . : ° - 0-98926 1-9044 1-7454 
tion obtained by operating on this age 12- 0-98572 17259 15648 
distribution with the maternal frequency 13- 0-98107 1-4918 1-3332 
fi A; rs that oi 14— 0-97511 1-2415 1-0885 
gures 1s precisely the same as that given 15- 0-96748 0-9522 0-8141 
by the above integral. 16- 0-95797 0-7141 0-5907 
17- 0-94631 0-4618 0-3659 
. . 18- 0-93247 0-2518 0-1888 
(3) N umerioal values of the matrix elements f 19- 0-91649 0.0901 0.0630 
The numerical elements of the matrices 20- 0-0035 0-0022 
A and B to which reference has been made 








in §§3 and 6 of the text are given in Table 5. 
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ERRORS (N THE ROUTINE DAILY MEASUREMENT OF 
THE PUERPERAL UTERUS 


By C. SCOTT RUSSELL, M.A., M.B., Cu.B., F.R.C.S.Ep., M.R.C.O.G. 
First Assistant, Nuffield Department of Obstetrics and Gynaecology, University of Oxford 


In both district and institutional midwifery, it is usual to measure each day the height of 
the puerperal uterus above the symphysis pubis. Though most observers would say that the 
purpose behind this custom is the detection of uterine subinvolution, there is no agreement 
as to the precise conditions that must be fulfilled or the criteria required for establishing 
such a diagnosis. The present paper is concerned with the systematic and the uncontrolled 
errors of this routine measurement, and with the significance of fluctuations in the observed 
height of the puerperal uterus. 

In institutions, the puerperal uterus is usually measured by whichever member of the 
nursing staff happens to be available: often the measurements are made by pupil midwives. 
It must be rare for a uterus in any period of 7 days to be measured by only one observer. 
The measurement, which is made between the upper surface of the uterus and the upper 
surface of the symphysis pubis, is recorded on the temperature chart; it is probably never 
recorded more accurately than to the nearest } in., and usually the nearest } in. is considered 
sufficiently accurate. 

NORMAL UTERINE INVOLUTION 

It is widely taught that the normal human uterus involutes }in.a day during the puerperium. 
This is only an approximation; if it were true, the puerperal uterus having shrunk in 1 day 
the last } in., would suddenly stop involuting. This would be contrary to clinical experience 
which tells us that the process starts fast and finishes slowly. The changing daily amount 
of normal involution would be too cumbersome a standard in practice; besides, great accuracy 
is not necessary when the uterus is measured only once a day and slight irregularities are of 
little clinical importance. A rate of } in. a day may therefore be taken as accurate enough 
for practical purposes. 


FACTORS INFLUENCING THE RECORDED HEIGHT OF THE UTERUS 


The factors that influence the recorded height of the puerperal uterus may be divided into 
two broad groups. First, there are those which act constantly in one direction, for example, 
the distended bladder or lower bowel. Secondly, there are the uncontrolled factors, such as 
the thickness of the abdominal wall, the mobility of the uterus, the uterine tone, the method 
of measurement, the points chosen for measuring, the individual differences between 
observers and so on. Though these many factors are just as likely to cause the reading to be 
too high as too low, they are important because together they constitute the error of the 
measurement; without knowledge of this error, logical deductions based on the measure- 
ments are not possible. The following five studies were made to obtain estimates of this error. 


METHOD 
First investigation 
Seven puerperal patients whose uteri could be felt were chosen, and six unpractised 
observers measured each uterus twice, the interval between measurements being about 
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15 min. For convenience, medical students were used as observers; it was thought that they 
would be neither more nor less accurate than pupil midwives. An effort was made to prevent 
conscious or subconscious bias on the part of each observer (1) by explaining the nature of 


Table 1. Duplicate measurements of puerperal uterus by unpractised observers 







































































Observers’ ——— Observers’ 
dupli if ers 
luplicate totals 
readings 1 2 3 4 5 6 7 
1 5-00 5-00 4-25 4-00 4-25 5°75 4-75 33-00 68-25 
5-50 4-25 4-75 4-50 4-75 6-00 5-50 35-25 
2 4-75 4-75 4-25 4-75 4-50 5:50 5-00 33-50 65-00 
5-00 4-25 3°75 4-00 4-25 5-25 5-00 31-50 
3 5-50 4-75 4-00 4-25 4°75 5°75 5-00 34-00 69-50 
5°75 4°75 4-50 4-50 4-75 5°75 5-50 35-50 
4 6-00 5-00 4-50 3°75 5°50 5-50 5-50 35°75 76-25 
6-75 5°75 5-25 4-50 6-50 6-00 5-75 40-50 
5 5-25 450 4-00 4-00 5-75 6-00 4-50 34-00 70-25 
5-75 4°75 5-00 4-50 5-25 5°75 5°25 36-25 
6 5-25 4-50 4-25 4-25 6-00 5°50 5-25 35-00 72-00 
6-00 5-00 5-00 4-75 5-75 575 | 475 37-00 
31-75 28-50 25-25 25-00 30-75 34-00 30-00 205-25 
oer 28-75 28°25 26°75 31-25 34-50 31-75 216-00 
Patients 
totals |— 
66-50 57°25 53°50 51-75 62-00 68-50 61-75 421-25 
Analysis of variance 
Sum of |Degreesof} Mean Variance a 
squares | freedom | square ratio Probability 
Main effects 
Patients 19-9349 6 3-3225 43-6310 Less than 0-001 
Observers 5 1-0169 13-3539 Less than 0-001 
1-3778 1 1-3778 18-0932 Less than 0-001 
First order interactions 
Patients-observers 5-6187 30 0-1873 2-4576 0-05-0-01 
Patients-readings 0-6798 6 0-1133 1-4879 Greater than 0-20 
Observers-readings 1-6891 5 0-3378 4-4360 Less than 0-01 
Second order interactions 
Patients-observers- 2-2845 30 0-0762 
readings [error] 
Total 36-6693 83 














the investigation, and (2) by ensuring that previous readings of other observers were not 


known until subsequent readings were recorded. Measurements to the nearest } in. were 
requested. In Table 1 the measurements are recorded together with the corresponding 
analysis of variance. This analysis and certain important facts arising from it must be 
considered in detail. 
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Analysis of variance 

The total variability of the eighty-four separate observations may be expressed by a 
“mean square’, i.e. the sum of squares of deviations of each observation from the arithmetic 
mean divided by the degrees of freedom. This variability, as expressed by a ‘mean square’, 
is divisible into mean squares due to the main effects and the first order interactions together 
with a remainder which is the best first estimate of uncontrolled error. Based on 30 degrees 
of freedom this has a mean square of 0-076, against which all the other mean squares are 
compared to give the different variance ratios. From these the probabilities are found from 
the appropriate tables (Fisher & Yates, 1943). 

The interaction ‘patients-readings’ as expressed by the corresponding variance ratio is 
not significant;* the interval between readings, therefore, appeared to have had a similar 
effect on the seven patients. The interaction ‘observers-readings’ is highly significant.* 
This is found, by an analysis of variance on the first and second readings taken separately, 
to be entirely due to the second readings. Whereas the variance ratio ‘between observers’ 
obtained from the first readings is not significant (0-933), that from the second readings is 
highly significant (11-03). A possible explanation is that the observers largely responsible 
subconsciously altered the technique of measurement between reading 1 and reading 2. The 
interaction ‘patients-observers’ is significant,* suggesting that the observers found the 
different uteri difficult to measure in different ways. j 

Passing to the main effects, the variance ratio ‘patients’ is highly significant, as would 
be expected from patients observed on different days of the puerperium. The variance ratio 
‘observers’ is highly significant, emphasizing the importance of the personal factor. The 
variance ratio ‘readings’ is also highly significant; this may be due to an accumuiation of 
urine between the first and second readings causing slight upward displacement of the 
uterus. 

The analysis of variance has been taken to the limit. Of the three first order interactions, 
the first and third are not controllable with unpractised observers; the second is not sig- 
nificant. All three may therefore be included in the estimate of uncontrolled error. The mean 
square of the adjusted error is now 0-1447. The variance of the difference between two 
observations is therefore 0-2894 and the standard error of the difference, the square root of 
this, i.e. 0-5380. If the effect of different observers is taken into account, the mean square 
‘observers’ must be included in the error estimate. The mean square for error is then 0-2021, 
and the standard error of the difference between two readings 0-6358. 


Second investigation 

A similar study was carried out with seven other unpractised observers and seven other 
patients. In addition to the duplicate readings being separated by an interval of time— 
approximately 15 min.—each patients was given 0-5 mg. ergometrine by mouth as soon as 
the first readings had been completed. The results are given in Table 2. All the main effects 
and all the first order interactions are highly significant, supporting the evidence already 
presented that the measurement not only is difficult but also is influenced by many factors. 
The three first order interactions would not be controlled in practice, so that all may reason- 
ably be included in the error estimate. This addition alters the mean square for error to 
0-1621. The standard error of the difference between two observations is therefore 0-5694, 

* The term ‘significant’ is used where the probability is between0-05 and 0-01. The term ‘highly significant’ 


is used where the probubility is less than 0-01. ‘Not significant’ refers to a probability greater than 0-05. 
Biometrika 33 1+ 
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which is in close agreement with the value 0-5380 obtained from the first investigation. The 
inclusion of the effect ‘between observers’ raises the standard error of the difference to 
0-6821, which also agrees closely with the corresponding figure 0-6358 obtained from the 
first investigation. 


Table 2. Duplicate measurements by unpractised observers before and after ergometrine 




































































Observers’ Patients 
before | Observers’ 
and after totals 
ergometrine 1 2 3 4 5 6 7 
1 5-25 5-50 3-25 4-75 6-25 | 6-50 | 4-50 36-00 68-00 
4-75 4-50 2-75 3-50 5-50 6-00 5-00 32-00 
2 6-00 5-50 3-00 4-25 6-50 6-50 5-00 36-75 72-25 
5-25 5-50 3°00 4-50. 6-00 6-00 5-25 35-50 
3 4°50 4-75 3-00 4°50 6-00 5°75 4-75 33-25 63-50 
4-50 4-50 2°75 4-25 5-50 4-50 4-25 30-25 
4 4-50 5-50 4-25 5-00 5°75 5-50 4-75 35-25 66-00 
4-50 5-00 3-00 4°25 5-00 450 | 4:50 30-75 
5 5-00 _ 550 2-75 4°75 6-00 6-50 | 5°50 36-00 72:75 
5-00 5:50 3-00 5-00 650 6-00 | 5-75 36-75 
6 4-50 5°25 2-75 4-50 6-00 575 | 425° 33-00 64-75 
4°75 5-00 2-25 4:75 5°75 4575 | 450 31-75 
7 5-50 5-25 3-50 4-00 6°25 625 | 525 | 36-00 | 73°50 
5°75 5°75 3°50 5-00 6-25 575 | 5:50 37-50 
| . } 
35°25 37-25 22-50 31-75 42-75 42°75 34°00 | 246-25 
34-50 35°75 20°25 31-25 40°50 37-50 34°75 | 234-50 
Patients’ 
totals 
69°75 73-00 42-75 63-00 83°25 80-25 68-75 480-75 
| 
Analysis of variance 
| | | | | 
|  Sumof | Degreesof | Mean | Variance | = 
| | squares | freedom | square | ratio Probability 
| | 
| | Fried sony | 7 
Between patients | 76-4707 6 | 12-7451 | 182-0707 | Less than 0-001 
Between observers | 73190 | 6 12198 =| 17-4257 Less than 0-001 
Between readings 1-4089 1 | 11-4089 20-1271 Less than 0-001 
Patients-observers | 7-3061 | 36 | 02029 | 2-8986 | Less than 0-01 
Patients-readings 1-5420 | 6 | 02570 | 3-6714 Less than 0-01 
| Observers-readings |  2-2473 6 | 03745 | 53500 | Less than 0-001 
| Exror | 25205 | 36 | 00700 | 
| 
| | | | 
| Totals | 98-8145 | 97 | 
| 








Error of practised observers (third investigation) 
The error of measurement with unpractised observers was of special importance for those 
institutions where the routine puerperal uterine measurements are made by pupil midwives. 
Elsewhere the routine is the responsibility of comparatively senior nurses. In order to 
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obtain an estimate of the error of trained observers another simple study was made. Ten 
puerperal patients each had the uterus measured by five observers of whom four were nurses 
and one was a doctor. The observers were experienced in the measurement; they were asked 
to use the method to which they were accustomed. Each observer measured each uterus 
once. In Table 3 the measurementsare tabulated with the corresponding analysis of variance. 

It should be noted that not one of the measurements was given more accurately than to 
the nearest } in., and that forty-four out of the fifty readings were given to the nearest } in. 

The two main effects, ‘between patients’ and ‘between observers’, are both highly signi- 
ficant, as would be expected. The remainder or error estimate has a mean square of 0-3308, 
which is consideravly higher than that of the unpractised observers. This apparent paradox 


Table 3. Measurements of the puerperal uterus by practised observers 




















Patients 
Observers Observers’ 
totals 
1 2 3 4 5 6 7 8 9 10 

7 | | 
1 5-00 | 500 | 500 | 300 | 400 | 450 | 250 | 400 | 400 | 625 | 43-25 
2 425 | 400 | 575 | 150 | 250 | 300 | 100 | 250 | 3-00 | 5-00 32-50 
3 450 | 450 | 500 | 200 | 250 | 200 | 300 | 1:50 | 400 | 600 | 35-00 
4 450 | 5:00 | 600 | 250 | 450 | 3:00 | 1:50 | 4.00-| 3:50 | 5-50 40-00 
5 450 | 5:00 | 550 | 300 | 400 | 4:25 | 385 | 3:50 | 425 | 650 | 42-75 

| | | 

Pationta’ 22-75 23:50 27:25 12:00 17:50 1675 1025 1550 1875 2925 | 193-50 

tota 

















Analysis of variance 





























Sum of Degrees of Mean Variance 1 
squares freedom square ratio Seely 
Between patients 70-2800 9 7-8088 23-608 Less than 0-001 
Between observers 9-0925 4 2-2731 6-87 Less than 0-001 
Error 11-9075 36 0-3308 
Total 91-2800 | 49 











is probably largely due to the fact that the unpractised observers were asked to measure to 
the nearest } in. and thereby encouraged to be accurate, whereas four out of five of the 


practised observers made the measurement to the nearest } in.—except on three occasions 
—as this was their usual practice. 


Error of individual observers (investigations 4 and 5) 
Two studies were now made to obtain estimates of the error of individual persons. 
An experienced midwifery sister tutor made five replicate readings to the nearest } in. 
on nine puerperal patients. The results are shown in Table 4. The mean square between 
readings was not significant and was therefore included in the error estimate. The mean 


square for error was now 0-1639, which is not appreciably different from 0-1447 and 0-1621 
which are the corresponding figures previously obtained. 
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Having had special experience of the measurement, I repeated the above study. Of twelve 
unselected puerperal patients only eight were measured, the remaining four being thought 
unsuitable because of abdominal tenderness (one case), recent delivery (one case), and 
difficulty in palpating the uterus (two cases). This partial selection of cases left for analysis 
only those whose uteri could be felt fairly easily. The results are shown in Table 5. The mean 
square between readings was again not significant and was included in the estimate of error. 
The new value was 0-0813, which is lower than the previous best. The improvement, though 
partly the result of greater care in measurement, is probably largely due to some selection 
of material. 


Table 4. Replicate measurements by midwifery sister tutor 





























Patients 
Replicate Readings’ 
readings totals 
1 2 3 4 5 6 7 8 9 
1 5-00 5°75 3-00 2-00 3-75 3-00 4-00 4-50 6-25 37-25 
2 5-00 6-00 2-75 2-75 3-75 3-00 5-00 4-50 6-00 38-75 
3 5-00 5-00 3-00 3-00 4-75 3-00 4-75 4-50 6-00 39-00 
4 4-75 6-00 3-50 4-00 4-75 3-25 4-50 4-25 6-00 41-00 
5 4-50 5-00 3-00 | 3-75 4-50 3-00 4:75 4-50 6-00 39-00 
Patients’ | 24-25 27-75 1525 15:50 21-50 15:25 23:00 22-25 30-25 195-00 
totals 

















Analysis of variance 























Sum of Degrees of Mean | Variance oi 
squares freedom square ratio Probability 
Between patients 47-9750 8 5-9969 | 37-565 Less than 0-001 
Between readings 0-7916 4 0-1979 1-240 Greater than 0-20 
Error 5-1084 32 0-1596 
Readings + error 5-9000 36 0-1639 
Totals 53-8750 44 

















EFFECT OF THE BLADDER ON THE PUERPERAL UTERUS 


Attention must now be directed to the factors which bias the recorded height of the uterus. 
The only factor which constantly makes the recorded size smaller than reality is a tendency 
for the uterus to fall backward into the pelvis towards the end of the first week of the puer- 
perium. There are, however, two important influences tending to make the uterus seem larger 
than reality. The one is the urinary bladder, of special importance in the immediate post- 
partum period: the other is the pelvic colon and rectum, of more importance a few days later 
The effect of the bladder is well known, and I have even found the uterine fundus displaced 
so far upwards by the bladder that it could be felt under the costal margin and ballotted 


like the foetal head in a breech presentation. To show this upwards displacement the 
following short study was made. 
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Table 5. Replicate measurements by author 
































Patients 
Replicate __| Readings’ 
observations totals 
1 2 3 4 5 6 " 8 
1 4-75 | 350 | 250 3-00 5-00 | 525 | 7-00 4-25 35-25 
2 425 | 3-25 | 2-00 3-00 5-25 5-25 6-50 4-75 
3 4-50 3-50 | 2-50 3-00 5-50 5-25 6-00 4-50 34-75 
4 4-50 3-50 | 2-00 3-25 5-00 4-75 7-00 4-50 34-£0 
5 3-75 3-25 | 2-25 3-50 | 475 | 4-75 | 6-50 4-75 50 
38 | 
Patients’ 21-75 17-00 11-25 15-75 25-50 25-25 33-00 22-75 | 172-25 
totals 





ast of variance 

















7 ] 
| Sumof | De f | | 
grees 0 Mean | Variance | F 
| squares c freedom square | ratio Probability 
eresene | |. | 
Between patients | 64-4609 7 92087 | 107:83 | Less than 0-001 
Between readings | 0-2093 4 0-0523 
Patients-readings [error] | -* — 28 0-0854 | 
Readings + error ites 32 0-0813 
Totals 67-0609 39 

















Fourteen unselected patients were used. Each patient’s uterus was measured by myself 
just before a morning bed-pan round. Following the use of the bed pan each uterus was 
again measured. In Table 6 the duplicate readings and the differences are tabulated. Being 
concerned with the displacement of the uterus by the bladder the appropriate statistical 
test is based on the difference column, in which it should be noted all fourteen figures bear 
the same sign. Using the ¢ test (Fisher, 1941, p. 117), the probability is found to be less than 


Table 6. Measurements of the puerperal uterus before and after micturition 









































Measurements Measurements 
Difference Difference 
[y] ‘ [y] 

First Second First Second 
7-50 6-25 1:25 6-25 5-25 1-00 
4-25 3-75 0-50 5-00 2-00 3-00 
7-00 4-50 2-50 6-75 5-00 1-75 
6-00 | 4-50 1-50 7-50 5-25 | 2-25 
8-00 4-25 3-75 425 — 2-50 1-75 
3°75 2-25 1-50 6-50 4-25 2-25 

00 4:00 1-00 5-50 4-00 | 1-50 | 

| 
— #2 
Sum oat 1 = 0-7445 ores = 0-0532 = (—0-231)*. 


g=1-8214 t,=783 P<0-001. 
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0-001. As the condition of the bladder was the only important factor that varied between 
reading | and reading 2, there can be no doubt that the uterus had been displaced upwards. 

Because one would expect a high correlation between the displacement of the uterus and 
the amount of urine passed, each difference was compared with the corresponding output 
of urine. The regression coefficient of fall in uterine height on urinary output is 0-046. The 
mean square for error is very high—further evidence of the difficulty in measuring the 
puerperal uterus—and the regression coefficient is, in fact, not significantly different from 
zero (P between 0-1 and 0-2). If further cases were studied, one would expect from clinical 
experience that a significant coefficient would be obtained. 


Discussion 


Though we cannot tell in advance precisely how much a given uterus will involute in a given 
time, it is known that the average involution rate of the healthy organ is, as a close approxi- 
mation, }in. per day during the early days of the puerperium. Therefore, the expected height 
of a uterus on any one of these days may be estimated as } in. less than that recorded for the 
previous day. 

The records of the uterine height are influenced by many factors outside clinical control; 
these are of real importance because together they constitute the error of measurement. 
Without knowledge of the magnitude of this error, it is impossible to judge the significance 
of measurements and differences between measurements. 

Five investigations have been made to determine the magnitude of this error; two were 
with unpractised observers and three with practised observers. The standard error of the 
difference between two readings of a single average careful observer is approximately } in. 
Unless, therefore, the difference between the recorded height of the uterus and the expected 
height is more than 1 in., i.e. twice the standard error of the difference, it is not safe to 
conclude that anything other than the error of the measurement has been responsible.* But, 
according to our estimate, the uterus should involute only } in. per day. Therefore, unless 
the recorded height is at least } in. greater than that recorded for the previous day, we should 
not be surprised nor assume anything to be wrong, as such a difference is likely to be due to 
this error only. This argument may be extended to cover comparisons between measurements 
at intervals up to 2 or 3 days. If taken much beyond 3 days, however, a serious fallacy might 
result. Our estimate of the average daily uterine involution is an approximation, and as this 
is multiplied, so will the error of the approximation increase. 

By exercising great care, measuring to the nearest } in., and selecting only those patients 
whose uterus can be felt easily, it should be possible to reduce the standard error of the 
difference between two measurements to the region of } in. Any investigation designed to 
compare involution rates under different conditions would be helped by such increase in 
accuracy. 

There are two common controllable influences which elevate the puerperal uterus. The 
one, more important in the early days of the puerperium, is a full bladder; the other, 
exerting its influence later, and of !ess importance, is a full bowel (Moir & Russell, 1943). 
It has now been shown, not only that the full bladder raises the uterus, but also that in 


* Strictly the comparison between normal and abnormal uterine involution involves a test of significance 
between two regression lines. In practice, however, it.would be unreasonable to expect any medical attendant 
to calculate the regression line on every occasion that the size of the uterus was larger than expectation. The 
method described should be adequate and, requiring only mental arithmetic, is suitable for the bedside. 
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thirteen oui of fourteen patients with normal bladder distension the extent was up to or 
greater than twice the standard error of the difference between two measurements of 
uterine height. With pathological bladder distension, the upward displacement of the 
uterus will be still greater and correspondingly more than the error of measurement; in 
contrast, alteration in uterine height due to subinvolution will be less than the error of 
measurement. 

The practical conclusion is that careful routine daily puerperal uterine measurement is 
of value, more because it helps in the recognition of abnormal states of the bladder or bowel 
than those of the uterus. Distension of the bladder or bowel should always be suspected if, 
from any two measurements over a period of not more than 3 or 4 days, the observed height 
of the uterus is more than | in. greater than the expected height. It must, however, be 
emphasized that this figure of 1 in. has been obtained from measurements carefully made to 
the nearest } in.; it must be increased to 14 or even 2 in. if measurements are perfunctorily 
made by different persons each day. A diagnosis of uterine subinvolution is only permissible 
after due allowance has been made for the effect of all factors, controllable and uncon- 
trollable, that influence the recorded height. 


I wish to record my thanks for help in the preparation of this paper to Mr P. H. Leslie, 
Dr W. T. Russell and Professor Chassar Moir. 
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ON A METHOD OF ESTIMATING FREQUENCIES 
By J. B. 8S. HALDANE, F.RS. 


In a very large variety of investigations we desire to estimate the frequency of an attribute 
in a series of populations, each of which is so much larger than the sample taken from it that 
it may be regarded as infinite. If p be the frequency of the attribute, q = 1—>p, and n the 
number in the sample, the standard error of the estimate of pis of course ,/(pq/n). Provided. 
p does not vary much from one sample to another, it is desirable to keep n approximately 
constant. But when p varies greatly this is unsatisfactory. Thus if n = 1000, and p = 0°3, 
its standard error is 0-015, but if p is 0-01 its standard error is 0-0031, so that we could not 
distinguish between populations where p was 0-01 and 0-005. In such a case it may be desired 
that the standard error of each value of p should be roughly proportional to p rather than to 
its square root. 

Such cases arise in haematology. My friend Dr R. A. M. Case has been investigating the 
frequency of siderocytes, an abnormal type of red blood corpuscle described by Griineberg 
(1942), in a number of bloods. Their frequency ranged from about 0-2 to 20 % or more. He 
has adopted the method of counting stained red corpuscles in a film until he had counted 
some definite number m, usually about 20, of siderocytes. If, in order to count this number, 


he had observed a total of n red = he took mn as his estimate of p. It will be 


shown that the correct estimate is ——— = , So that his error was negligible; and the standard 
error will be calculated. feake 
Let p be the frequency of abnormal cells, and g = 1—p, 
m be the number of abnormal cells counted, 
n be the total number of cells counted, 


at . 
c= = i be the estimate of p. 





Clearly n may have any positive value exceeding m—1. Let w,, be the probability that 
exactly n cells are counted before m »’»normal cells are observed. Two things are necessary 
and sufficient if +his is to be the case. ‘he first n — 1 cells must include m— 1 abnormals, and 


. —T "— n—1 
the nth cell must be abnormal. The prol.abilities of ti.28e two events are ( :) rr. 


ae m 
and p. Hence w, = jee 1)? mg”"-™. This is the coefficient of ¢” in (; i) ‘ 


The mean value of z is Z = (m—1) Wy, 
n=m 2 — 1 
he 2 
= MaAn—m 
ane ‘a 3)? ig 
or ifn—m=r p™ > — ‘\e 
r=0 
= p™(1—q)*™ 
= p- 


Thus z is an unbiased estimate of p. 
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The modal value of x is the value given by the value of n which makes w, maximal. 
es -1 
Set! . Hence w, exceeds w,_, and w,,,, if n lies between > a ek 
Wn n—m—1 


Pp 

2 

Therefore the modal value of x ranges between p and p+ _F—. Thus the distribution 
r—l—p 

of x is slightly asymmetrical, as is obvious from the fact that since n can have any -value 

exceeding m — 1, the range of possible values of x is from 0 to 1. If we used m/n as our estimate 


of p a bias would be introduced. In fact, the mean value of m/n can ve shown to be 








MyM a —#\)-—m . m!r!qr 
mpg {" (1-—t)-“dt, or p> sae 
On the other hand, the mean value of n/m is easily shown to be p-'. Thus if we were in the 
habit of expressing frequencies in such a form as ‘1 in 20’ rather than ‘5 %’ or ‘0-05’, the 
method of counting up to a fixed number m of abnormals, and dividing this into the total 
number n, would give us an unbiased estimate, while the method of counting the number of 
abnormals in a sample of fixed size n and dividing n by m would not do so. For, of course, 


the reciprocal of an unbiased estimate of a parameter is not an unbiased estimate of the 
reciprocal of the parameter. 








ts —1)2 
Similarly = > si ae 


ce) pa n—1~ 
(m—1)prgm ¥: (5) 4 





tm \m—2}n—1 
; a2 (n-2 
= (m—1) m pas ( )e-tat 
ali nn m—2 
qa 
= (m— prgtm {ema — ama (1) 


When m is small this can be integrated directly, e.g. for m = 4, 
< 2 2p? los: 
a. 2 (2- 2p-"2- 8?) 


For larger values of m it is better to expand the integral in an infinite series. Put u = 











q(i—t)’ 
so that ¢ = 2 , dt = ae ; 
ptqu (p+qu)? 
_ 1 ym-2 dy 
Then x? = (m—1): a) 
( )P 0 ptqu 


1 
= (m=1)p* | wm-81—g(l urdu 


ro) 1 
= (m—1)p* So {we —uy du 
r=0 0 
© (m—2)!r'iqr 
snian—-tiat 2 Se 
(e—P" & taee—i)! 
© (m—1)!r! > ipl £ 
= 7? uetiptintiom= 2 oP r 
, 2%, (m+r—1yl! . r=0 r ' 





q 2!q? 3!q8 
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Hence the variance of z is 


ee 3!¢° 
o? = 22 p? = a [teat miners}: (2) 





If in this equation we insert the estimated values of p and q, namely, 2-2 and oe 





na — + 
we find Rr re 
_ m(n—m) 2(n —m) (n— 3m 
Thuso? = ont is sufficiently accurate for all purposes and the classical value 0? = “ie m) 


sufficiently for most. These values can also be taken as the variances of p when m and n are 
given, though of course the exact values will depend on the prior probability distribution 
of p. It must be noted that if the variance is calculated from the estimated value of p it is 


approximately ~~; Pa —~ not —— Fr "4 , 80 long as p is small. Since o approximates to p J 1 , the 


standard error of pisa alate constant fraction of the value of p when p is small, provided 
m is kept constant throughout a series of determinations. The higher moments can be 
obtained by a similar method, but as they involve multiple integrals, the series expansions 
are somewhat complicated. 

Suppose now that the population sampled consists of several classes, and that a count is 
made until the number of the smallest class is m. Then the remaining n—m consist of the 
other classes, and if their frequencies are q,, 72, 7s, etc., the sum being q, the expected numbers 


a - , etc. Hence the following rule may be laid down: 





‘If a sample is counted until the number of members of one class is m,, those of the others 
being then mg, mg, etc., and the total n, then the estimate of the frequency of the first class 





, a , etc. The variances of these frequencies are approxi- 
mately m,(n—m,)n-.’ 

A formal proof presents no difficulty. 

Similar problems have of course been discussed, and there is a close analogy with problems 
concerning the duration of play when gambling (cf. Fieller, 1931). Here we have to find the 
probability that a player starting with a stake of £m, and with a probability p of winning 
£1 and a probability 1 — p of losing it per game, will be ruined after n games. This could be 
generalized to cover the case where the probability of a win was r, of a draw q, and of a loss p; 
and the problem here considered would be the degenerate case where r = 0. 

There is also an analogy with one of the methods in use for estimating the frequency of 
abnormal conditions in human families. These methods have been reviewed by Haldane 
(1938). If we find m abnormals in a family of n, the expected frequency (on certain assump- 


tions) in a very large family produced by the same parents is not ~ but = . For families 


of this type containing no abnormals would not be recorded, and therefore it is reasonable, 
as Weinberg (1927) pointed out, to take one of the abnormals as being merely a guarantee 
that the family is of a type including abnormals, and to base our frequency estimate on his 
or her sibs. 


I have discussed this problem in some detail because I believe that if it is realized that 
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frequencies can be estimated just as accurately by counting up to a certain number of the 
rarest type as by counting a certain total, haematologists and others will be saved a good 
deal of needless effort. I have to thank Drs Kestelmann and Spurway for suggestions. 


SUMMARY 


If the frequency of an attribute is estimated by counting a sample until m members possessing 
m 





-. . : : . —1., , : 
this attribute are found, and if the total number in the sample is n, then 7 san unbiased 

: : g : : m(n—m) 
estimate of the frequency, and its variance is very approximately n(n—1)" 
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THE MATHEMATICS OF A POPULATION COMPOSED OF & 
STATIONARY STRATA EACH RECRUITED FROM THE STRATUM 
BELOW AND SUPPORTED AT THE LOWEST LEVEL BY A 
UNIFORM ANNUAL NUMBER OF ENTRANTS 


By H. L. SEAL 


Consider a population subject to specified stochastic decremental forces operating at each 
attained age and supported by a uniform annual number of entrants distributed between 
ages a and # according to a known probability law. It is assumed that the population has 
existed for at least w — a years where w represents the age, not necessarily finite, at which the 
expected number of individuals in the population becomes identically zero. The individuals 
of the, now stationary, population are subdivided into k strata by titular or other distinc- 
tions conferred without reference to characteristics influencing the incidence of the decre- 
mental forces. The numbers in each of the strata are determined by the action of a stochastic 
selective force operating on the individuals in stratum g (g = 1,2,...,k—1), with an in- 
tensity which depends only on the length of time the individual has spent in that stratum; 
if selected the individual moves up into stratum (g+1). It is assumed that the total decre- 
mental and selective forces acting at each age and duration, respectively, are invariant with 
the passage of time. ; 

Write 1, for the total decremental force applicable stochastically throughout the popula- 
tion at exact age x; this function is supposed continuous and subject to 0<e<y,. If ,p, 
equals the probability that an individual now aged z survives as a member of the population 
for at least ¢ years then, 


t Ria, 
Pz = exp| — moredé | = ZS say (0<t<w—z), 


=z 
where A, > 0 may be chosen arbitrarily. 

Let AY dg, (g = 1, 2, ...,&) denote the expected number of entrants at any given moment 
into stratum g at exact age £; in particular, AB is assumed to be a normalized function of 
bounded variation with Af} = 0 (0<é<a;&>). Write p for the probability that an 
individual, if alive, remains unselected for at least ¢ years after entry into stratum g, 


(g = 1,2,...,4; p = 1); p is thus bounded and monotonically decreasing, and at points 
of discontinuity it is to be defined by 


PP = {po + vo}. 


Writing q = 1-—p, g is a normalized function of bounded variaticn in (0, w). 


It is intended to derive expressions in terms of A,, Af? and pf” (g = 1, 2, ..., &), supposed 
known, for: 


(i) AW, dx (g = 1,2,...,k), the expected number of entrants at any given moment into 
stratum g at exact age z, 


(ii) A (g = 1,2,...,%), the expected total number of individuals in stratum g at any 
moment of time. 














a2aooaorre 
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The expected number of entrants at age x into grade (g + 1) is defined to be the aggregate 
of all selectees for this grade deriving from entrants into grade g at all ages below x; hence, 
formally, 


AGT? = fia Oe ePx— pq? (g = 1,2,...,k—1), 


; dey» =) 
i.e. > = RS cae é) (1) 


zx z—£ 


The integral on the right-hand side of (1) exists when g = 1 and, after defining A?) suitably 
at points of discontinuity, A{?, is a normalized function of bounded variation in (0,w) (Widder, 


1941, Ch. 11). By induction A® (g = 3, 4, ..., &), are normalized functions of bounded variation 
in (0, @). 


The expected total of individuals in strata g + 1, 9+ 2, ..., k,is composed of all the survivors 
from entrants at various ages into grade (g + 1); thus 


at: NgiPde | ‘peat 
r= 2 0 


way» 7 
Se a Mae [ A, dé 


-|° A, dé { Agi” dx by Dirichlet’s integral theorem, 


k w 3 Pe 
or > Am = { eds | (€-z)d—2_.. (2) 
0 0 az 


r=9g+1 


Now write, s = o +17, 


[Joe e-“dgf (r = 1,2,...,k), 
FAs) = 
Lf etal (r = 0), 
0 


then (Widder, 1941) if the integrals defining f,(s) all converge for o = o,>0, the rule of 
multiplication of Laplace transforms results in 


(3) 


o AgtD 
[rewe™® = Type) = 940) say (9 = 0,1,2,..k—1) 


A, 


and, inverting, (c >>), 
Agt 1 (ett? esx 
oe ae (4) 


x To 277% e-iT 
In a similar manner, if 


Qy+» 7 [ig-mae (g = 0,1, 2,...,4-—1) 


[rewaape? = Th fle) [eat = >44(0) 
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d Qt = mil. wae mS )ds (c>a) 
oat =. 2m r 8? Pols as 
k 1 fe ct+io est 
so that 3 Mn fn rede} abule) ds. (5) 


Example (i). Let (r = 1, 2,...,k—1), 





1 dp? {c, a,<t<b,, 
~ pf) ~ \0 t<a,,t>b,, 


dt 
r) 
then adj = dpi? 


, os 





c,e-r!-4) a, <t<b,, 
~ \0 t<a,, t>b,, 
and thus (r = 2,3, ...,k—1), 


- 


by 
fs) = [; edi”? = c ven | es+eN dt = 


ay 





c 
“ (e-84r— Re-%r), 
C, 


where R, = e~¢rr—a), 
It is further assumed that 
Ap 1 t=a, 
A ={0 t#a, 


and since in this case relation (1) does not hold for g = 1, f,(s) must be replaced by unity and 
fils) given by 


fils) = [ie eae 





where AB, = AD) Pa = (a<a<w), 
2 _ >a 
ie. = 
A, dx 
so that 


1 (say a R e781). 


fils) zs [rea - , — a ses {* eda? = 
0 . 


! 8+, Cy 
C,8e- " 9 

Thus $,(8) = —_——— Tl (e*% — R,e-*r), where C, = J] ¢, 

II (s+¢,)"~" 

r=] 

-sa @ 
= phew, 
II (¢ +¢,) 
r=1 


where W,(v) assumes in succession all the possible sum-combinations of g of the quantities 
a,, b, (r = 1,2,...,g), without using an a and a 6 with the same suffix. The coefficient 7, is 


equal to the product of those — R,’s corresponding to the 6,’s used in the formation of W,(v) 





ee 











wD! 





ee 


— —— 
a ~— - 
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Hence 


+ E ihE—a—W fv) 
aif... arteledds = $2 5 T,et-<-wone (” SSE __y 


2m J c—iw (c+ it) tl (c, +e +it) 








= se ST, 1 eé—-a- sre | ethE-2-Wyr) > @, dt (cy=0,¢,#c;; r,j = 1,2 


nay r=0C, +o+tt 
G3 F Qe = WI) (E>a+W,(r)); 
where re 9, = { Wtes— a) (r = 0, 1,2, ...,g). 
(A similar type of expression may be found for A{gt”/A,.) 
And thus, by (5), 
s AM = [2 tee: 1, > Q, mie Co 


r=g+1 


=C, x % x Q,ecdet# I) NO wy yy, 


fe- a—W,(v) 


where N® w 7 kv) a J ‘ e—ChatW fr)+) A.W, )+t at. 


The above relation holds provided c,#c,; (r,j = 1, 2, ...,g); if, on the other hand, 


C, = Cg = Cz = ... = Cy, it. may be shown that 


g ctl 


: w—a—W,(v) , 
> AM = — > Tf LN ww z= (r—-)! i. hi e7t1 Res W fot at 4 


r=g+1 v=1 


The case where some of the c’s are equal presents no particular difficulties. 


The results obtained have a close analogy with those appropriate tw a problem in the 


theory of radioactive transformations (Bateman, 1910). 
Example (ii). Let 
a<t<f, 
t=a,t= £, 
10 O<t<a,t>Af, 


Aa 
r 


and (r = 1,2,...,4-—1) 





_ apf _ = _ (res e-al—e(t—¢,)r-t t>c, (a,b,>0), 
dt i 





Then Sols) = e-** —e-* 


and fs) = 
and thus (g = 0,1, 2,...,4—1) 


,(8) = TL f(s) = 


e-8r 
ee (r = 1,2,...,k—1), 
e8Cy? _ e—8Cgt 


(s+a)3o ° 
where 


C=a+ Fo, B= yo, and Ch=C8-a+f (Obs. C2= a, Ch= 8, By =). 
r=1 r=1 


lo t<c, (a+¢,>A). 
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= 
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Hence 














1 [etim ess ds 1 fetim exE—Cy% _ exXE—Cy") 
Sail uo PO = Se ie He+a)%o 
efé —Cx)c o eiKé —C,) efé —Co)ec bed etl€ —C;,") d 
= : : - ; : t 
20 ecm (a+c+it)®s 27 |" yexmp (a+c+it)¥o 
— 
But ere rm ettz Pe ra| ev yBo-(x—y)dy x>O0, 
27) _.(c+it*(a+c+it)o a 
x<0, 


(see, e.g. § 13-8 of Bochner, 1932) and 
o z—C,° oC, € 

[ Aede | eae (x — Co—y) dy = i) Agroedé e-av yBo- HE — y) dy 
J, 0 0 0 


wo—C,° w—C,° 
, = | eA yBg-1 ay| (E-y) Ag,004§ 
y 


0 
° —_ w—Zz 
Writing ;= | tA,,,dt 
0 


(cp. Steffensen, 1934), there results from (5) 


1 — a 1 w—Cy sa 
y Ar) = ——_ —ay yBy-150 day ——_— ~ay 4) By-130 
ax TB eran Sysco UY I(B,) i; eran ye y+ 07 4Y- 


I am very grateful to Dr Stefan Vajda for a number of useful suggestions in connexion 
with the preceding analysis. 
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MOMENTS OF r AND ,? FOR A FOURFOLD TABLE IN THE 
ABSENCE OF ASSOCIATION 


By J. B. 8. HALDANE, F.RB.S. 


Fourfold contingency tables are in constant use, and it is of interest to calculate the moments 
of x* derived from them when one or more of the expectations are small, even though the 
exact method of Fisher (1936) and the table of Fisher & Yates (1938) render this less important 
than would once have been the case. Further, Kendall’s (1942) use of the product-moment 
correlation of a fourfold table in connexion with rank correlation has given it a new interest. 
As the third moment of x? and the sixth moment of r are readily derivable as special cases 


from the formulae of Haldane (1945), it seems worth while to give them, along with the 
third and fifth moments of r. 


Consider the table 




















a b L 

c d l 

M m 8 
_ (ad—bc)*? 8 _ ad—be _ a 
rele" ae 


If L and / or M and m are samples of the same population falling into two classes, then 
7 = 0, and x? = 1 approximately. The exact values of the first three moments of x* are readily 
derived from equations (1) of Haldane (1945), putting 
S? S S? 2 
TD? = im? * = int Mm’ 
Let L1Mm = A?, Ll+ Mm = p, (L—1)(M—m) = v. Then 


a=2, k= 


8 
Pie: sieekeseas 
x | a 


— S8%{3(S+6)A?—6S%u+ 88+ 1)] 
‘ A2(S— 1) (S—2)(S—3) : 








: S° =72 C2 a 2 2 
v= WS) (S—B (S—3) (GS-4 (Gs) HS + 86S + 120) A*— 1082(13S + 60) A2a 
+ 1208*yu? — 5.8?(5S? + 878 + 60) A? — 3085(S+ 3) 4+ SS +1) (S?+ 15S —4)], 








Ke = [ly = WSL SLD Ga US— N7 (S" + 108 — 12) AP 6S + 84S + 1)], 


Ss 1 i 
Ps = XSL TLS H GLH Gasp US= (St + 1S" + 225"— 3088 + 240) A8 


Ke => 





— 8(S — 1) 814.8? + 79S — 120) A2u + 120842 
— 2(S— 1)-1 $3(1483 + 193S?— 518 — 120) A2— 3085S +3) nu 


+ S5(S +1) (S24 158 —4)]. (1) 
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The odd moments of r can be calculated as follows. To obtain the mean value of an odd 
power of (ad — bc) we use the operator A. This permutes the indices of a and d, and also those 
of 6 and c, without change of sign. It also permutes the indices (normal and factorial) of a 
and d with those of b and c, the sign being changed. No term is repeated if obtained more 
than once, and all are added together. Thus: 


A(a*d?) = a*d? + a*d* — b'c? — b*c*, 
A(a®bd?) = a°bd? + a®bd* + acd? + a*cd? — ab3c? — b8c?d — ab*c* — b*c*d. 
Now the mean value of a product of factorials of a, 6, c, d is 
Ela@bOcrd®] = Ela'=+9be—-oelr-oge—2)] 
= Lia+Ayrt) Met nm 6+) Sath+r+4), 
where a =a(a— 1) (a—2) ...(a—a+1), and so on (Haldane, 1940). It follows that 
EA{a™bOceMd®) = 0, if a+fP=y+é or a+y =f+4+4, 
and EA{ad®)} = S2+ LY — LOY) (M@m® — Mm), 
EA[a@bd®) = SAe+8+0[{ Le+H]@ — LOYe+D} {M@mE+D — MED 
+ {L@O+D — [8+0]@} {Mem ® — MOme+3}), 
EA{a™bd™) = — EAf{a+%d@), 
Thus (ad—be)® = A(a%d° — 3a%cd*) 
= A[a®d® — 3a™bed® + 3a%d® — 3a™bed + ad + 9a™d® + 3a@d + ad). 
So E{(ad —bc)?] = EA{ad + 3a@d}, the other terms vanishing, 
= SL — LI®)| [Mm — Mm) + 3S8S-@ LE — LI) (Mm — Mm) 
= (S—1)-A?2p. 
In general, a* must be expanded in factorials as a polynomial whose coefficients are the 
initial differences of the powers of integers, divided by the appropriate factorials. 
(ad —bc)® = A(a'd5 — 5a*bcd* + 10a°b*c?d3) 
= Afad® — 5abcd® + 10ab2e2qs) 
+ 10{ad® — Zabed® + ZaMHBA2q2 + aDp2%eq} 
+ 5{5a9d® + 20a%d® — Tabed® — 34a™bed® 4 JaMb®%c2d + 6aMb@%cd® 
+ 18a™H@B%l2q2} 
+ 5{3ad® + 50a%d® — abcd — 36a™bed® — 12a™b@c2%d + 2a™b@cd} 
+ {ad + 150a%d® + 625a%d® — 20a™bed — 165a™bed® + 30abcd} 
+ 5{2ad + 75a%d® — 3a™bed} + 25{ad + Gad} + 15ad]. 
So E{(ad—bc)*] = EA[10{a®d® + a®d® + aMbd®} + {ad + 130a%d®} 
+ 10{ad + 36ad®} + 25ad + 15a%d] 
= 108% { LOL — LAO {MOm® — Mm; + {LOS — LAV} {MOm® — Mm} 
+ {LOW — LQ} {MOm® — M@m4} + {LO — LO {MOm® — Mm) 
+ S8({ LZ — LI} {Mm — Mm} + 130{LO1® — LE} {M®m® — M@m}) 
+ 10S-9[{ Ll — LI} {M@m — Mm} + 36{ LOl® — L219}{M9m® — M@m}) 
+ 258-O[L1— LI) [Mm — Mm] + 15{LL— LI) (Mm — Mm) 
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= 10S-A2p[A2— (S— 1) 4+ (S— 1)?] (S — 5)? (S— 6)? 
+ S~)2y[134A2 — 2(8? + 60S — 55) w+ S4— 10S? + 175S?— 3608 + 230] (S—5)? 
+ 10S-®A2y[37A2 — (S? + 30S — 25) w+ S*— 128% + 948? — 2048 + 157] 
+ 25S—®A2p(S — 3)? + 1582p 
(S— 1)-“®A2r[.2(5S8 + 12) A? — 12827 + S3(S + 5)]. 
Hence the moments of r are 


My = 0, 


fe: = S-1° 
v 
= ($= 1)(8— 2)’ 
_ 3(8 +6) A?— 6S2%z + S3(S + 1) 

















Me W(S=1)(S=2)(S=3)” 
sc yf 2(5S8 + 12) A? — 12S°n + §3(8 + 5)] 
Bs = ~)8(§— 1) (S—2)(S—3)(S—4) ”’ 
[ig = [5(3S? + 868 + 120) A*— 10S2(13.8 + 60) A2z + 120S*n2— 5S*(5S2+ 87S + 60) A2 
— 30,85(S +3) w+ SS +1) (S24 158 —4)] +A4(S— 1) (S—2)(S—3)(S—4)(S—5), 
ne 2 
S68 — O14" — 68%. + 548+ 1) 
“a = 329 — 1) (S—2)(S—3) (S—4)’ 
—_1]9 2 
| — 12824 + S3(S+ 5) 
K, = 





A3(S — 1) (S—2)(S—3)(S—4) 


The odd moments vanish if Z = 1, or M =m, i.e. if the samples are equal, or the class 
frequencies 50%. «x, is negative if ZL and 1, or M and m, are nearly equal, but becomes 
positive if both a class frequency and a sample are small. 

These expressions may be used for accurate tests of the significance of observed values 


of y? or r, or, which is more difficult by the methods mentioned, above, the significance.of a 
series of values of these parameters. 
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THE USE OF ,? AS A TEST OF HOMOGENEITY IN A (nx2)-FOLD 
TABLE WHEN EXPECTATIONS ARE SMALL 





By J. B. 8S. HALDANE, F.RS. 


The x? test has proved so useful that any extension of its field of application is likely to be 
of value. Its use is at present deprecated when the expectation in any group is small. This 
may occur either because samples are small, or because one of the classes expected in each 
sample is rare. Fisher (1941) recommends 5 as a lower limit for expectations. It will be shown 
that the test may still be applied even when the expectation is less than unity. Where the 
number of classes is large, we may still sometimes use the ordinary tables. Otherwise the 
value of P must be calculated in each particular case. Haldane (1937) showed how this 
could be done when ,? is used as a test of goodness of fit. But, particularly in genetical work, 
it has found its greatest use as a test of homogeneity. The formulae required for a (m x n)-fold 
table are very cumbrous. Those for a (n x 2)-fold table are developed below. The value of 
the variance has already been given (Haldane, 1940). That of the third moment is now 
calculated. 

Consider a set of n samples of s,, 89, ...,8;, ...,8, individuals. In the ith sample let there 
be a, individuals of class X and 6, of class Y. Let Za; = A, 2b; = B, Xs; = S, so that the 
table is 


Bu 








_— 


So 





by he 


a, eee a; eee a, 
a as 





8, 8 ot 8; web ty 8 














Let Zsji=R,, Zej2=R,, S8*/4B=k. 
Then it can readily be shown that x? = s(1 - iB) , where x = 2’ (“**). Ii also follows 
i 


at once from the argument of Haldane (1940, p. 347) that in a homogeneous population the | 





expectation of ab is S+2+f A@BYs%+M, where y™ = y(y—1)(y—2)...(y—a+1l). ™ 
Hence 5; — S-+®) AB(S—n), 
x? = S-® 4*B% §?— 2(n + 2) S+n(n+ 10)—6R, + kS-* 
x {—(n? + 2n— 2) S+n(n—2)+ R, S(S+ 1}. 


The proof is similar to that which follows for 23, in which multiple summations are taken 
over all unequal sets of values of i, 7, and i, j, k. 


He 


a? = Lsz*azb} + 32287 20th} s; a,b; + CLILs7 a,b, 87 1a,b, 8, 1a,,b, 
= Ls7[aP bP +.3aP bP + 3aP bP + ab; + 9a bP + a, bP + 3a; + 3a, bY +.a;b,] 
+ 32287 * 87 [aPbPa,b, + aPb,a,b;+a,bPa,b;+a,b,a,b,] 
+ 622287 'a,b,871a,b,8;1a,b,. 
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Hence 23 = SABO L573 + 3 TT sz 245728 + CLIZ8z 19% 57149, 1 52) 
+3S-9[A@B® + AMBO) [T9739 + TLs7 2835159) 
+ S+9[A@B + 9A@B® + ABO) Fez + 38-OAMB® TL 672 3 352 5M 
+ 38-8 A@B + AB) S97 3s9 + s7@AB Zaz 35 
= S+®AB(AB— 8 +1) (AB—28 + 4) [2(s3 — 1583 + 858, — 225 + 274872 — 12087?) 
+ 322 (s?— 6s; + 11 — 6871)(s;— 1) + 6Z22(s;—1)(s;—1)(8,—1)] 
+3894 B(AB— S +1) (S—4)[2(s?— 10s, + 35 — 50872 
+ 2487 *) + TE(s,—3+ 2871) (s;—1)] 
+ S-“®AB(7AB + S*— 128 + 13) 2(s;— 6+ llsz4— 687?) 
+38-%A B(AB— 8 +1) Z2(1—s71) (s;—1) 
+38-® A B(S — 2) 5(1 — 3873+ 2872) + S-®ABL(sz— sz). 
But 2(s}—3s?+3s;—1)+322(s?—2s;+ 1) (8; — 1) + 6222(8;—1) (8; — 1) (s,—1) 
= [2(8,— 1) = (S—n}, 
2 (sj — 28; + 1) + 2Z2(s;— 1) (s;—1) = [2(s,;— 1)? = (S—n)?*, 
2(8;— 2 +871) + 2Z2(1 — 873) (8;— 1) = Z(s,— 1) (1-371) = (S—n) (n—R)). 
So 2° = SA Bl A®B*— (3S—5) AB+2(S—1)(S—2)] 
x [(S—n)§ — 12(S —n)? + 18(S —n) (n— R,) + 40S — 17€n + 256R, — 120R,] 
+3S-®AB(AB-— 8 +1)[(S—n)*?—(S—n) (n—R,) — 78 + 32n—49R, + 24K) 
+ S+®A B(7AB + S*— 128 + 13) (S—6n+ 11R, —6R,) + SA B(3n —8R, + 5R,) 
= S486 4% B35 S4{ 83 — 3(n + 4) S* + (3n? + 42n + 40) S—n(n> + 30n + 176)} 
— kS*{3(n? + 2n — 2) S?— (3n* + 51n? + 72n — 80) S + 5n(n* + 6n — 40)} 
— 2k*{(n? + 9n? + 14n — 12) S?—3(n— 2) (n+ 5) S— 2n(n?— 4)} 
+ 32R,{6S* — kS*(S— 1) (S + 10) +%(S—1)2(S+4)} 
— R,{284*(9.S — 128) — kS*(3,83 — 43,82— 168S — 120) — 2k2(11S* + 238? + 108 —4)} 
— R, S{120S% — 30k82(.S + 3) + 2.8 + 1) (S?+ 158 — 4)}}. 


= 28% S'zt\ — 38Z 3S%x* S%3 

= 2 ei wk oe —_ es ee ee ae ee 

= a= (1 AB ae): % = (1 AB* AB 7B): 
Hence xy? = S(S—1)-1(n—1), 


2 

xt = S(S—1)~*[(n? — 1) S?+ 6(2n— t) S— 6R, S*— k{ —(n?+ 2n—2)8 
+n(n—2)+ R, S(S + 1}, 

x° = S(S—1)-[ S2{(n — 1) (n+ 1) (n + 3) S2+ 2(30n? + 69n — 43) S + 120(3n— 1)} 
— kS{(3n3 + 21n2 + 24n — 26) S?— (Sn ~ 57n2 — 266n + 120) S— 60n(n— 2)} 
+ 2k*{(n3 + 9n2 + 14n— 12) S*—3n(n— 2) (n + 5) S + 2n(n?— 4)} 
— 3nR,{6S* — kS*(S — 1) (S + 10) + (S— 1)? (S+4)} 
— R,{2S%(47S + 180) — kS*(19.8? + 201.8 + 180) + 2k°(118* + 238? + 108 — 8)} 
+ Ry S{120S% — 30k82(S + 3) + k2(S + 1) (S?+ 158 — 4)}]. (1) 
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These equations are rarely the most convenient. It is usually better to subtract the classical 
values, and write 


x? = n—-14+(S—1)7*(n- 1), 
xt = n?—14 (S—1)-[(k— 6) S{R, S(S + 1) — (n? + 2n— 2) S+-n(n—2)} 
+ 6R, S?—(5n? + 12n — 11) S + 6(n?—1)], 
= (n—1)(n+1) (+3) +(S—1)-[S5{(5— k) (3n3 + 21n? + 24n — 26) 8 
+ (3n—1) (S + 120)}—(n—1) (n + 1) (n +3) (85S — 2258? + 2748 — 120) 
+ kS*{(5n3 — 57n? — 266n + 120) S— 60n(n —2)} 
+ 2k28{(n? + 9n2 + 14n— 12) S?—3n(n—2) (n+ 5) 8+ 2n(n2—4)} 
— 3nR, S{6S* — kS*(S— 1) (S+10)+k(S—1)2(8+4)} 
— R, S{289(47S + 180) — kS2(19.8? + 2108 + 180) + 2k2(1183 + 2382+ 108 — 8)} 
+ R, S*{120S? — 30k8(S + 3) + k2(S + 1) (S?+ 158 —4)}]. (2 
( 


| 


) 
To calculate the moments about the mean, we write 
x=n-l+a, xt=n®?-14+8, x = (n—1)(n4+1)(n+3)+y, 
whence Kk, = fg = 2(n—1)—2(n—1)a+ fP—2?, 
Ks = fs = 8(n—1)+3(n—1)(n— 3) a—3(n— 1) (B— 2a?) + y— 3a 8 + 2a. 3) 

The full algebraical expressions for these moments are rather cumbrous. However, by 

expanding equations (2) in descending powers of S, we have 
Ky = fy = n—14+(n-—1)S+..., 
Ke = fla = 2(n—1)+(k—6) R,—[k(n? + 2n — 2) — 2(2n? + 8n—7)] S2+..., 
Ks = fg = 8(n—1) + 2(11k— 56) R, + (k*— 30k + 120) R, 
— 2(3n — 2) [(3n + 8) kK— 4(3n + 11)] S-*+.... (4) 

In these equations the comparatively small terms involving R, and R, are omitted in the 

coefficients of S. The easiest forms for computation are 
Ka = 2(n—1)+ (k—6) R, —[(k —4) (nm? + 2n — 2) — 2(4n + 1)] S, 
ks = 8(n—1) + 2[11(k—5) — 1] R, + [(k—15)?— 105] R, 
— 2(3n — 2) [(k— 4) (3n + 8) — 12] S—. 
These equations are sufficiently accurate for most practical applications. 

The terms in equations (4) not involving negative powers of S are identical with the 
moments of x* for n — 1 degrees of freedom, used as a test of goodness of fit, derivable from 
Haldane’s (1937) equations (3), if the expectations of a, and 6, are ps, and qs,, and pg = k-. 
The effect of using x* as a test of homogeneity rather than goodness of fit is to reduce the 
variance and the positive skewness, the reduction being considerable if k is large. There can 
also be little doubt, by analogy with the simpler case, that the fourth cumulant of y* as a 
test of homogeneity is 

Ky = 48(n— 1) + 96(4k— 19) R, + 16(7k? — 125k + 420) R, 
+ (k3 — 126k? + 1680k— 5040) R, + O(S-). 

It will be seen that when k—4 is small, that is to say, the two types are almost equally 
frequent, both variance and skewness are reduced if samples are small. But when one class 
is fairly rare, so that k is large, they may be considerably increased. 
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In any particular case, one of three things may happen. It may be clear from equations 
(4) that the value of x? is not significantly above that expected on a basis of homogeneity. 
This will certainly be so if x? —(n— 1) is less than x}, and in this case there is no need to cal- 
culate x;. It may be clear that the corrections to the moments are unimportant, in which 
case the ordinary tables may be used. This is especially frequent for large values of n, where 
the distribution of x? tends to normality. Or finally it is necessary to take x, into considera- 
tion. If so we calculate x, and x, from equations (2), and make the transformation (Haldane, 

K}K 


1938) E- 2 ee hs (5) 
ote fn hx? 
where 4 = 1—-—= 


3,2 Then € is almost normally distributed with mean zero and variance 
2 





unity. The probability that x? should exceed its observed value is the probability that & 
should do so. 


MoMENTs OF x? WHEN ALL SAMPLES ARE EQUAL 
If every s, is equal, then SR, = n?, S?R, = n*. So equations (1) become 
x? = (n—1) 8(S—1)>, 
xt = (n—1) S(S—1)-®[(n— 1) S?-6(n— 1) S—2k(S—n)], 
x® = (n—1)(S—1)-©[83{(n + 1) (n + 3) S?— 2(9n? + 26n — 43) S + 120(n—1)} 
—- 2kS*(S —n) {(n + 13) S— 60(n — 1)} — 4k2(.S — n) {(n — 6) S? + 9nS — 4n}}], 


py = (n= 1) 8(8—n)(8-1)- (= 8), 





8 
fg = 2(n—1)(S—n) (S— 1) [484(S — 1)-*£.8? — (4n — 5) S— 2n} 
— kS3(S—1)-{(n — 8) S— 17n + 10} — 4k{(n — 6) S? + 9nS — 4n}}. (6) 


These equations, and those of the following paper, are a useful check on the accuracy of 
equations (1). They show that when the numbers in the samples are equal, or nearly so, 
the rarity of one class will always reduce j, and will reduce j, unless n is less than 8. However, 
Hg is always positive. 

A NUMERICAL APPLICATION 


Spurway (1945) investigated the frequency of crossing-over between a number of sex-linked 
genes in Drosophila subobscura, and used x as a test of the homogeneity of different families. 
In most cases no expectation was less than 5, and the classical method was used. However, 
for 9 pairs of genes the expectations were often very small, usually because crossing-over 
was infrequent. Her results are tabulated in Table 1. In the first column C denotes that the 
genes were in coupling, R that they were in repulsion, » is the number of families, S the total 
flies in them, A the number of flies showing crossing-over. Table 2 shows the actual data for 
the genes bg and ct, tabulated in the last row. 

K, was calculated from formulae (2) and (3), so the values are exact. It would have been 
quite sufficient to use formula (4). For example, in the case of wi and y this gives 234-9 
instead of 235-8, and in that of bgct, 17-570 in place of 17-507. The values of x, are con- 
siderably greater than the classical value 2(n — 1) in some cases, slightly less in others. x, is 
calculated from formula (4) throughout. It may also be increased up to about 12 times, or 
slightly decreased. The accurate formula is not required, since, when n is large, a com- 
paratively large change in x; does not greatly affect P. 
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The values of P are calculated from equation (5) except when P is nearly 0-5, when the 
distribution is taken as normal. It will be seen that none of the values of x? are significantly 
larger than their expectation, nor is their total. If no correction had been made, but the 


Table 1. Data on crossing-over in Drosophila subobscura. Explanation in text 











Genes n b 8 A R, R, Ky Ke Ks P 
wiyC 63 83-245 | 1342 18 | 4-7266 | 0-56584 62-047 | 235-838 | 5915 0-10 
cvlC 25 24-563 776 35 | 1-4129 | 0-12783 24-031 56-034 454-7 0-47 
cwlR 66 64-662 | 1656 88 | 4-2961 | 0-49543 65-040 | 149-245 | 1106 0-51 
lsnC 39 47-259 | 1167 | 386 | 2-2898 | 0-28976 38-033 72-162 266-1 0-14 
len R 22 20-410 694 | 283 | 1-2712 | 0-12315 21-030 39-764 143-3 0-54 
cp en C 40 39-959 | 2359 | 113 | 1-1336 | 0-06469 39-017 83-467 499-8 0-46 
epsn R | 103 | 121-970 | 5617 | 300 | 2-9840 | 0-17773 | 102-018 | 214-916 | 1220 0-09 
sn bg C 10 6-523 273 45 | 0-8752 | 0-16451 9-033 17-993 83-43 | 0-25 
bg ct C 10 9-231 273 54 | 08752 | 0-16451 9-033 17-507 74-85 | 0-48 












































Table 2. Numbers of male flies, s,, and cross-overs, a,, between the loci of bg and ct, in ten 
families. Xs>1 = 0-8752, Xs? = 0-1645. Smallest expectation of a, is 0-593 





8, 23 68 28 12 28 7 7 38 3 59 273 
a, 2 18 3 2 7 0 1 10 1 10 54 











usual distribution of x* had been employed, the values of P for wi y, cp sn R, and the total, 
would have been 0-035, 0-06 and 0-04 respectively. Thus the data would have been judged 
significantly heterogeneous. 


SUMMARY 


x? can be used as a test of homogeneity, even when expectations are less than unity, by the 


use of the formulae here given. In many cases the approximate formulae (4) are quite 
sufficient. 


I have to thank Dr Spurway for the use of her numerical data, and for help with the 
computations. 
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THE TREATMENT OF TIES IN RANKING PROBLEMS 
By M. G. KENDALL 


1. When a number of objects are presented for ranking by an observer there sometimes 
arise cases in which he is unable to express a preference in regard to certain of them and 
‘ranks them equal’ or regards them as ‘tying’. The effect may arise either because the objects 
really are indistinguishable, so far as the quality under consideration is concerned, or 
because the observer is unable to discern such differences as exist. Ties of this character are 
by no means uncommon—and indeed may be more the rule than the exception in some 
classes of work—and it is desirable to have a systematic method of dealing with them. In 
this paper I consider the effect of ties on coefficients of rank correlation, the estimation of 
rankings and the measurement of concordance in judges. 


RANK CORRELATIONS 


2. The method of allocating ranking numbers to tied individuals in general use is to 
average the ranks which they cover. For instance, if the observer ties the third and fourth 
members each is allotted the number 3}, and if he ties the second to the seventh inclusive, 
each is allotted the number }(2+3+4+5+6+7) = 4}. This is known as the mid-rank 
method and is the only one I shall consider. In fact I have seen only two other courses 
mentioned: 

(a) ‘Student’ (1921) refers to a suggestion by Karl Pearson, as an alternative to mid- 
ranks, that the ties should all be ranked as if they were the highest member of the tie. 

This is subject to the obvious disadvantages that it gives different results if one ranks 
from the other end of the scale and that it destroys the useful property that the mean rank 
of the whole series shall be }(n+1). So far as I know it has never been used in problems 
involving the calculation of ranking coefficients. 

(b) According to Woodbury (1940), DuBois (1939), in a paper which I have been unable to 
consult, has suggested allotting the ties an equal rank but proposes to determine it so that 
the sum of squares of the ranks shall be that of an untied ranking, namely, of the first n 
integers, 4n(n + 1) (2n + 1). This is rather troublesome, and it is not clear to me what advan- 
tages it possesses over the mid-rank method. 


3. The effect of ties on the calculation of Spearman’s p was considered by ‘Student’ 
(1921). p may be regarded as the product-moment correlation between two variates given 
by the two rankings. Since the variance of a ranking of n is given by ;},(n*— 1), p is given by 


12 Si Mat WH Hnt a) 


ns — 





p= 


where X;, and Y, represent the two rankings. This is easily reduced to the simpler and more 
familiar form 
i- 62(d?) 


ne—n’ 


where d; = X; aie Y,. (3) 





(2) 
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Pursuing this analogy with the product-moment correlation ‘Student’ shows that, on the 
mid-rank method, the effect of a tie of t consecutive members is to lower the variance of the 


ranking by = (# —t). This is additive for any number of sets of ties in either the X- or the 
Y-variate, and if we write Tx = 7:2(#—1), (4) 


the summation being over the ties of the X-variate, and 7, for the similar sum for Y, we 
find for the product-moment correlation, say pg, 


_ lvarX+var Y —var(X — Y) 








Pa= 3 (var X var Y) 
_ __a(n?-n)—(Tx+Ty)-2@) (5) 
V{h(n? —n) — 27x} {3 (n> — n) — 2Ty} 


$(n* —n) — (Ty + Ty) — 2(d*) 











= (6) 
Tx.-Ty ] 
n—n)—(Tx+T. Jp- 2 mah 
{3( ) ( x y)} {4(n3 —n) — (Tx +Ty)}? 
‘Student’ notes that if 7’, — 7'y is small, we have approximately 
i 2(@) ‘ 
Ps =~ Uren) (Tx Ty) i 
It is also useful to note that if 7’, and 7, are small compared with }(n?—n), we have 
62(d?) 
Ps = i- ne—n’ 


so that the correction to be applied to the ordinary formula is negligible for many practical 
purposes. 
Example 1. Consider, for instance, the two rankings of 10: 

X: 1 2 23 44 43 63 63 8 9} 93 

Y: 1 2 4 44 44 44 8 8 8 #10 
In the first ranking there are four tied pairs and hence 

Ty = 74(2° — 2) = 2. 

In the second there is one set of four ties and one of three, and hence 


Ty = 33s(60 + 24) = 7. 





We also have 2(d?) = 13. 
165 — 22 1 
Hence, in accordance with (5), 0; = aemeraa = ae 
V(161.151) 155-92 
= 0-9171. 
Calculation on the basis of (2) gives p= 1— ,48; = 0-9212. 
, ? l 
The value giyen by (7) is Ps =1- : 
165-9 
= 0-9167 


4. There is another way of looking at this problem which ‘Student’ did not mention. 
Suppose we regard any set ¢ of tied ranks as due to inability on the part of the observer to 
distinguish real differences; i.e. we assume that there does exist a set of integral ranks 
although we are ignorant of it on present evidence. Then we may ask, what is the average | 
value of p over all the ¢! possible ways of assigning integral ranks to the tied members? 
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5. If the ¢ corresponding members in the Y-ranking are held fixed, then the average co- 
variance for all ¢! arrangements of the X-members is the covariance of the fixed Y-members 
and the average of the X-members. But this latter gives the mean ranks of the tied members, 
and consequently the mean covariance of the two rankings is 


= faiy(n® —n)— 320) — (Px + Ty) (8) 


the effect of various sets of ties being additive. If now we divide by the actual variances of 
X and Y we arrive at equation (5). Thus ‘Student’s’ formula may be regarded as giving a 
mean value of the coefficient which would be obtained if the ties were replaced in all possible 


ways by the integral.ranks which they cover; always bearing in mind that we have not 
averaged the variances. 


6. A similar point of view has been adopted by Woodbury (1940) who does not seem to 
have been aware of ‘Student’s’ results; but Woodbury takes as his variance the quantity 
zs(m*— 1), that is to say, he determines the average p which would be obtained if the ties 
were replaced by appropriate integral rankings in all possible ways, the variance in each 
case being that of the first » integers. This results in py, say, where 


a he a + Ty} | (9) 





the difference from pg of equation (5) lying, of course, in the denominators in the second term 
on the right. 
Example 2. For example, in the illustration considered above Woodbury’s value would be 


6(13 + 2+ 7) 
= i-_— 
Pw= 990 


= 0-8667, 





against p, = 0-9171. 


7. The question then arises, which is the better measure of rank correlation, pg or py? 
It is useful in the first instance to consider some special cases. 


(a) Suppose that the two rankings are both completely tied, i.e. that each rank is $(n + 1). 
We clearly have 


p=1, 
indicating complete correlation. For the ‘Student’ form we have 


0 


Ps = (0 x 0)’ 


an indeterminate form which, however, may be regarded as unity as a limiting case in virtue 


of the next subsection. For Woodbury’'s form we find 


mm 0 
Pw = Tak —mn) 
= 0, 


indicating zero correlation. In short, Woodbury’s ‘correction’ has reduced p from 
1 to 0. 
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(6) Suppose that both rankings are the same, that the last member in each is ranked n 
and that the others are all tied and hence have rank 4n. Then it will be found that 


p=, 
Ps = 1, 
3 
Pum >) 


The crude form of p and ‘Student’s’ corrected form are in agreement that the correlation 
is unity. Woodbury’s form differs entirely and gives a correlation which is small for large n. 


(c) Generally, if the two rankings are identical and there are ties giving a T'-number of 
T in each 
P= ps=1, 


12T 
ne—n° 
(d) Suppose that one ranking is in the natural order 1, ..., and has no ties, so that 7'y = 0. 
If the other ranking has the last member ranked n and the others completely tied we find 


2(d?) = jyn(n— 1) (n—2), 


Pw = 1- 





BSS a* 
~ 2m+1) 2?’ 
- 3 
Ps n+l’ 
— 3 ° 
Pw atl" 


For large n, p tends to 0-5, whereas pg and py tend to zero, the latter faster than the former. 


8. It appears to me that the decision as to which of the coefficients p, or py, is preferable 
can only rest on the use to which they are to be put. 

(a) Let us suppose in the first instance that the objects have a definite ranking 1,...,” 
determined in some objective way. The purpose of correlating the order assigned by an 
observer is then to determine the observer’s accuracy, not the real ranking. ‘Student’s’ 
form of the coefficient would measure the product-moment correlation of ranks, giving 
weight to the fact that if the observer produces ties the variance of his estimates is reduced. 
Woodbury’s form would measure the average correlation of all the results obtained if the 
observer allotted to the tied groups integral ranks determined at random. For instance in 
a ranking of 8 1s be sé e's 


44444 4 4 8 


Pg = V4 = 0-577, py = 4 = 0-333. There does not seem to me to be much to choose between 
the two, but on the whole Woodbury’s form gives figures nearer to what one would expect. 
We may suppose the observer to be in a genuine state of indecision when considering the 
tied members, and the average of all the values given by guessing integral ranks at random 
seems a fair measure of his ability. ‘Student’s’ form gives higher values because he divides 
the product-moment by the actual standard deviations, and hence gives the observer credit, 
so to speak, for clustering his values in spite of the fact that he ought not to do so because 


there really is an objective order. In a case of this kind I should therefore favour Woodbury’s 
form. 











M. G. KENDALL 243 


(6) The situation is quite different if no objective order is given and we are measuring the 
concordance between two judges. In this case Woodbury’s form seems to me to give the 
wrong answer. In the case where two rankings are identical, for instance, one is entitledto 
expect that a measure of correlation should be unity—agreement could not be better. Both 
judges may be wrong, but that is not the point. We are measuring their agreement, not their 
accuracy. It has been shown above that if all members are tied Woodbury’s form would give 
a zero correlation between the two rankings, which on the face of it seems ridiculous. We are 
no longer entitled to assume an objective order, or, even if there are real differences in the 
- objects, to suppose that they fall above the threshold of the discriminatory power of the judges. 

‘Student’s’ form appears to be far better. 


9. It is, of course, undeniable that ‘Student’s’ form is more troublesome to calculate. 
This is unimportant if only two or three rankings are to be compared, but might be more 
important if there were large numbers of rankings. In such a case, however, it is more usual 
to work out a single measure of joint correlation rather than many pairs. The problem of m 
rankings for tied variates is dealt with below. 


10. I proceed to consider the appropriate method of dealing with ties in calculating the 
alternative coefficient of correlation known as 7 (Kendall, 1938). In an elegant synthesis of 
the rank correlation problem Daniels (1944) has recently pointed out that 7 may be defined as 

2 (445545) 
coe 10 
{2a}, 2b3;} a 


where a3; = — a, bi; = —b;;, 


a,; is a score allotted to each pair of ranks X;, X;as +1 ifj>iand —1 if i<j, 6,; similarly 
relating to the Y-ranking. In the ordinary ranking case, of course, 


Daniels’ form has the advantage of revealing 7 as analogous to an ordinary product-moment 
coefficient. 


11. To extend this definition to the case of tied ranks we have only to define the score 
a,; for equal ranks, and this is easily done by defining it as zero, midway between the scores 
of + l and — 1 taken when the ranks are unequal. This, it will be noticed, affects the denomin- 
ator in the definition of equation (10) as well as the numerator. 


Example 3. Let us consider the rankings of Example 1, namely, 


X: 1 Qh QE 4h 42 6} 6 8 9} 9} 
Y: 1 2 4) 4) 43 43 8 8 8 10 


Considering the first member in association with the other 9, we see that the score in both rankings 
in each case is +1, so that the total score is 9; the second and third members in the X-ranking are tied 
and this pair therefore scores 0, whatever the Y-position. The score from pairs associated with the 
second member will be found to be 7; in the Y-ranking members 3-6 are tied and therefore the only 
non-vanishing scores arise from association of the third member with the seventh, eighth, ninth and 
tenth members, score 4. The total score will be found to be 


94+74+4444443414140=33. 


The sum Za}, is found to be 41 and Yb}, is 36 and hence, writing 7, for the corresponding quantity 
to ps, we have 


33 
= ———. = 0-8589. 
7s = (41.36) 
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The value of ps, is 0-9171 but the difference need cause no concern as the coefficients do give rather 
different results, having different scales. 

The general rule for the formation of Zaj, will be clear. If there is a tie of extent ¢ we calculate 
4¢(¢— 1) and sum for all ties. If this sum is U then 


Za}, = 4n(n—1)-U. (12) 


12. If we replace any tied set by the corresponding integral ranks in any order and 
average for all the ¢! possible orders we get the same result as by replacing a,; for the tied 
members by zero; for in the ¢! arrangement each pair will occur an equal number of times in 
the order X Y and the order Y X, so that the allocation of + 1 in the first case and — 1 in the 
second is equivalent to the allocation of zero on the average. Thus we may regard our score 
for tied ranks as the mean of the values obtained by allotting integral ranks in all possible 
ways. On the analogy of Woodbury’s treatment of p we could then define 


7 = _800re_ 
% ~ yn(n—1)° 
The choice between Tg and 7,7 is precisely the same as in the case of Spearman’s p; that is 


to say we might use the latter where an objective order is known but the former where it is 
a question of measuring the concordance between judges. 


(13) 


13. For the purposes of comparison with the special cases considered in § 7 it may be worth 
while giving the corresponding values of 7’: 
(a) Both rankings completely tied: 
Tg indeterminate, Ty = 0. 
(b) Rankings identical and all tied except last member: 
2 
Tg=1, Ty = a. 
(c) Rankings identical, ties giving U-member U in each: 
2U 
Tg = i. Ty = 1-— n(n—1) ° 


(d) Oneranking equal to thenatural order 1, ...,n, the other all tied except the last member: 
2 2 
Tg = . > r= = : 


THE ESTIMATION OF A RANKING 


14. In a previous note (1942) I considered the problem of estimating the true ranking 
(or the ranking on which there was the greatest measure of agreement) for m-rankings of n 
individual exemplified by 





Object A A; A; sie A; A A, A, 
4 2 1 7 6 3 5 8 
7 2 1 6 4 5 3 8 
7 4 2 6 5 3 l 8 
Sum of ranks 18 8 4 19 15 ll 9 24 


It was shown that a reasonable estimate was to be obtained by ranking according to the 
sums of ranks, beginning with the lowest, e.g. in this example the ranking would be 


LA & & MELA he 
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This is the best in that it minimizes the sum of squares of deviations from what they would 


be if the m-rankings were identical; and it also maximizes the average p between the observed 
and the estimated rankings. 


15. The above method may also be regarded in this way: if an object is ranked r, it is 
preferred to n—r members but r— 1 members are preferred to it. Allotting as usual + 1 for 
the first type and — 1 to the second we see that the individual scores n +1 —2r in its own 


ranking. Summing over the m-rankings we see that an individual ranked X,, Y,, Z,, etc. has 
altogether a score of 


m(n+ 1)—25(X). (14) 


If then we rank the individuals according to their total scores, beginning with the highest, 
we arrive at exactly the same result as by ranking according to 2(X) beginning with the 
lowest. Thus our method arranges the objects in the order of numbers of preferences; a 
further argument in its favour. It is also easy to see that the method minimizes the sums of 
squares of deviations of preferences from what they would be if there were complete agree- 
ment. In fact, denote the estimated ranking by X,, ..., X,, and let the corresponding sums of 
preferences be £,, £5, ...,&, this being a permutation of m(n+1)—mX;,, j = 1,...,n. If the 
actual preferences are given in sum by S,, ...,S,, we have to minimize 


= (S,—£))® = D8*+DE-2(SE). (15) 


The first two terms on the right are constants and we have therefore to maximize 2(S£). 
This is clearly done by multiplying the largest S by the largest £, that is to say the largest S 
by the smallest X, and so on. In other words, we allot X, to the largest S and so on in order. 


16. To complete the story one would like to be able to prove that the method maximized 
the average 7 between observed and estimated rankings. 
Unfortunately the proposition fails, as is shown by the following example: 





A, Ay A; A 
2 3 1 4 
1 2 4 3 
1 2 4 3 
4 7 9 10 


The estimated order here would be that running from left to right across the page as written, 
and the total score between that order and the three observed orders will be found to be 
2+4+44= 10. Bui if we interchange the last two columns it becomes 0+ 6+6 = 12. Such 
a situation, however, is of rare occurrence and can only occur when there is substantial 
disagreement between judges on the two objects interchanged, in which case no ranking is 
very reliable. I do not think, therefore, that the failure of the result in extreme cases is 
important. 


17. Suppose now that some of the ranks are tied. Does the method of summing ranks 
apply unchanged to give a good estimate ? 

(a) In the first place, the method continues to give an answer which appears reasonable 
on the face of it. Moreover it may be regarded as an average result for all the ways of 
permuting the tied ranks when replaced by the appropriate integral ranks. 

(b) If the question is regarded as one of ranking according to preferences, the replacement 
of a pair of integral ranks by a tie does not affect the preferences with other members and 
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merely cancels a preference between the tied pair; and so for any set of ties. In consequence 
the method preserves the property of ranking according to the number of preferences. 

(c) If we measure the aVerage p with the estimated ranking in Woodbury’s form of cor- 
rected p the method provides a maximum average p unless the estimated ranking itself 
contains ties, in which case the result might conceivably fail, though it is unlikely to do so. 

In fact, let the estimated ranking be X,, ..., X,, and the rank of the jth object in the kth 
ranking be Y;,. We shall maximize the average p by maximizing 


V= SSX Hat VD} Lads 1) 


= F(X Hn+ HS — bmn + D), (16) 


which reduces to maximizing 2(X 8). If, however, there are ties in the estimated ranking 
our problem is to minimize something of the form 


V 
1 2” 
3(n? — 1)— > =(Tx) 





(17) 


and variations in 7'; may upset the result. This, however, is not likely to be serious unless 
there are many ties in the estimated ranking, in which case estimation of any kind is 
unreliable. 

(d) If we measure the average p with the estimated ranking in ‘Student’s’ form the result 
again may fail for (16) then becomes of the form 


Z2{X;— }(n+1)} Ap {Vj — (n+ 1)}, (18) 


where the coefficients A ,, differ from our ranking to another because they depend on differing 
variances. 

(e) Similar considerations apply to the proposition that the method minimizes the sums 
of squares of deviations from what the sums of ranks or preferences would be if all rankings 
were alike. Apart from complications due to ties in the estimated ranking, the minimal 
properties continue to obtain. 


18. To sum up, therefore, the method of estimating the ranking according to sums of 


ranks appears to give satisfactory results when ties are involved. 


Example 4. Consider the three rankings 


1 2 3 44 43 6 a 1) Dot 

1 23 23 44 44 6 63 8 9} 94 

' s 44 43 44 4) 8 8 8 10 
Sums of ranks 3 6 10 133 #133 17 22 233 263 293 
Estimated ranking 1 2 3 44 4) 6 7 8 9 10 


The sums of ranks give an estimated ranking as shown. There is one case here where the sums of 
ranks are equal and the individual ranks yielding those sums are also equal. There seems no better 
course than to tie them. Had the sums for these two been 





44 34 
4} 4} 
4} 5} 


13} 134 


an 222 ot a oe id 
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nce we might perhaps have ranked the former as 4 and the latter as 5, on the ground that the variance of 
the former is less and the group therefore ‘cluster’ better than the other. This is a new principle 
20r- ‘ deriving no support from the various minimal principles already introduced. It is the usual practice, 
self I think, to regard an estimate as better when it is based on more closely grouped observations; but 
ae hére the resulting estimate of the mean ranking is 4} so it can also be argued that the tie should 
SO. remain. On the whole this seems to me the better course.* 
kth 
TESTS OF SIGNIFICANCE FOR m RANKINGS 
19. I proceed to consider what modifications, if any, are required in significance tests 
16 when ties can appear in the rankings. Babington Smith & I (1939) have discussed the case 
(16) when the ranks are integral. The algebra required for the more extended discussion has been 
kin, given by Pitman (1938) in considering a similar problem in the analysis of variance. 
Consider an array of m rows 
a! Gece, (19) 
(17) S {- = 
ky ke k,, 
nless } 
nd is If S is the sum of square of column totals about their mean and S’ the sum of squares of all 
values about their mean, we define 
S 
esult W =— (20) 
mS 
(18) as the coefficient of concordance. In the case when the a’s, b’s, etc: are permutations of the 
oti natural numbers | to n we have 
ering g 
W =>: (21) 
zem?(n* —n) 
sums , j 
: W can vary from 0 to 1, attaining the latter value only if all rankings are identical. 
kings . g : & 
nimal 


20. Let us in the first place consider what happens to formula (21) when some of the 
integral ranks are replaced by ties. If the 7'-numbers for the various rows are T,, T,, etc., 
ms of the formula becomes 
8 


ee a. ne 
ysm*(n? —n)—m2(T') 


(22) 


This is as simple a form as we require. 
In the data of Example 4, for instance, we find 


S = 682-9, 
682-9 
sums of W = 3390.5—30 
» better 
= 0-958. 


* In my note (1942), which dealt only with integral ranks, I suggested that, where the sums of ranks are equal, 
precedence should be given to the one with the smaller variance; but I was there considering only an estimated 
ranking which itself was integral. When ties are permitted I should, as stated above, use them where sums of 
ranks are equal. 


Biometrika 33 16 
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21. Denoting by a, the rth moment of the a-row in (19), and similarly for £,, y,, etc., and 
vy a; the rth k-statistic, we have for the moments of W, from Pitman’s results, 


1 























W = EW) =~, (23) 
_ 4 a,f? 
EW-W)}= mn—1) Sa,’ (24) 
wy — 28  2%Aave , 8(m—1)(n—2) Zag fs 
Said tar m*(n— 1) Ya, 4 mint 23a, ” (25) 
aw-wy-—* Dabs 96 ZakP} 1152 Sap Bayeds 
~ m(n—1)? Sag — m(n—1)?(n+1) Day | m(n—1p> Sha, 
16(n—1) (n—2) (n—3) ZaiAy , 252(n—2) ZasBsre 196) 
m4n*(n + 1) 24a, m*n* 2a, 
In the case when the numbers are permutations of the first n integers these expressions 
reduce to 1 
W= ae (27) 
WW) = (28) 
BW - Wy =e. (29) 
B(W — Wye = 12M=—VF , 48(m—1)(m—2)(m—3)___48(m—1) (20) 


m5®(n — 1)? m?(n—1)8 m7(n— 1)? (n+1)° 


If m and n are moderately large, these expressions are approximately the same (exactly so 
for the first two) as the moments of 


1 


dF = ———. W?-'(1- W)s"dW, 31 
Bp.) ' "e 
1 
where p = t(n—-1) “*| (32) 
q = (m—1)p. 
It follows that W can be tested in Fisher’s z-distribution by writing 
2. (m—1)W 
z= flog “ae (33) 
v1 = (n-1)-2,| 


(34) 


How far does this require modification for tied ranks? 


22. For the purposes of an accurate test we can, of course, work out the first four moments 


of W from (23) to (26) in individual cases and fit an ad hoc curve; but this is a tedious process 
and some approximation is necessary. 








4) 


5) 


32) 


33) 


34) 
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The first two moments of (31) are 
Pp 





about zero) = ——, 
Hs ( ie 
Hy = es 
2 (p+)? (p+q+1)’ 


and if we identify them with the first two moments of W, ~ and ,(W), say, we find 


ye 4) 1 
Pm mip(W)’ (35) 
so that approximately W can be tested in the z-distribution with 

», = 2 2(m=1) 

te min(W)’ (36) 
Vs = (m—1)7. 

23. We have, as in (24), 
4 af. 





BW) = m*(n—1) Lax, - 
Writing A for 2a, and B for La? we have 
2 B 
so that the appropriate degrees of freedom are 
2 (n—1)(m—1) | 





71m m(l— BAY) ° 


Ve = (m—1)?,. 


(37) 


If the 7'-numbers appropriate to the various rankings are small compared with ;4,(n*—n) 
we can approximate further. In fact write N for ;(n’—n). Then 


A=mN-3X(T), 
B = mN*-2N3(T)+2(T), 
and to the first order in 2(7') 
B__mN*—2N2(T)_ 
A? m?N?—2mN2Z(T) 
_ 1, 227) (,_ 22(7)7* 
~m\ mN |\ mN | 








2 
On substitution in (27) Ry set (nx —1), (38) 


Ve = (m—1)¥. 


This, of course, is the same as (34) so that, if the number or extent of the ties is small, the 
test for untied ranks requires no modification (other than that necessary in the calculation 


of W itself). 
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24. It thus appears that we can apply the usual test unchanged unless the ties are 
numerous enough to render 2(7') not small compared with mN. If the ties are numerous we 
can work with the modified degrees of freedom given by (35), but in such a case it would 
probably be as well to verify by direct calculation that the third and fourth moments of W 
were in reasonable agreement with those given by the f-approximation implicit in the use 
of the z-test. If it happens that one or two rankings contribute the major part of 2(7') we 
may perhaps reject them on the grounds that the judges are very bad, but the rejection of 
observations has to be done with some care and only after we are satisfied that they really 
are exceptional and not merely outlying members of a continuous range. 


ADDITION OF EXTRA MEMBERS TO A RANKING 


25. There is one class of case in which I have found the coefficient 7 to have definite 
advantages over p. An example will illustrate the point best. Suppose I send out an inquiry 
to a number of firms asking for some information which they may or may not wish to dis- 
close; and suppose that the information is of a type for which one would expect that the 
non-participants might differ from the participants. By a certain date a number of replies 
have been received and it is then necessary to close the inquiry and to summarize the 
results. How far can | assume that the replies to hand are representative of the population 
to which the inquiry was addressed? Is there any evidence to suggest that those who reply 
earlier to the inquiry differ from those who reply later? 


26. To simplify the illustration suppose that I receive 15 replies in the form of a per- 
centage figure which occur in the following order: 


Order of receipt: 123 46567 8 9 10 11 12 13 14s O15 
Percentage: 165 13 12 16 25 8 9 14 17 11 18 20 10 21 19 
Rank of percentage: 8 6 5 91512 7 W 4 11 13 3 14 12 


I have chosen percentages which are all different so as not to complicate the example, but 
if ties are present there is no essential modification. 

Now if there is some relation between the order of reply and the magnitude of the per- 
centage, it ought to be shown up by the rank correlation between the order of reply and the 


order of magnitude of the percentage. The latter is shown in the last row of the example 
above and we find E(d*) = 392, 


p = 0300. 


This in fact, is barely significant, but I am not for the moment concerned with signi- 
ficance. Suppose that after we have completed this calculation two more replies arrive with 
percentages 7 and 23. We now have to calculate a revised value of p by re-numbering nearly 
all the replies and working ab initio. In practice the continual arrival of stragglers is quite 
common and to work out a new value of p each time is a great arithmetical nuisance. The 
point I wish to make is that 7 is not subject to this disability, extra values being capable of 
addition as required. 


In the above example for 15 members the value of 2(a,;6,;) is easily seen to be 


0+3+4+1-—10+94+8+4+34+24+3+4+24+1+42-1 = 27, 
80 that T = yy = 0-257. 
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If now we add a 16th member with the value 7, the contribution to Z(ab) is obtained by 
considering this new member in conjunction with the other fifteen, and is seen to be — 15. 
Similarly, a further member valued 23 adds 13. The new value of Z(ab) is thus 25 and the 
new T is given by 1 = 25, = 0-184, 

In this way a kind of running value of 7 can be ascertained without re-ranking at each stage 
as is necessary witn p. Thus 7 has a decided advantage in this class of case, namely the 
calculation of ranking coefficients for time series which may be extended in length. 
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THE PROBABILITY INTEGRAL OF THE MEAN DEVIATION 
EDITORIAL NOTE 


1. About 3 years ago a need arose to obtain the probability levels of the mean deviation 
in random samples from a normal population. The requirement was in a field of production 
quality contiol where it was customary, as well as convenient, to use the mean deviation as 
a measure of dispersion in 4 sample, rather than the standard deviation or the range between 
extreme individuals. The need was pressing, and it appeared that the quickest answer for 
practical purposes would be obtained by using the known expressions for the mean and 
variance of the M.D. and getting a measure of the departure of the distributions from normality 
by a sampling experiment with random numbers. An investigation on these lines was under- 
taken by Dr E. H. Sealy and Mr C. D. Bates of the Advisory Service on Quality Control, 
Ministry of Supply. Their results, in the form of a table of factors for control limits for sample 
sizes varying from n = 5 to 20, were‘issued in 1943. 


2. These limits, though adequate for the immediate requirement, were of course not 
exact and are now superseded in the range n = 2 to 10 by the tables printed below in 
Mr Godwin’s paper. 

3. At the time when the planning of the earlier investigation was discussed with Dr Sealy, 
I had overlooked the fact that R. C. Geary’s paper of 1936 on the distribution of the ratio 
of the M.D. to the s.D. in samples from a normal population contained formulae from which 
the higher moments of the M.D. could be derived. When this oversight was realized and 
before Mr Godwin’s work was undertaken, I had asked Dr Geary to develop from his earlier 
work, expansions for the 3rd and 4th order semi-invariants of the M.D. in terms of inverse 


powers of vy = n— 1 (where n is the sample size). It seems of interest to take this opportunity 
of putting these results on record. 


4. Dr Geary’s expansions (population standard deviation as unit). 
Mean deviation in a sample of n observations: 





2. wi 
m=— &|%—-2|. (1) 
i=1 
Expectation of m: 
m= /(20=0) ;' 
Jar (2) 


Semi-invariants of the sampling distribution of m when the population is normal: 





y \i {0-218 014 0-074170 0-057313 0-040457 0-023601 \ 
s=(—} |—_,— - —__,_ +- |; - + a 
v+1) | Yr v v4 ys vs | 
— vy \? j0-114771 _ 0-068 509 4 0-033 371 0-120 003 4 
«~\ps1) | yA sien bfietacive Sole (4) 
The frequency constants £, and f, can be obtained from 
By = AS/AZ, By = 3+A,/AG, (5) 





where A, = 0%, = ) td +./{n(n —2)}-—n+sin-! al ; (6) 
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It should be noted that Geary (1936) has taken n’ for the sample size and written n = n’— 1, 
but to be consistent with the paper which follows, I have written n for his n’ and vp for his 
n=n'—1.* 

5. The accompanying table shows the 3rd and 4th semi-invariants of m computed from 
the expansions given above and also the moment ratios £, and f,; cor n = 4 the figures for 














n f, and £, may be compared with R. A. Fisher’s (1920) values of £, = 0-297, 2, = 3-280. 
n Differences between the expansion and true values will become rapidly less as m increases. 
S 
- The sampling moments of the mean deviation 
or P Rn HINTiNGS Mi yitvenlorodw hue am walcyeriae | 
d 1 07, As A, A, By 
y } 3. oy a ere z pies 2 
' | | 
-_ | 4 0-014 32 0-001 900 0209 | 3244 ~~ | 
1, | 5 0-009 057 0-000 9785 0-230 3-194 
le | 6 0-006 244 0-000 5636 | 0-187 3-160 | 
| 8 | 0-003 483 0-000 2351 | O16 | sis | 
| 10 0-002 218 0-000 1193 0-106 3-093 | 
rt | 12 0-001 536 0-000 06864 0-088 3-076 | 
, | | 
- a 15 0-000 9800 0-000 03492 0-069 3-060 | 
| 20 0-000 5495 | 0-000 01461 9-051 - 3-045 | 
y; - ee ire Sr eae a* a ae 
io 
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* In equations (8) and (22) of Geary’s paper (pp. 296 and 300), for d read d’ = = | x, —%|/,/{n’ (n’ — 1}. 
E.S.P. 
1) 
2) 
(3) 
(4) 
\ 
(5) 
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ON THE DISTRIBUTION OF THE ESTIMATE OF MEAN DEVIATION 
OBTAINED FROM SAMPLES FROM A NORMAL POPULATION 


By H. J. GODWIN, Advisory Service on Statistical Method, Ministry of Supply 


The relative merits and demerits of the mean deviation and standard deviation as measures 
of the dispersion of a population have been discussed by Fisher (1920): though the balance 
is rather in favour of the latter, the mean deviation is widely used, especially in experimental 
work where many small samples are taken and where saving in computation is a considera- 
tion. The distribution of the estimate given by random samples from a Normal population 
has not previously been obtained, save for the special cases of sample sizes four (by Fisher, 
1920) and five (by Jones): Helmert (1876), and later Fisher, found the second moment of the 
distribution, and Geary (1936) found the third and fourth moments and showed that /, and 
8, were O(n-*) and O(n-") respectively for large sample size n. Thus the distribution may be 
approximated to by a Normal distribution, and this approximation improves as increases. 
The distribution was estimated empirically, for small sample sizes, by Sealy & Bates (1943). In 
the present paper an expression for the distribution for general n is found, suitable for calcula- 
tion by quadratures and Table 1 gives the resulting probability integral to 5 decimal places 
for n = 2to 10. From this table certain percentage points have been calculated and are given 
in Table 2. The Normal approximations to the percentage points for sample size ten are 
given for comparison with the true values: Sealy’s values were much closer, being least good 
for the extreme percentage points of the smallest sample sizes. 


Let the sample values be 2,,%2,...,%,, Where the z’s are distributed according to the 


frequency function 7 om e~+* dz (there is no loss of generality in taking the population mean 





and standard deviation to be zero and unity respectively). The numbers 1, 2, ..., can be 
assigned to the members of the sample in n! ways: we suppose that x; <2, <a... <2x,. We 


consider separately the cases when the mean falls between x, and 2,, 2, and 2g, ..., ,_, and 


22; 
«,. Suppose a, < Stet 
Then the mean deviation 


(Sega t...+%,)— (e+... +—)-=—— 2x. 
m = — 


; 2a 
i.e. MOM = Matyyy t+. +2) —(n—k) (ay +... +a). 





n 


The frequency function of m is found by evaluating 


(2m)-* |. fe-teen dx, ...d%, (1) 

over the region defined by Ly S%y_<K... <y, (2) 
| 

rRS :s SX ps1) (3) 


2m 2 
and er Stags t os tq) — (= b) (yt on ay) <M), (4) 
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lance 
ental 
dera- 
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The various functions so obtained are summed over & from 1 ‘to n—1, and the whole 


multiplied by !. The integral is evaluated by a transformation of the quadratic form 22?: 
this is most easily done in two stages. 


First put iyi — Ui = Yee 
Then (2) becomes y;>0, (5) 
(3) becomes Yt 2y2+.. -+ (k al 1) Yy_ 1 < (n— k) Yet +e FYp— 1\ (6) 
and Yr + 2yat ... + hyp <(m—k—-U) yp t---+9na J 
and (4) becomes 
n2m : n*(m + dm) 
= < (n—k)y, + 2(n—k) ygt ... + R(n—k) y, t h(n —b— 1) yg yt --- + hyp = 3 ) 
; (7) 
Now put U; = Y,+2yo+...+jy; (F< k—1), 
= (n—J— 1) Yjyat---+Yna (j > 4). 
Then (5), (6) and (7) become 
O<u,<u,<... oe. ee 
(m+dm) (8) 
and O< ty _9<Wy_g<...< Uy <~ - 
Oy, Xa» --+» Fn) ¥ 


ie, 8,...,%9) Ela B)l 
n-1 n-2 n-1 


n—1 
and Eat= nett Vex (n—jyjt VO-)y+2 DD m—Dyyy 
= j= j= 








F nm\2 *-1 uz n—2 uj m2n3 
a a ee 
= n(x + z IGin* sr) 7 IG+)) +2 a ae - 1) * dk(n—k) 


We now define a series of functions G(x); such that 


zx 2 
G(z)=1, G(x) = [exp | - set " G,_x(#) dt. (9) 


The integral (1) now appears as the product of n simple integrals and, after integration 
subject to the restrictions (8), gives 


n® 


n(2a mn nm nm 
Ener a P| ~ atte -b] ("gee “Sm 


The frequency function of m is wnt 
n® (2a m*n8 nm nm 
= D7z)—tn VA 
Fn (m m)dm n! = 5 2k'(n— k)! (2 7) /n ) ox xp| - 8k(n — “5 | G,_ ( 2 (S)¢. k— iG 9 >) am 


ni "S . (= 
~ Qn+ D7 ia—D | ‘*cyexp| - ite 6 |¢ k-1 () Ga_x- Tbe : )\am, (10) 


The calculation of the G-functions, the distribution function of m, and the percentage 
points was done under the direction of Dr H. O. Hartley, whose care and assistance I grate- 


fully acknowledge: a note by him on the method of computation appears below as an 
Appendix. 
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Although the theoretical part of the work was done privately, the computations were 
carried out, and the work prepared for publication, under the auspices of the Ministry of 
Supply (S.R. 17) from whom permission to publish has been obtained. 
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APPENDIX 


NOTE ON THE CALCULATION OF THE DISTRIBUTION OF THE ESTIMATE 
OF MEAN DEVIATION IN NORMAL SAMPLES 


By H. 0. HARTLEY, Scientific Computing Service, Ltd. 


The formula used for the computation of this distribution function is the finite series (10) 
of the above paper. This expression involves the functions G(x) which are defined by the 
recurrence formula (9). The numerical work therefore consists of: 

(a) The calculation of the G,(x) by a recurrence of numerical quadratures. 

(6) The calculation of the distribution functions f,,(m) from formula (10). 


me 
(c) The numerical quadratures f,(m)dm yielding, for each n, the probability integral 
0 
for the mean deviation m. 


(a) The essential feature of the numerical quadratures is a new method on the National 
Accounting Machine. With this metnod, of which it is hoped to publish details in due course, 
it is possible to produce a table of the integral | f(x)dx from the 4th differences of the 

Ja 
integrand f(x) in a single operation. 

The starting point of the recurrence was the table of 


1 . f* 4. 2 re, 
=~ G(x) = -- | et dt = = e—** da, 
V7 V7 J 0 NL 0 


given in T'ables of the Probability Integral, vol. 1, W.P.A., New York. 
This integral was multiplied by the ordinate 


2 2 
—_e-tz =_ r eazy" 
va v7 


which is tabulated next to it in the W.P.A. Tables. Prodycts were formed at interval 0-05 
in x and these were differenced. The function* 

(480) 2-1 2-4 G(x) 
was then obtained by the mechanical method of numerical quadrature in accordance with 
formula (9). 


This process was then repeated for r = 3,...,8; producing the functions G(x) with 
increasing constant factors and for ranges as shown below: 


r Factors of G,(x) Range 

2 480 n-!2-2 x=0(0-05) 13-5 
3 480? n-? 2-1 x=0(0-05) 19-0 
4 480° 1-2 2-2 x =0 (0-05) 25-0 
5 480! n-32-3 x=O0(01) 31-0 
6 480° n-3 2-3 x=O0(01) 36-0 
7 480° 7-2-6 x=V0(0-1) 39-0 
8 480? n-* 2-** x=0(0-1) 46-0 


* The factor (480) is necessitated by the method of quadrature. The same applies to those shown in the table 
below. 











258 Appendix 


Seven significant figures were accurate in the maximum value of G,, but this accuracy 
gradually decreased to 5 significant figures for the maximum value of G,. The ordinate 


functions 1 2 
pa? | —a9ret| 


which occur as multipliers in formula (9) were obtained by interpolation in the Tables of z 
(Table II, Tables for Statisticians and Biometricians, vol. 1). 


(6) For convenience of computation formula (10) was rewritten as follows: 
n 
f,(m) = ni 2Het DA HOD & "Oe Ge-al®) In—e-1()s (11) 


where x = nm and g,(x) = G,(x) exp | - sean |: 


Using the symmetry in k the number of terms may be halved. 
The G-functions were first converted to g-functions through multiplication by the ordinates 





c x? : 
Jan? | - r+ 5 | (c, = suitably chosen constant). 


These ordinates were obtained by interpolation in the z-tables of Table II of Tables for 
Statisticians and Biometricians, vol. 1. The x-interval was 0-05 for go, ...,g, and 0-2 for 
Js» «++» 9g. Formula (11) was then applied to obtain f,,(m) at the following x- and m-intervals: 


r 3 4 5 6 7 8 9 10 
z-interval 0-075 0-10 0-125 0-15 0-175 0-2 “0-225 0-25 
m-interval 0-05 0-05 0-05 : 0-05 0-05 0-05 0-05 0-05 


In certain cases it will be seen that the z-interval is not a tabular interval for the g-func- 
tions. In such cases Lagrangian Interpolation had to be applied. Finally the f,,(m) functions 
were subtabulated to the final m-interval of 0-01. 


m 


(c) The published tables of the probability integrals | f,,(m)dm were then obtained by 
0 


the process of mechanical quadrature on the National Machine. The tables of percentage 
points were obtained by inverse interpolation. 
Checks were as follows. Apart from the usual checks by differencing, the first and second 


m 
moments “4; and “ were calculated from the final tables of f,(m) as well as | f,(m) dm. 
0 


These were compared with check values calculated from the theoretical formulae 


»_ (2(n-1)\* 
«= (5 


2(n—1)(,, ,/ » \* 18 
“7 ee (+5) +nt(n ay. 


n27 


fy = 
Five-decimal agreement was obtained throughout. 


Grateful acknowledgement is made to Mr M. Sumner for the expert help rendered in 
the calculation of these tables. 
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Table 1. The probability integral of the mean deviation (m) in normal samples 
of n observations. 


(Population standard deviation as unit) 



































| ] 
el tee age aes See” 6 a 9 | 10 
m\ | | 
| | | 
| | | 
0-00 | 0-00000 | 0-00000 | 
01 | 01128 00019 | 0-00000 
-02 02256 00074 | -00003 | 0-00000 
03 | 03384 | -00167 | -00009 | -00001 
04 | -04511 00297 | -00022 | -00002 
| i j | 
0-05 | 0-05637 | 0-00464 | 0-00042 | 0-00004 | 0-00000 | | 
06 | -06762 00668 | -00073 | -00009 | -00001 | | 
07 | -07886 00908 | -00115 00016 | -00002 | 0-00000 | | 
08 | -09008 01184 | -00172 | -00027 00004 | -90001 | 
-09 | -10128 01496 | -00244 -00042 -00007 | 00001 
0-10 | 0-11246 | 0-01843 | 0-00333 | 0-00064 | 0-00012 | 0-00002 | 0-00000 | 
‘M1 | -12362 | 02226 | 00442 00093 | 00019 | -00004 | -00001 | 
12 | -13476 -02644 00571 | -00130 -00030 00007 | -00002 | 0-00000 
13 -14587 | -03095 00723 | -00178.| -00044 | -00011 00003 | -00001 | 
| 14 -15695 -03581 00899 | -00237 00064 -00017 00005 | -00001 | 0-00000 
| 0-15 | 0:16800 | 0-04100 | 0-01101 | 0-00310 | 0-00089 | 0-00026 | 0-00008 | 0-00002 | 0-00001 
16 -17901 04651 01329 | -00398 | -00122 | -00037 | -00012 -00004 -00001 
17 -18999 05234 01585 00503 | -00163 00053 | -00018 -00006 -00002 
18 20094 | -05849 | -01871 00626 | -00214 | -00074 | 00026 | -00009 -00003 
19 ‘21184 | -06495 | -02187 | ‘00770 | -00277 | -00101 | -00038 | -00014 -00005 
| | | | 
0-20 | 0-22270 | 0-07171 | 0-02534 | 0-00936 | 0-00354 | 0-00135 | 0-00053 | 0-00021 | 0-00008 
21 23352 -07876 02914 -01126 00445 | -00179 00073 -00030 -00012 
-22 -24430 08610 | -03327 ‘01342 | -00554 | -00232 00099 | -00042 | -00018 
‘23 | -25502 09371 | -03773 01585 00682 | -00297 00132 | -00059 | -00026 
‘24 | -26570 -10160 04254 01858 | -00830 | -00377 00173 00080 | -00037 
0-25 | 0-27633 | 0-10974 | 0-04769 | 0-02161 | 0-01002 | 0-00472 | 000225 | 0-00108 | 0-00052 
| -26 -28690 “11814 05320 | -02497 01199 00585 | -00289 | -00143 -00072 
27 -29742 -12679 | -05907 -02867 01423 -00717 -00366 00188 -00097 
-28 -30788 -13567 06528 03271 | -01677 | -00873 00459 | -00244 -00130 
-29 -31828 -14478 07186 | -03713 | -01962 ‘01052 00571 00312 00172 
0-30 | 032863 | 0-15410 | 0-07879 | 0-04192 | 0-02280 | 0-01259 | 0-00703 | 0-00395 | 0-00224 
31 | -33891 -16364 08608 | -04709 02634 01495 00858 00496 | -00289 
-32 -34913 -17337 09372 05266 ‘03025 | -01763 | -01039 | -00616 | -00368 
-33 -35928 -18330 ‘10171 05864 | -03455 | -02065 | 01248 | -00759 | -00465 
34 -36936 -19340 -11005 06503 | -03926 | -02404 | 01488 | -00928 | -00582 
0-35 | 0-37938 | 0:20367 | 0-11872 | 0-07183 | 0-04439 | 0-02783 0-01762 | 0-01124 | 0-00722 
36 | -38933 ‘21410 | -12773 | -07905 | -04996 | -03202 | -02073 | 01352 | -00887 
37 ‘39921 | -22469 | -13706 | -08670 | -05598 03666 | -02424 | -01615 | -01082 
| .38 -40901 ‘23541 | -14671 | -09476 | -06247 | -04175 | 02817 | -01915 | 01309 
-39 -41874 24626 | -15667 | 10325 | 06942 | -04731 03256 | 02257 01573 
0-40 | 0-42839 | 0-25724 | 0-16693 | 011215 | 0-07686 | 0-05338 | 0-03742 | 0-02643 | 0-01877 
‘41 | -43797 | -26832 ‘17748 | -12147 | -08478 | -05995 | -04279 | -03076 | -02224 
-42 -44747 | -27951 | -18831 | -13120 09319 | -06705 | -04869 | -03560 | -02618 
-43 45689 | -29079 | -19941 | ‘14133 | -10209 | -07468 05513 04099 | -03063 
44 46623 | -30215 | -21075 | -15186 | -11148 | -08286 | 06215 | -04693 03563 
| | | | | | 
0-45 | 0-47548 | 0-31358 | 0-22234 | 016277 | 0-12136 | 0-09160 | 0-06975 | 0-05347 | 0-04121 
-46 | -48466 | -32507 23416 | -17405 | -13172 ‘10089 | 07796 | -06063 | -04741 
‘47 | -49375 | -33661 ‘24620 | -18569 | “14255 | -11074 | -08677 | -06843 | -05425 
-48 50275 | -34820 ‘25843 | -19768 15386 | -12115 | -09621 | -07689 | -06177 
49 ‘51167 | -35982 -27086 -21001 | -16562 | 13212 | -10627 | 08602 | -06998 
0-50 | 0:52050 | 0-37146 | 0-28345 | 0-22265 | 0-17783 | 0-14364 | 0-11697 | 0-09584 | 0-07892 
| eG 
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Table 1 (cont.). The probability integral of the mean deviation (m) in normal samples 


of n observations. (Population standard deviation as unit) 








a | 5 6 v 8 9 






































10 
e neeRe: BV ME Soe A 
0-50 | 0-52050 | 0-37146 | 0-28345 | 0-22265 | 0-17783 | 0-14364 | 0-11697 | 0-09584 | 0-07892 
-51 52924 | -38313 | -29620 | -23559 | -19047 -15569 -12829 -10635 -08860 
-52 | -53790 -39479 | -30909 | -24881 -20352 | -16828 -14023 -11756 -09903 
-53 54646 | -40646 | -32211 | -26230 “21697 | -18137 -15279 -12947 -11022 
54 | -55494 | -41811 | = -33525 | ‘27604 | -23079 | = -19497 -16595 -14207 -12218 
0-55 | 0-56332 | 0-42975 | 0-34847 0-28999 | 0-24497 0-20904 | 0-17970 | 0-15536 | 0-13491 
-56 | -57162 | -44135 36178 | -30416 -25948 | -22357 | -19402 -16931 -14840 
‘57 | 57982 | -45293 | -37516 | -31851 | 27430 | -23852 | -20889 | -18392 | -16264 
-58 | -58792 -46446 | “38858 | -33302 | -28941 | -25389 | -22427 -19916 17761 
59 | -59594 | -47593 | -40204 ‘34768 | -30477 | =-26963 | -24016 21501 | -19329 
| 
0-60 | 0-60386 | 0-48735 | 0-41552 | 036245 | 032037 | 0-28572 | 0-25650 | 0-23144 | 0-20965 
‘61 | -61168 -49871 -42901 -37733 | -33617 | -30213 | -27328 -24840 -22667 
-62 “61941 -50999 | -44249 | -39228 | -35215 | -31882 | *29046 -26588 -24431 
63 | -62705 52120 | -45594 | +40729 -36827 | -33576 | 30799 -28383 -26252 
-64 -63459 -53232 -46936 -42233 -38452 | -35293 32585 | 30220 28127 
| 0-65 | 0-64203 0-54335 | 0-48273 0-43739 0-40087 | 0-37027 | 0-34398 | 0-32095 0-30050 
66 | -64938 -55428 -49603 45244 | -41727 -38777 -36235 | -34004 -32016 
| -67 -65663 ‘56511 | -50926 | -46746 | -43372 | -40537 | -38092 | 35941 | -34021 
| .68 | -66378 -57583 | -52240 -48244 | -45017 -42306 -39965 37902 -36058 
| 69 | -67084 58644 53544 | -49734 -46661 | -44078 | -41848 | -39882 | -38122 
’ | | } > | | 
| 0-70 | 0-67780 | 0-59693 | 0:54836 | 0-51217 | 0-48300 | 0-45852 | 0-43739 0-41875 | 0-40206 
71 -68467 -60729 ‘56117 | -52689 -49932 -47623 -45632 43877 -42306 
-72 -69143 -61753 ‘57384 -54148 51555 | -49388 | -47523 45882 | -44414 
-73 -69810 -62764 -58637 -55594 -53166 51143 | -49409 | -47886 | -46525 
‘74 | -70468 62762 -59874 57025 54762 52886 | -51284 | -49883 | -48634 
| 0-75 0-71116 0-64745 0-61096 0-58439 0-56342 | 0-54614 | 0-53147 0-51868 | 0-50735 
-76 *71754 -65714 -62300 -59834 -57903 -56324 | -54992 -53838 -52821 
77 -72382 “66669 -63487 -61210 -59443 -58012 -56816 -55788 -54888 
‘78 | -73001 -67609 -64655 -62564 -60961 -59677 -58616 57713 | -56930 
| “7 73610 -68534 -65804 -63897 -62454 -61316 -60388 -59609 -58943 
0-80 0-74210 | 0-69443 | 0-66934 | 0-65206  0-63922 0-62926 | 0-62130 | 0-61474 | 6-60922 
| -81 -74800 -70337 -68043 -66492 -65362 64506 -63839 -63303 | -62864 
| -82 -75381 -71215 69131 67752 -66773 -66054 65512 65092 | -64763 
83 75952 -72078 -70198 -68986 68154 67567 67147 66840 | -66617 
-84 -76514 -72924 -71243 -70193 -69503 “69045 -68742 -68544 -68422 
| 0-85 0-77067 0-73754 0-72266 071373 | 0-70821 0-70486 | 0-70295 | 0-70201 | 0-70175 
| +86 ‘77610 -74568 -73266 -72525 | = -72105 ‘71888 | -71805 “71810 | 71875 
| ‘87 | -78144 -75366 -74244 -73649 -73355 -73252 -73270 -73368 | 73519 
88 -78669 -76147 -75199 ‘74744 | -74571 | -74575 -74689 -74874 | -75105 
-89 -79184 -76912 -76132 “75810 -75752 *75857 -76061 -76327 | -76632 
0-90 | 0-79691 | 0:77660 | 0-77040 | 0-76847 | 0-76898 | 0-77097 | 0-77386 | 0-77727 | 0-78039 
91 -80188 -78392 ‘77926 | -77854 ‘78008 | -78296 -78663 -79072 -79505 
-92 -80677 -79107 ‘78789 | -78832 ‘79082 | -79453 “79891 -80363 -80851 
-93 *81156 -79806 79628 | -79780 | -80120 | -80567 -81071 81599 | = -82135 
94 | -81627 “80489 "80444 “80698 “81122 -81640 -82202 *82780 | 83359 
0-95 | 082089 0-81155 | 0-81237 | 081587 | 0-82089 | 0-82671 | 083286 | 083907 | 084522 
| +96 82542 -81805 ‘82007 | -82447 83021 -83660 -84321 -84981 85626 
| 97 -82987 “82439 ‘82754 | -83277 “83917 “84608 “85310 “86001 -86671 
-98 83423 -83057 ‘83478 | -84079 84778 85515 “86252 -86970 87659 
-99 83851 83659 84180 “84852 85605 86382 87149 ‘87887 | -88591 
hem 0-84270 | 0-84245 | 0-84860 | 0-85597 | 0-86398 | 0-87210 | 0-88001 | 088755 | 0-89468 
‘ | j | 
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Table 1 (cont.). The probability integral of the mean deviation (m) in normal samples 

















0-89612 
89910 | 
-90200 | 
-90484 
90761 


0-91031 
-91296 
91553 
91805 
*92051 





0-92200 
92524 
9275! 
92973 
93190 


0-93401 
93606 
-936,07 
94002 
94191 


0-94376 
94556 
94731 
-94902 
-95067 


0-98229 





98335 
95538 
95686 
95830 


0-95970 
96105 
96237 








of n observations. (Population standard deviation as unit) 





0-91275 
91635 
91984 
92321 
92647 


0-92961 
93264 
93556 
93838 
94109 








0-87935 
88487 
89020 
89533 
-90027 


0-90502 
-90959 
-91398 
91820 
-92224 


0-92613 
92985 
93341 
93682 
-84009 


0-94321 
4620 
4905 
5177 

195437 


0-95685 
95921 
-06146- 
96360 
96564 


; 096758 


-96942 
‘97117 
97283 
97441 


0-97591 
97733 
97867 
97995 
“98115 


0-98229 
98337 
98439 
98535 
98626 


0-98712 
98793 
98869 
98941 
99009 


0-99073 











5 6 7 8 9 10 
0-85597 | 0-86398 | 0-87210 | 0-88001 | 0-88755 | 0-89468 
86315 87158 87999 “88809 89574 -90292 
-87005 “87885 88751 89575 90347 -91065 
-87668 88581 89465 -90299 91074 -91788 
*88305 *89245 90144 -90984 91756 92464 
0-88916 | 0-89878 | 0-90788 | 0-91629 | 0-92397 | 0-93095 
89502 90482 91398 92237 92997 93681. 
90063 91057 91976 92809 93557 94227 
-90600 91604 92521 93346 94081 94733 
‘91114 92123 93037 93850 94568 95202 
0-91605 | 0-92616 | 0-93523 | 0-94322 | 0-95022 | 0-95635 
92074 93084 93980 94764 95444 96034 
92521 93527 94411 95176 95835 96403 
92947 93947 94815 95561 -96198 -96741 
93353 94344 95195 -95920 96533 97052 
0-93740 | 0-94718 | 0-95551 | 0-96253 | 0-96843 | 0-97337 
94108 -95072 95885 96564 -97128 -97598 
94458 95406 -96197 96851 97391 97836 
94790 95720 -96488 97118 97633 98054 
95105 96016 -96760 97365 97855 98252 
0-95404 | 0-96294 | -0-97014 | 0-97594 | 0-98058 | 0-98432 
95687 96555 -97250 97805 98245 98595 
95955 -96800 97470 97999 98415 98743 
96209 97030 -97675 98178 98571 98877 
96449 97246 97865 98343 98713 -98998 
0-96676 | 0:97448 | 0-98041 | 0-98495 | 0-98842 | 0-99107 
96890 97637 98204 98634 98958 -99206 
97092 97813 98355 -98761 99066 99294 
-97283 97978 98495 98878 99.162 99373 
97462. 98132 -98624 98985 99250 99445 
0-97631 | 0-98275 | 0-98743 | 0-99082 | 0-99329 | 0-99508 
97790 98408 98852 99172 99401 -99566 
-97940 98533 98953 99253 99465 -99616 
98081 98648 99046 99327 99523 99662 
98213 98755 99132 99394 99576 -99702 
0-98337 | 0-98855 | 0-99210 | 0-99455 | 0-99623 | 0-99738 
98453 98947 99282 99510 99665 99770 
98562 99033 -99348 99560 99703 99798 
98664 99112 99409 -$9606 99736 99823 
-98759 99186 99464 99647 99767 99845 
0-98849 | 0-99254 | 0-99515 | 0-99684 | 0-99794 | 0-99865 
98932 99316 -99561 99717 -99818 99882 
-99010 99374 99603 99748 99839 99897 
-99083 99427 99641 99775 99858 99910 
99151 99477 -99676 99799 99875 -99922 
0-99215 | 0-99522 | 0-99708 | 0-99821 | 0-99890 | 0-99932 
99274 99564 99737 99841 4 99941 - 
-99329 99602 99763 -99859 -99915 -99949 
99380 -99637 99786 99874 -99926 99956 
99427 -99669 99808 99889 99935 99962 
0-99472 | 0-99699 | 0-99828 | 0-99901 | 0:99943 | 0-99967 


























Table 1 (cont.). The probability integral of the mean deviation (m) in normal samples 
of n observations. (Population standard deviation as unit) 





















































a“ 2 3 4 5 6 7 
m 
1-50 | 0-96611 | 0-98385 | 0-99073 | 099472 | 0-99699 | 0-99828 
51 -96728 98471 99133 -99513 -99726 99845 
52 -96841 98552 99190 99551 99751 99861 
53 -96952 -98630 99243 -99586 99774 99875 
54 -97059 -98704 99293 -99619 99795 99888 
1-55 | 0-97162 | 0-98774 | 0-99340 | 0-99650 | 0-99814 | 0-99900 
-56 -97263 -98841 -99384 -99678 -99831 -99911 
“57 -97360 -98905 99425 -99704 99847 -99920 
-58 97455 -98966 -99464 -99728 99861 -99929 
-59 -97546 -99023 -99500 99750 99875 -99936 
1-60 | 0-97635 | 0-99078 | 0-99534 | 0-99771 | 0-99887 | 0-99943 
-61 97721 -99130 99566 -99790 -99898 -99949 
-62 -97804 -99179 99596 -99808 -99908 99955 
63 97884 99226 99624 99824 -99917 -99960 
-64 -97962 99270 99650 | 99839 -99925 99964 
1-65 | 0-98038 | 0-99312 | 0-99675 | 0-99852 | 0-99932 | 0-99968 
-66 -98110 99352 99698 99865 -99939 -99972 
-67 98181 -99390 99719 99877 99945 99975 
-68 -98249 99425 99739 99887 -99951 -99978 
69 98315 99459 99758 99897 -99956 -99980 
1-70 | 0-98379 | 0-99491 | 0-99775 | 0-99906 | 0-99960 | 0-99983 
71 98441 99521 99791 99915 -99964 -99985 
72 -98500 99550 -99806 99922 -99968 -99986 
73 98558 99577 -99820 -99929 99971 -99988 
74 98614 -99603 99834 -99936 -99974 -99989 
1:75 | 0-98667 | 0-99627 | 0:99846 | 0-99941 | 0-99977 | 0-99991 
76 98719 -99650 99857 99947 -99979 -99992 
77 -98769 -99671 99868 99952 -99982 -99993 
78 98817 -99691 99878 99956 -99983 -99993 
79 98864 99710 99887 -99960 99985 -99994 
1-80 | 0-98909 | 0-99729 | 6-99895 | 0-99964 | 0-99987 | 0-99995 
81 98952 -99746 -99903 -99967 -99988 99995 
82 98994 99762 99911 -99970 -99990 -99996 
83 -99035 99777 99917 -99973 -99991 -99996 
84 -99074 99791 99924 99976 -99992 -99997 
1-85 | 0-99111 | 0-99804 | 0-99930 | 0-99978 | 0-99993 | 0-99997 
86 99147 99817 99935 -99980 -99993 -99997 
87 99182 99829 -99940 -99982 -99994 -99998 
88 99216 99840 99945 99984 -99995 -99998 
89 99248 99850 -99949 -99986 -99995 -99998 
1:90 | 0-99279 | 0-99860 | 0-99953 | 0-99987 | 0-99996 | 0-99998 
91 998 99869 99957 -99988 -99996 -99999 
-92 -99338 -99878 -99960 -99989 -99997 -99999 
93 99366 99886 -99963 -99990 -99997 -99999 
94 -99392 99894 -99966 -99991 -99998 -99999 
1-95 | 0-99418 | 0-99901 | 0-99969 | 0-99992 | 0-99998 | 0-99999 
96 99443 -99907 -99972 -99993 -99998 -99999 
97 -99466 99914 99974 -99993 -99998 | 1-00000 
98 99489 -99920 99976 -99994 -99998 
99 99511 99925 -99978 99994 99999 
2-00 | 0:99532 | 0-99930 | 0-99980 | 0-99995 | 0-99999 
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of n observations. (Population standard deviation as unit) 








(cont.). The probability integral of the mean deviatian (m) in normal samples 





























5 6 m 2 3 
0-99532 | 0-99930 | 0-99980 | 0-99995 | 0-99999 2-50 | 0-99959 | 0-99999 
99552 -99935 -99981 -99995 -99999 51 -99961 -99999 
99572 -99940 -99983 -99996 -99999 -52 -99964 -99999 
99591 -99944 -99984 -99996 -99999 -53 -99965 -99999 
99609 -99948 -99986 -99997 -99999 -54 -99967 -99999 
0-99626 | 0-99951 | 0-99987 | 0-99997 | 0-99999 2-55 | 0-99969 | 0-99999 
-99642 -99955 -99988 -99998 -99999 -56 -99971 -99999 
-99658 -99958 -99989 -99998 | 1-00000 -57 -99972 -99999 
-99673 -99961 -99990 -99998 -58 -99974 -99999 
99688 -99964 -99991 -99998 59 -99975 -99999 
0-99702 | 0-99966 | 0-99992 | 0-99998 2-60 | 0-99976 | 0-99999 
-99716 99969 -99992 -99999 61 -99978 -99999 
99728 -99971 -99993 -99999 -62 -99979 | 1-00000 
99741 -99973 -99994 -99999 -63 -99980 
99753 99975 -99994 -99999 -64 -99981 
0-99764 | 0-9997 0-99995 | 0-99999 2-65 | 0-99982 
99775 -99979 -99995 -99999 -66 -99983 
99785 ‘99980 -99996 -99999 -67 -99984 
99795 -99982 -99996 -99999 -68 -99985 
99805 -99983 -99996 -99999 -69 -99986 
0-99814 | 0-99984 | 0-99997 | 0-99999 2-70 | 0-99987 
-99822 -99985 -99997 | 1-00000 71 -99987 
99831 -99987 -99997 72 -99988 
99839 -99988 -99997 73 -99989 
99846 -99989 -99998 74 -99989 
0-99854 | 0-99989 | 0-99998 2:75 | 0-99990 
-99861 -99990 99998 | 16 -99991 
99867 -99991 99998 | 77 -99991 
99874 -99992 -99998 | 78 -99992 
99880 -99992 -99999 | -79 -99992 
0-99886 | 0-99993 | 0-99999 | 2:80 | 0-99992 
-99891 -99993 99999 ‘81 -99993 
-99897 99994 | -99999 -82 -99993 
99902 -99994 -99999 | -83 -99994 
-99906 99995 -99999 ‘84 -99994 
0-99911 | 0-99995 | 0-99999 | 2:85 | 0-99994 
99916 -99996 -99999 -86 -99995 
99920 99996 | -99999' -87 -99995 
99924 99996 | -99999 -88 -99995 
99928 99997 | -99999 -89 -99996 
0-99931 | 0-99997 | 1-00000 | 2:90 | 0-99996 
99935 ‘99997 | 91 -99996 
99938 99997 | 92 -99996 
99941 99998 | -93 -99997 
99944 99998 | 94 -99997 
0-99947 | 0-99998 | 2:95 | 0-99997 
99950 | -99998 96 | -99997 
99952 -99998 97 -99997 
99955 -99998 | | -98 -99998 
99957 -99999 -99 -99998 
0-99959 | 0-99999 | | | 3-00 | 0-99998* 





i 


























* 0-99999 is reached for m = 3-07; 1-00000 is reached for m = 3-23. 
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Table 2. Percentage points of the probability integral of the mean deviation (m), 
with the population standard deviation as unit 





(a) Lower percentage points 






































) oe ers 
Size of | 
sample 0-1% 0-5 % 10% 25% 50% 10-0 % 
n 
|— | | 
2 0-001 0004 | 0009 | 0022 | 0044 0-089 
3 0-022 0052 | 0073 | O116 0-166 0-238 
4 0-066 0114 | 0145 | 0-199 0-254 0-328 
5 0-112 0-170 | 0203 | 0-260 0-315 0-386 
6 0-153 0-215 | 0-250 0-306 0-360 0-428 
7 0-190 0-252 | 0-287 0-342 0-394 0-459 
8 0-220 0283 | O318 | O872 | 0-422 0-484 
9 0-247 | 0310 | 0344 | 0396 | 0-445 | 0.504 
10 0-271 0-333 | 0366 #=%| O417 | 0-464 0-521 
Normal approximation: 
10 | 0-171 | 0269 | 0-316 | 0-386 {| 0-445 | 0-514 
| | 
(6) Upper percentage points 
| | | ] | 
Size of | | | 
sample 10-0 % 50% 25% $| 10% | OS % 01% 
| | | | 
n | | 
| enon 4 | Kelme: 
| | | 
.. 7 1-163 | 1-386 | 1585 =| 1-82] 1-985 | 2-327 
3 may | 1276 =| 1-417 | 1-586 | 1-703 | 1-949 
4 1-089 | 1224 4 86| #1844) «| 1-489 1590 =| 1-806 
Bon 4 1-069 =| 1-187 1-292 1-419 | 1507 | 1-693 
iy 2 1-052 1158 | 1-253 1-366 1-445 | 1-613 
7 | 10388 | 1135 | 1-222 1325 | 1397 | 1/550 
S| 1.026 | 1-116 1-196 1-292 | 1:358 | 1-499 
9 1016 | 1-100 1175 | 1-264 | 1-326 | 1-457 
10 1007 | 1-086 1-156 1-240 1-299 | 1-492 
Normal approximation : 
1-000 1-069 1:245 | 1-342 


i} 
== 


| 1128 | 1-198 


| 
i 
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BOOK REVIEW 


The Advanced Theory of Statistics. Vol.1. By Maurice G. Kenpati. London: Charles 
Griffin and Co. Ltd., 1943. Pp. 457. Price 42s. 


It is difficult to review the present volume without knowing precisely how the author will deal with 
the topics reserved for its promised successor. For although Mr Kendall expresses the hope that this 
first instalment can profitably be read before the publication of the second, it is clear; nevertheless, that 
the two parts will be complementary and that full justice can only be done to this first part after the 
two volumes have been considered together. 

Mr Kendall defines his objective as the provision of a systematic treatment of statistical theory as it 
exists at the present time. The work is encyclopaedic and will receive little criticism on the grounds of 
what is omitted. It is not an elementary book, the various topics being all carried to an advanced stage 
and at times requiring of the reader considerable mathematical powers. As with advanced theoretical 
work in most of the sciences, the practical problems which originally suggested the discussions have 
often receded well into the background. This is not mentioned as a criticism, for the sole value of 
scientific work does not necessarily lie in it beitig of immediate or even of ultimate practical importance. 
It is, however, proper to point out that Mr Kendall is here, in the main, content to present us with a 
picture of statistical theory as he finds it. The more controversial job of assessing the value of the 
different parts of the structure, whether from a purely practical or from an aesthetic viewpoint, he 
leaves to the reader. Within these self-imposed limitations he has scored a notable success. 

The first six chapters deal in some detail with the properties of frequency distributions. Chapter 3, 
entitled ‘Moments and Cumulants’, is here particularly satisfying. It develops concisely the general 
relationships between the various families of power statistics. The simple presentation of the trans- 
formations giving moments in terms of cumulants and vice versa will be welcomed, and the listing of a 
large number of the resulting formulae will enhance the value of the book as a reference work. In 
addition, the familiar corrections for grouping, due to Sheppard, are derived, and the rather subtle 
distinctions between the conditions necessary for their application and the conditions under which the 
so-called average corrections for grouping may be applied are clearly drawn. 

Chapter 4, entitled ‘Characteristic Functions’, begins with a proof of the Inversion Theorem which 
states that the characteristic function uniquely determines the distribution function. It then discusses 
at length various theorems connected with the limits of infinite sequences of distribution functions, 
and with the so-called Problem of Moments, i.e. the problem of specifying conditions under which the 
moments determine uniquely the distribution function. As in discussions of the convergence of infinite 
series in pure mathematics, the interest in the limits of sequences of distribution functions is almost 
exclusively theoretical, and this chapter should appeal to anyone who enjoys himself in this type of 
work. The more practical question of determining approximations with the aid of only a few moments 
is considered in later chapters. 

Chepter 5 introduces the simpler distributions which are of central importance in statistical practice. 
Chapter 6 continues with a description of the Pearson system of curves and of the series developments 
of the Normal and Poisson distributions, associated with the names of Gram, Charlier and Edgeworth. 
One welcomes the fact that both the Gram-Charlier and the Edgeworth developments from the normal 
distribution are given, for although these series may be but rearrangements of each other, it is the 
order and grouping of the terms that are all important when any practical applications are intended. 

Indeed, although this point is well made, it is of such importance that, even at the risk of appearing to 
labour it, the addition of some further numerical illustrations might have been useful. 

The next five chapters deal with Theories of Probability, with Sampling and with Sampling Dis- 
tributions. On the first topic the author rightly does not dogmatize, judging that the rest of the subject 
appears to go forward in much the same way, whatever are one’s basic concepts as to the meaning of 
probability. This determination not to take sides in the present book is extended also to the problem of 
induction, where he minimizes the too violent contrasts which have in recent years been drawn between 
Bayes’s Theorem and the Principle of Maximum Likelihood. Far from being diametrically opposed, 
Mr Kendall observes that, if some account is taken of the limiting processes by which continuous 
distributions are defined, and if Bayes’s Postulate is introduced in an appropriate manner, the two 
principles have a very strong resemblance. 

The chapter devoted to the so-called ‘exact’ sampling distributions follows familiar lines. Perhaps 
of greater interest is Chapter 11, which is entitled ‘Approximations to Sampling Distributions’. Here 
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the problems associated with the distribution of the k-statistics are given 2, very full treatment. It is, 
moreover, an authoritative treatment since Mr Kendall has himself contributed so much to the elucida- 
tion of the methods which are here required. The k-statistics, which were introduced by R. A. Fisher 
in 1928, have the property that their expected values in repeated samples are the cumulants of the 
populations sampled. In general the exact distributions of the k-statistics cannot be derived, but the 
cumulants of the k-distributions can be obtained by following out certain rules which were given by 
Fisher. Mr Kendall describes these rules and gives full proofs of their validity. A large number of the 
formulae derived from applying the rules in particular cases are quoted for reference. This is a very live 
subject, and one feels that there is still some scope for development and simplification of these methods, 
for it must still be admitted that their successful application requires a high degree of virtuosity. 
Aesthetically, however, they are far in advance of the heavy algebraic manipulation demanded by the 
earlier approach. The mathematician will perhaps find most in the chapter to satisfy his sense of what 
has come to be called elegance in such studies. 

After this the reviewer found the remainder of the book dealing with the x? distribution and with the 
Theory of Correlation something of an anti-climax. Most of all at this stage one feels the need to refer 
to the promised second volume, to which has been assigned the general theory of regression analysis. 
Correlation and regression are so closely allied that even a temporary separation is distressing. 

While the above remarks may give some idea of the nature of the subjects discussed in this first 
volume, they do not do justice to the thoroughness with which the author has accomplished his task. 
It is by no means a book to skim through. Mr Kendall notes rather ruefully in his preface that statis- 
tical theory is essentially mathematical and suggests that it is not easy to keep the mathematics from 
getting on top. He says, however, that he intends his work to be one on statistics and not on statistical 
mathematics. The distinction is a fine one, and I must confess to some difficulty in appreciating 
Mr Kendall’s point here. Readers who are familiar with Dr Aitken’s small work entitled Statistical 
Mathematics may perhaps ask whether his book ought not, for the same reasons, to be termed one on 
statistics and not on statistical mathematics. However, by whatever name we speak of this subject, 
anyone who can follow a mathematical argument, and who has also at least some small experience of 
practical statistical problems, will firid plenty to reward him in a study of Mr Kendall’s book. 


B. L. WELCH 
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